Hostname: page-component-586b7cd67f-vdxz6 Total loading time: 0 Render date: 2024-11-25T22:28:09.912Z Has data issue: false hasContentIssue false

Influence of tides and tidal current on Mertz Glacier, Antarctica

Published online by Cambridge University Press:  08 September 2017

Benoît Legrésy
Affiliation:
Legos (CNES-CNRS-UPS-IRD), 18 Av. Edouard Belin, 31401 Toulouse Cedex, France E-mail: [email protected]
Anja Wendt
Affiliation:
Institut für Planetare Geodasie, Technische Universität Dresden, D-01062 Dresden, Germany
Ignazio Tabacco
Affiliation:
Dipartimento di Scienze della Terra, Sez. Geofisica, Università degli Studi di Milano, Via Cicognara 7, I-20129 Milan, Italy
Frédérique Rémy
Affiliation:
Legos (CNES-CNRS-UPS-IRD), 18 Av. Edouard Belin, 31401 Toulouse Cedex, France E-mail: [email protected]
Reinhard Dietrich
Affiliation:
Institut für Planetare Geodasie, Technische Universität Dresden, D-01062 Dresden, Germany
Rights & Permissions [Opens in a new window]

Abstract

Mertz Glacier, East Antarctica, is characterized by a 140 km long, 25 km wide floating ice tongue. In this paper, we combine a large number of remotely sensed datasets, including in situ global positioning system measurements, satellite radar altimetry, airborne radio-echo sounding and satellite synthetic aperture radar imagery and interferometry. These various datasets allow us to study the interaction of the ice tongue with the tides and currents. However, the inverse barometer effect needs to be applied to sea-level variations affecting the tongue. We find that the tide-induced currents exert a small lateral pressure on the tongue which, when integrated over the large surface of the tongue, induce a flexure of up to 2 m amplitude per day. Simple elastic modelling of the flexure confirms the observations and helps validate the boundary conditions necessary to explain different eastward and westward tongue deflections. In addition, the along-flow velocity of the tongue does vary daily from 1.9 to 6.8 md-1 depending on the tidal current. When the current pushes the tongue toward the eastern boundary of the fjord, the tongue is retarded by the drag and the velocity decreases. The accumulated stress is released, allowing the tongue to flow very rapidly when the current pushes the tongue westward. These forcing and boundary conditions on the floating ice flow are important and must be taken into account when studying glacier discharge and calving.

Type
Research Article
Copyright
Copyright © International Glaciological Society 2004

1. Introduction

The George V Coast of East Antarctica has three main glacier outflows: Mertz and Ninnis Glaciers and the Cook Ice Shelf. They have the peculiarity of producing large ice tongues several tens of km long which sometimes break up, as occurred at Ninnis Glacier in 1999/2000. These tongues regularly produce large icebergs, which ground on the shallow waters of the continental margin.

These large ice tongues were first reported by the Australian expedition of 1911−14 led by Sir D. Mawson. A Soviet expedition in 1957 resurveyed the region and showed large changes in size and shape of the Mertz and Ninnis tongues. The Deep Freeze 79 expedition helped map the bathymetry of the area more precisely (Reference Domack, Anderson, Oliver, James and JagoDomack and Anderson, 1983). The bathymetry of the area is characterized by the presence of a deep valley running parallel to the coast, the George V Basin where the water depth reaches 1000—1400 m. A large continental margin extends 150 km north of the present-day coast, with water depths of a few hundred metres. The continental margin also comprises banks of shallow water, such as the Mertz and Ninnis banks, which are preferential resting areas for large icebergs calved by the ice tongues. This cyclical release of large ice masses is found to play a major role in the sedimentation and shaping of the continental margin (Reference Barnes, Eittreim and HamptonBarnes, 1987; Reference Domack, Eittreim and HamptonDomack, 1987).

There have been several attempts to measure the ice discharge of these glaciers, including that by Reference Frezzotti, Cimbelli and FerrignoFrezzotti and others (1998), who estimated the accumulated discharge of Mertz and Ninnis Glaciers to be 62 Gt a-1 and the accumulated ice in the drainage basin to be 58 Gt a-1. In this study, we concentrate on the Mertz ice tongue (Fig. 1), which extends over 140 km from the grounding line to the front, with a width of 30 km at the front edge.

Fig. 1. ERS SAR mosaic of the Mertz Glacier area. European Space Agency single look complex data were used to make this high-resolution mosaic. Straight white lines correspond to altimeter tracks, and the straight black dashed line to the radioecho sounding (RES) profile. The frames of the SAR scenes are indicated, and section A–B corresponds to the profile used to fit the tidal vertical bending. The curved white line corresponds to the grounding-line position recovered from the interferograms, the curved black line to the coastline and grounding line given in the Antarctic Digital Database, and the straight black dashed line to the RES flight-line. The black dotted lines correspond to the boundary of the ice tongue considered as effective for the mechanical bending.

Reference Wendler, Ahlnas and LingleWendler and others (1996) estimated a flow velocity of 1km a-1 on the ice tongue by picking details on repeated Japanese Earth Resources Satellite 1 (JERS1) synthetic aperture radar (SAR) images and comparing the tongue front position with former observations. Potzsch and others (2000) derived the grounding-line position and the topography of the tongue basin using SAR interferometry, and reported good agreement between the observed and modelled tides.

The study of tidal interaction with the floating ice has a long history. Reference HoldsworthHoldsworth (1969) set the formalism of the vertical flexure of a floating ice tongue due to tidal forcing and questioned the type of interaction to use (elastic or plastic). More recently, simultaneous measurements of tide and ice flow show that there must be some link between the tide level and the ice velocity of glaciers on the floating part of the ice shelf. Reference Riedel, Nixdorf, Heinert, Eckstaller and MayerRiedel and others (1999) used the global positioning system (GPS) to study the ice flow on the grounding zone of Ekstromisen, Dronning Maud Land, Antarctica. They found a clear influence of the ocean tides on the ice movements on the floating part and at least 1 km upstream of the grounding line. They found that high/low tides imply a landward/seaward movement.

Reference DoakeDoake and others (2002) used 30 days of continuous GPS observations on the Brunt Ice Shelf and found that tides influenced the lateral movement of the ice shelf with a delay between the vertical tidal movement and the horizontal tide-linked movement. Reference Anandakrishnan, Voigt, Alley and KingAnandakrishnan and others (2003) found a very strong modulation of the flow speed of Ice Stream D due to the tide beneath the Ross Ice Shelf. Finally, Reference Bindschadler, Vornberger, King and PadmanBindschadler and others (2003a) observed a link between tides and the flow of Ice Stream B, invoking a stick-slip driven motion. They observed that the motion events did not correspond to high tides and did not relate generally to tides in a simple way. They came to the conclusion that an elastic strain storage-and-release mechanism, depending on the tide level (Reference Bindschadler, King, Alley, Anandakrishnan and PadmanBindschadler and others, 2003b), helps explain the relation between tides and flow variations. Bind- schadler and others (2003a) mention the possibility that tidal currents play a role in the floating ice movement by drag below the ice.

In the present paper, we concentrate on studying the interaction between the Mertz ice-tongue flow and the ocean tides and currents. Firstly, we describe our complete set of observations. We study the effect of an ocean current intersecting laterally with the tongue flow direction, using the formalism of an elastic tongue flexure. We then show that our observations can be explained by a lateral tidal current, and we discuss the implications for the ice-tongue flow and break-up.

2. Datasets

We used four types of remotely sensed data to estimate the ice displacements. Vertical displacements were estimated using the FES99 ocean tide model (Reference Lefevre, Lyard, leProvost and SchramaLefevre and others, 2002) and the inverse barometer effect (IBE) based onpres- sure measurements from the nearby automatic weather station (AWS) at Penguin Point (http://uwamrc.ssec.wisc.edu/aws/awsdata.html). All of the vertical displacements will be compared to this (Tide+IBE) reference. The locations of the available data are shown in Figure 1.

2.1. Radio-echo sounding (RES) data

In 2000 a survey of the George V Coast was made by the Italian Antarctic Expeditions. In particular, a profile was acquired along the Mertz tongue from the ice-shelf front up to the grounded area (dashed line in Fig. 1). The 60 MHz central frequency and 16 MHz bandwidth radar were mounted on an airplane. Laboratory tests give a 10 m error for the radar measurements. A noisy signal was obtained over the tongue, perturbed by the large crevasse patterns, and this can increase the error by up to 100 m in some places. The profile is, however, continuous along the ice tongue. Cross-checking between altimetric height and RES data, assuming hydrostatic equilibrium and a mean ice density of 0.9, gave us confidence in both datasets and revealed that, at least up to longitude 145° E, the tongue can be considered as composed of pure ice. This can easily be understood given the strong winds which flow down the coast in this area and erode the snowpack. The accumulation may be more important over the northern part of the tongue where the snow transported by wind is deposited on the tongue as the wind spreads over the flat ocean.

2.2. Radar altimetry data

The European Remote-sensing Satellite (ERS-1 and -2) radar altimeters have recorded data since 1991. Two tracks of the 35 day repeat mission cross the Mertz tongue (track numbers 145 and 417, marked in Fig. 1). Two other track segments acquired during the ERS-1 geodetic mission are also shown in Figure 1. The raw altimeter data have been processed within the OSCAR project using the retracking algorithm described by Legresy and Remy (1997) and used with the name Ice2 on the ENVISA Tradar altimetry processing chain. These data were corrected for ionospheric and tropospheric range delays. An average value was computed over a 10 km segment for each repeat measurement. The comparison for each date with the (Tide+IBE) data shows fairly good agreement (Fig. 2).

Fig. 2. Comparison of the different remote-sensing measurements of sea-level variations vs the tide (FES99 model) + IBE (computed from Penguin Point pressure gauge). All values are in m. The comparison is made with GPS data (a), ERS altimetry (b) and ERS SAR interferometry (c). m is the regression coefficient, r is the correlation and σ is the rms data dispersion.

2.3. GPS data

The main goal of the German Antarctic North Victoria Land Expedition GANOVEX VIII in 1999/2000 was the geological and geophysical investigation of the region. During this expedition, GPS measurements were carried out on the tip of the Mertz tongue (Reference Korth, Bozzo and DamaskeKorth, 2001). Due to unfavourable weather conditions, only 4 days of data, with a 1 day gap, were collected. The data were analyzed with the Bernese GPS Software version 4.2 (Reference BeutlerBeutler and others, 2001). The 85 km long baseline between the Mertz Glacier tongue front and Penguin Point (a fixed rock site) was solved pseudo-kinematically. An iterative approach was applied to account for the movement of the receiver within the 20 min time slices. For this reason, there was no need to introduce any a priori information on the horizontal and vertical movements of the ice tongue. The accuracy of the height determination was assessed by solving a static baseline in the same way. Since only the 1300 km long baseline (Penguin Point−McMurdo) was available for this assessment, the resulting standard deviation of 8.5 cm can be regarded as an upper limit. Taking into account the vicinity of the GPS observations to the magnetic pole and the maximum of solar activity resulting in ionospheric disturbances, this 8−9 cm upper limit for the accuracy seems reasonable. Comparison with the (Tide+IBE) data shows very good agreement (Fig. 2).

2.4. SAR data

We computed four interferograms using ERS-1 and -2 data acquired in 1996. Three interferograms have been shown in Potzsch and others (2000). Table 1 summarizes all the information related to the SAR scenes that we used. All interferograms are from the tandem period (this means using 24 hour exact repeats between ERS-1 and -2), and they all correspond to descending tracks with a look direction from the southeast perpendicular to the frame direction (Fig. 1). One fringe in an interferogram corresponds to a 28 mm displacement in the satellite line of sight. Figure 3 shows one of these single-difference tandem interferograms. In order to retrieve the grounding line, we formed double-difference interferograms in which the remaining fringes can be attributed to the residual vertical tidal flexure of the ice beam (Reference RignotRignot, 1996). The fringes observed at the grounding line on these double differences do correspond to such vertical tidal flexure. Fitting this flexure as described in section 3 gives the corresponding vertical displacement. According toTable 1, the comparison with the (Tide+IBE) shows fairly good agreement (Fig. 2). We will see in the following sections that the double-difference interferograms also compound other non-steady displacement. GPS, and SAR data are also useful for studying horizontal displacements.

Table 1. Dates and conditions of acquisition of the SAR data

Fig. 3. Interferogram corresponding to the May (158) ERS SAR scenes. The many parallel fringes visible on the tongue reveal the horizontal flexure of the tongue in between the two acquisitions. The grounding line is drawn in white.

3. Ice-Tongue Mechanics

Hereafter we will consider elastic bending. Reference HoldsworthHoldsworth (1969) provides the formulation that the second derivative of the bending is proportional to the bending moment. We observed horizontal bending of the tongue with both GPS and interferometric SAR (InSAR) data. Reference Holdsworth and JacobsHoldsworth (1985), while looking for physical processes affecting floating glacier tongues, mentioned one that can apply to the specific situation of the Mertz: ocean current action on the ice. Effectively, the currents exert a pressure on the floating ice; this situation is investigated in this section.

We take Reference HoldsworthHoldsworth’s (1969) formulation and apply it to the case of a pressure force acting horizontally, orthogonal to the tongue flow, and solve it on the neutral axis running along the centre of the tongue.

If y is the horizontal position perpendicular to the main direction of the tongue,

(1)

where E b is the elastic modulus of the ice, y is Poisson’s ratio, taken as 0.3, M (x) is the bending moment and Ia is the moment of inertia. We assume that the tongue has a rectangular section with variable width (W) and thickness (H).

On the segment A−B shown in Figure 1 near 145° E, we computed an ice thickness of 440 m from hydrostatic equilibrium and altimeter measurements, which show a flat, stable surface over a few kilometres. We computed the bending parameters of Reference HoldsworthHoldsworth (1969) using double-difference interferograms, according to the method described by Reference RignotRignot (1996). We computed the tide-damping factor for several parallel transects on the double-difference interfero- grams I 1 (April (430)—May (158)) and I 2 (April (430)− May (430)) (see Table 1). We chose this area because the grounding line follows a straight line, so there cannot be curvature effects on the parameter estimates, such as those described by Reference Rabus and LangRabus and Lang (2002). We did not observe any shift of the grounding-line position depending on the tide level, indicating that the bedrock slope must be large here. We thus consider the grounding line to be fixed throughout the tidal cycle. We estimated the grounding position, damping factor, amplitude of the tidal deflection and reference level by non-linear least-squares fitting. We estimated parameters for numerous parallel transects in the area and found the damping factor to be stable at 0.78 km-1which leads to a realistic value for the elastic modulus of 0.88 GPa (Reference VaughanVaughan, 1995). Large errors can arise for this parameter, both because the fitting minimum is not well defined by the procedure (Reference VaughanVaughan, 1995) and because a small error in thickness has a strong impact on the recovered E b value (H = 440 ± 10 m leads to E b = 0.88 ± 0.06).

To complete the scheme we computed I a according to

(2)

and M according to

(3)

where ρi is the ice density which we take as 900 kg m-3 and v is the component of the water current orthogonal to the tongue causing the excess pressure. The ice thickness over which the current acts is Hρi/ρw. The origin of the × axis is at the front of the tongue and oriented towards the grounding zone. An additional force should come from the drag of the current below the tongue. This additional force would be proportional to Wv 2, and it can then be taken into account by adding a drag coefficient in front of M(x) in Equation (1).

The horizontal flexure of the tongue due to the pressure exerted by the ocean current is found by integrating Equation (1) twice. We set M(x) at 0 on x = 0, which means that no lateral force is exerted on the tongue for x 0 since there is no tongue in this case. We also set y(x)=0 on the grounding line (x = Xg ), and both conditions are sufficient boundary conditions for the numerical integration. This latter boundary condition is less restrictive than the former and does not presume any ice deformation upstream of the grounding line, which might occur (the latter could be assumed as negligible in front of the flexure along the 140 km long ice tongue).We do not investigate here the effect of coupling between vertical and horizontal flexures; however, the impact on horizontal flexure would be very small. In the case of the Mertz tongue, F. Lyard (personal communication, 2003) finds that maximum tidal currents are of the order of 0.25 m s 1 during the ascending and descending tides, using the FES99 model. This current acts on the tongue between the main shore and the tongue edge, i.e. over the first 80 km of the tongue.

The pressure force from the usual tide differences is of order 1 kPa, which is much larger than the pressure force exerted by the current (about 30 Pa). The main difference between the horizontal flexure and the vertical flexure under tidal forcing is linked to the dimension which resists the bending. This dimension is the ice-section width for the horizontal flexure instead of the height for the vertical flexure, and ice tongues are usually much wider than they are thick. The other important difference is that there is no limitation to the horizontal bending, whilst vertical bending is limited by the tidal amplitude. Our ice flexure model will be limited by the poor knowledge of the current velocity as a function of distance along the tongue and by the unknown water/ice drag coefficient below the tongue. However, our results demonstrate the important effects which must apply. As a first case, we model the flexure assuming the thickness measured along the RES flight represents the whole crosssection. We assume the section width, which is mechanically effective, is the one that follows western and eastern flowlines at the grounding zone which are marked in Figure 1. Finally, we take into account the width variation when the ice spreads (from 18 to 26km) (I a = I a(x)). Figure 4 shows the results of our model for a 1m s 1 westward current applied to the first 80 km. Since alms-1 current is used in the model, changing the current to v, one must scale the result by v2 Since the thickness increases in the second part of the tongue and no more pressure is applied, the momentum does not increase significantly from the shoreline to the grounding line. Figure 4a shows the width and thickness profile used for the flexure computation. Figure 4b shows the momentum per thickness unit profile computed according to Equation (3). Figure 4c shows the bending computed for both the eastward and westward current. We obtain a 2.1 m lateral displacement of the edge of the tongue for a 1ms-1 westward current or 0.13 m for the 0.25 m s-1 maximum current obtained with FES99. Note that significant bending arises mainly from the shoreline to the edge, where the current is most effective.

Fig. 4. Modelling the horizontal elastic flexure of the tongue for a 1m s-1 lateral ocean current.

4. Observations of the Tongue Bending

4.1. Horizontal bending from InSAR observations

Looking at InSAR double differences Ii and I2 (Fig. 5), we can see that fringes corresponding to lateral displacements are present. Sixteen fringes can be counted on Ii while only three are seen on I2. According to the projection of the displacement in the satellite line of sight, we find 1.15 m horizontal bending for Ii and 0.21 m for I2. This difference can be explained by the fact that the May (track 158) interferogram corresponds to a 17 cm tide difference with the ERS-2 acquisition at a tide level of −17 cm (Table 1), the nearest to tide level 0 of all of the eight acquisitions. Tide level 0 corresponds to the largest tidal current (linked to the divergence of the tide-level time variations). Thus, the horizontal bending of the ERS-2 orbit 05523 acquisition will be the largest of the eight available scenes. The observed parallel fringes in Ii are thus in agreement with a horizontal bending following tidal currents, with values comparable to the model for a 0.74 m s lateral current. This is three times greater than the expected maximum current speed, and the tide differences of Ii may not lead to the maximum tidal- current speed difference. Either the current drag under the tongue must be at least twice the lateral pressure force, or another source of deformation must lead to this large number of fringes. The linear fringe pattern reveals that the tongue is flexing in phase, even though it is 20 km wide. It also has a linear trend (regular fringes), which implies a linear bending which is effectively modelled (Fig. 4). In Figure 4c, we put the profile of effective extracted deformation from the interferogram along a longitudinal line on the ice tongue. It compares to the kind of deflection expected by the model (westward and eastward flexure on the same graph).

Fig. 5. Double-difference interferograms showing the horizontal flexure of the tongue. One fringe corresponds to 28 mm distance variation in the satellite line of sight. (a) Double difference I i. Both vertical flexure at the grounding line and lateral flexure of the tongue are observed. (b) Double difference I 2. Mainly vertical flexure near the grounding line is present.

However, the double-difference interferogram corresponds to the difference of four epochs, each having its own flexure. In this study, we do not have actual ocean current values for each epoch. It is thus impossible to reproduce a reliable flexure model comparable to this particular line. None of the tide levels for the I2 acquisitions corresponds to very significant tidal currents. The April interferogram shows both acquisitions have similar tide levels separated by 1 day, so there must be a similar current, whereas the May interferogram has two high-tide acquisitions, which have a near-zero tidal current and thus a small current difference. It is thus logical that there are few fringes on this double difference, although a few fringes are always observed on individual interferograms which compensate each other.

4.2. Horizontal bending from GPS observations

In Figure 6 we show different representations of the GPS data acquired during 3 days of observations. Although short, this record is sufficient for several important conclusions to be drawn. Figure 6a shows the kinematic positions of the GPS beacon in along-flow and across-flow directions. According to our processing scheme, and given that we are in a polar region and the horizontal movements we find are larger than in the vertical, these positions cannot result from processing artefacts (Reference King, Coleman and NguyenKing and others, 2003). The right panel is the same with mean velocity removed. One can see there is both an across-flow movement and an along-flow movement modulation and they are linked so the along-flow movement increases when a westward shift occurs, and decreases when the tongue shifts eastward. Figure 6b shows the vertical position and velocity (or derivative of the height position) as well as the horizontal along-flow ice velocity. It clearly appears that the along-flow speeds are in agreement with the water-level variations. The ice flows at up to 7 m d 1 during ascending tides and slows down to 2 m d during descending tides.

Fig. 6. GPS measurements on the ice tongue. The measurement point is situated at the edge of the tongue (Fig. 1). (a) The horizontal position along the 4 days of observation. Abscissa is along-flow distance from the first measurement; ordinates are across flow distance (dark lines) and height (light lines).The right-side plots correspond to the same along-/across-flow position, but with the average velocity component removed. (b) The height level (light grey lines), along-flow velocity (black dots) and height change or time derivative of the height (triangles).Do Y, day of year. (c, d) The across-flow position vs the height level for ascending (c) and descending tides (d). +/− values correspond to eastward/westward respectively. Different dots shades correspond to different days of observation. (e, f) The along-flow speed vs the height level for ascending (e) and descending tides (f).

4.3. Across-flow movement

From the GPS observations, we can draw the conclusion that ascending tides correspond to a westward movement while descending tides lead to an eastward movement of the tongue. This is confirmed by the representation in Figure 6c and d. A positive derivative of the vertical position corresponds to ascending tides (Fig. 6c), and a negative derivative to descending tides (Fig. 6d). Positive across-flow displacement corresponds to eastward movements. For incoming tides, a large westward displacement occurs (with a peak around medium tide) while smaller (and more widely spread) displacement occurs eastward during descending tides.

To explain the difference between eastward and westward movement, we must go back to the horizontal bending modelling. Looking at both the interferogram from Figure 2 and the amplitude image in Figure 1, one can see that there are a few kilometres of pack ice or dead ice between the western side of the valley in which the tongue develops and the western side of the ‘active tongue’ (the one defined with the two flowlines originating from the grounding line at the back of the glacier). This dead ice must be fractured by the movements of the tongue, and the group of visible large cracks indicates that it must be very weakly resistant to the tongue pressure. By contrast, the eastern side of the tongue seems to join the coast wall, which runs straight over >20 km, and in particular the eastern flowline seems to “touch” the valley side just at the end before the shoreline. If so, the lateral movement is stopped by this wall. This is what we modelled in Figure 4 as the second negative solution. We apply the same 1ms lateral tidal current but now in the eastward direction and we add boundary conditions y = 0 and dy/dx = 0 at x = 75 km to simulate the contact of the tongue with the western side rock. In this case, the horizontal bending is more limited (five times less than in the westward direction). This new boundary condition, which may not be so strict in reality, explains very well the difference between eastward and westward movement observed with the GPS data (factor 2-3).

5. Along-Flow Movement Modulation by Tides

Although Figure 6a shows more peaky westward movements, Figure 6b shows the horizontal velocity is much larger during ascending tides. Figure 6e and f show the along-flow velocity variations during ascending and descending tides respectively. As in the across-flow bending case, the along-flow increase is large and peaky during the ascending tides while it is quite constant during descending tides. The flow-speed modulation by the tide is very important (a factor of 3). We do not observe a large impact on the interferograms since the satellite track is almost parallel to the tongue. The angle between the satellite track and the tongue flow direction is only 6°. Thus the projection of the along-flow velocity as seen by the satellite is small. Interferograms do measure the movement of the targets integrated over the 24 hour period. Thus, even if one of the two interferograms used to form Ii (Fig. 5a) compounds a longer period with high velocity, the difference in accumulated displacement may not be very important. It can be a few tens of cm. In this case, the projection in the satellite line- of-sight direction would be a few cm, i.e. a maximum of two or three fringes. At present, we do not have enough precise elements available to discriminate these fringes in I i from the lateral bending fringes.

What we observe on single- and double-difference inter- ferograms is the existence of ripples on the phase picture near 67.5° S, 144.75° E. These ripples are not always present at exactly the same position. They are in the middle of the tongue and correspond approximately to the limit of influence of the shoreline or to the east side position where the tongue is in contact with the rock as discussed above. They also correspond to the end of the large opening in the ice on the western side, and to the origin of the series of transverse crevasses, which can be seen in Figure 1. These crevasses are parallel to each other and limited by two main longitudinal fractures, which are visible up to the edge of the tongue. On the western side of the western longitudinal fracture and on the eastern side of the eastern fracture, Figure 1 shows crevasses developing at a 45° angle relative to the direction of main flow or of the longitudinal fractures. These last crevasses must come from the natural spreading of the ice, free of external stress.

The large variations of the flow speed in phase with tidal cycles may be linked to a variation in friction of the ice on the eastern rocky wall.

The next explanation is more speculative and would need more data to verify but it may well explain our observations correctly. During descending tides, the tidal current pushes the tongue toward the eastern rocks, which increases the friction so that the flow speed is reduced; stress accumulates on the ice with speeds at ~700 m a-1 instead of the ~1200 m a-1 of the natural regime. When the pressure is released, ascending tides push the tongue westward and eliminate the east side friction. The accumulated stress can be elastically released through an increased speed of up to ~2400 m a-1, much greater than the natural flow regime. These strike slip cycles can be assured within the elastic hypotheses since the displacement during the period of the ascending tide is within 1 m. This is the same order of magnitude as the flexure created either laterally by the current or vertically by the tides. The elastic assumption is thus able to explain the tongue behaviour we observe. If no elastic condition prevails, the tongue could be curved as suggested by the Reference HoldsworthHoldsworth (1982) measurements for the Erebus ice tongue following the non-symmetry between the eastern and western sides. The existence and development of crevasses indicates that the deformation must be plastic. Formulation in terms of effective thickness or width, as developed by Reference Lingle, Hughes and KollmeyerLingle and others (1981), should help conciliate the elastic behaviour with the important crevassing of the tongue.

Conclusion

A number of remotely sensed datasets, including in situ GPS measurements, satellite radar altimetry, airborne radio-echo sounding and satellite SAR imagery and interferometry, were used to study the Mertz Glacier tongue. With these data, we studied the movement of the ice tongue, both vertically and along and across the flow direction. We show that tidal currents, which operate along the George V Coast, must exert a significant pressure on the tongue to create its observed flexure. A simple model of the pressure forces and their integration helped us compute the elastic flexure resulting from the lateral pressure of tidal currents. GPS observations confirm that a boundary condition must be added on the eastern side of the fjord, where the lateral flexure of the tongue must be stopped by rigid rocks. Finally, observations of the temporal variation of the along-flow movement linked to the tidal signal suggest that this lateral pressure on the tongue changes the friction of the tongue on the eastern side of the fjord and makes the along- flow velocity vary temporally in phase with the tides.

This kind of study must be improved with a better specification of all of the forces and boundary conditions acting on the floating ice sections. Dedicated tidal height and current modelling would help improve our understanding of the action of the tidal currents on the tongue. Drag coefficients could be studied, as well as our assumption that the elastic modulus adopted is satisfactory. Further GPS observations, which would include an array of GPS receivers at several points on the Mertz Glacier tongue, would better describe the coupling between tides (or sea-level variations), boundary conditions and ice flow. They may explain the coupling observed between the tides and the ice flow in this area and on other floating ice shelves or ice streams of Antarctica. Coupled modelling of the floating ice as proposed by Reference Schmeltz, Rignot, Dupont and MacAyealSchmeltz and others (2002), if applied to Mertz Glacier including the lateral boundary and flow conditions we observe here, would help greatly to improve our knowledge of the Antarctic outlet glacier dynamics. Finally, this kind of study may help better describe the mechanics offloating ice and be a basis for studying more precisely the calving of giant icebergs.

Acknowledgements

We would like to thank R. Morrow of the Laboratoire d’Etudes en Geophysique et Océanographie Spatiale (Legos), France, and R. Coleman of the University of Tasmania, Australia, who helped improve the English and provided very helpful pre-review. Two reviewers made constructive comments which helped improve the manuscript.

References

Anandakrishnan, S., Voigt, D.E., Alley, R.B. and King, M.A.. 2003. Ice Stream D flow speed is strongly modulated by the tide beneath the Ross Ice Shelf. Geophys. Res. Lett., 30(7), 1361.(10.1029/2002GL016329.)Google Scholar
Barnes, P. W. 1987. Morphologic studies of the Wilkes Land continental shelf, Antarctica: glacial and iceberg effects. In Eittreim, S. L. and Hampton, M. A., eds. The Antarctic continental margin: geology and geophysics of of fshore Wilkes Land. Houston, TX, Circum-Pacific Concilium for Energy and Mineral Resources, 175-193. (CPCEMR Earth Science Series 5a.)Google Scholar
Beutler, G. and U others. 2001. Bernese GPS software. Version 4.2. Berne, University of Berne. Astronomical Institute.Google Scholar
Bindschadler, R. A., King, R.B. Alley, M.A., Anandakrishnan, S. and Padman, L.. 2003a. Tidally controlled stick-slip discharge of aWest Antarctic ice stream. Science, 301(5636), 1087-1089.Google Scholar
Bindschadler, R.A., Vornberger, P L., King, M. A. and Padman, L.. 2003b. Tidally-driven stick-slip motion in the mouth of Whillans Ice Stream, Antarctica. Ann. Glaciol., 36, 263-272.Google Scholar
Doake, C. S. M and 6 others. 2002. Tide-induced lateral movement of Brunt Ice Shelf, Antarctica. Geophys. Res. Lett., 29(8), 1226. (10.1029/ 2001GL014606)Google Scholar
Domack, E. W. 1987. Preliminary stratigraphy for a portion of Wilkes Land continental shelf, Antarctica: evidence from till provenance. In Eittreim, S. L. and Hampton, M. A., eds. The Antarctic continental margin: geology and geophysics of offshore Wilkes Land. Houston, TX, Circum-Pacific Concilium for Energy and Mineral Resources, 195-204. (CPCEMR Earth Science Series 5a.)Google Scholar
Domack, E.W. and Anderson, J. B.. 1983. Marine geology of the George V continental margin: combined results of the Deep Freeze 79 and the 1911-1914 Australian expedition. In Oliver, R. L., James, P. R. and Jago, J. B., eds. Antarctic earth science. Cambridge, etc., Cambridge University Press. Australian Academy of Science, 402-406.Google Scholar
Frezzotti, M., Cimbelli, A. and Ferrigno, J. G.. 1998. Ice-front change and iceberg behaviour along Oates and George V Coasts, Antarctica, 1912-96. Ann. Glaciol., 27, 643-650.Google Scholar
Holdsworth, G. 1969. Flexure of a floating ice tongue. J. Glaciol., 8(54), 385-397.Google Scholar
Holdsworth, G. 1982. Dynamics of Erebus Glacier tongue. Ann. Glaciol., 3, 131-137.Google Scholar
Holdsworth, G. 1985. Some effects of ocean currents and wave motion on the dynamics of floating glacier tongues. In Jacobs, S. S., ed. Oceanology of the Antarctic continental shelf. Washington, DC, American Geophysical Union, 253-271. (Antarctic Research Series 43.)Google Scholar
King, M.A., Coleman, R. and Nguyen, L. N.. 2003. Spurious periodic horizontal signals in sub-daily GPS position estimates. J. Geod., 77(1-2), 15-21.CrossRefGoogle Scholar
Korth, W. 2001. Realization of the geodetic and gravimetric project. In Bozzo, A. and Damaske, D., eds. GANOVEX VIII-ITALIANTARTIDE XV Antarctic Expedition1999-2000,103-108. (Terra Antarctica Reports No. 5.)Google Scholar
Lefevre, F., Lyard, F. H., leProvost, C. and Schrama, E. J. O.. 2002. FES99: a global tide finite element solution assimilating tide gauge and altimetric information. J. Atmos. Ocean. Tech., 19, 1345-1356.Google Scholar
Legrésy, B. and Remy, F.. 1997. Altimetric observations of surface characteristics of the Antarctic ice sheet. J. Glaciol., 43(144), 265-275. (Erratum: 43(145), p. 596)CrossRefGoogle Scholar
Lingle, C.S., Hughes, T. J. and Kollmeyer, R. C.. 1981. Tidal flexure of Jakobshavns glacier, West Greenland. J. Geophys. Res., 86(B5), 3960-3968.Google Scholar
Pötzsch, A., Legresy, B., Korth, W. and Dietrich, R.. 2000. Glaciological investigations of Mertz Glacier, East Antarctica, using SAR interferometry and field observations. In ERS-ENVISATSymposium “Looking down to Earth in the New Millennium’, 16-20 October 2000, Gothenburg, Sweden. Proceedings. Noordwijk, European Space Agency, CD-ROM/500poetz. pdf. (ESA SP-461)Google Scholar
Rabus, B.T. and Lang, O.. 2002. On the representation ofice-shelfgrounding zones in SAR interferograms. J. Glaciol., 48(162), 345-356.Google Scholar
Riedel, B., Nixdorf, U., Heinert, M., Eckstaller, A. and Mayer, C.. 1999. The response of the Ekstromisen (Antarctica) grounding zone to tidal forcing. Ann. Glaciol., 29, 239-242.Google Scholar
Rignot, E. 1996. Tidal motion, ice velocity and melt rate of Petermann Gletscher, Greenland, measured from radar interferometry. J. Glaciol., 42(142), 476-485.Google Scholar
Schmeltz, M., Rignot, E., Dupont, T. K. and MacAyeal, D. R.. 2002. Sensitivity of Pine Island Glacier, West Antarctica, to changes in ice-shelf and basal conditions: a model study. J. Glaciol., 48(163), 552-558.CrossRefGoogle Scholar
Vaughan, D. G. 1995. Tidal flexure at ice shelf margins. J. Geophys. Res., 100 (B4), 6213-6224.Google Scholar
Wendler, G., Ahlnas, K. and Lingle, C. S.. 1996. On Mertz and Ninnis Glaciers, East Antarctica. J. Glaciol., 42(142), 447-453.CrossRefGoogle Scholar
Figure 0

Fig. 1. ERS SAR mosaic of the Mertz Glacier area. European Space Agency single look complex data were used to make this high-resolution mosaic. Straight white lines correspond to altimeter tracks, and the straight black dashed line to the radioecho sounding (RES) profile. The frames of the SAR scenes are indicated, and section A–B corresponds to the profile used to fit the tidal vertical bending. The curved white line corresponds to the grounding-line position recovered from the interferograms, the curved black line to the coastline and grounding line given in the Antarctic Digital Database, and the straight black dashed line to the RES flight-line. The black dotted lines correspond to the boundary of the ice tongue considered as effective for the mechanical bending.

Figure 1

Fig. 2. Comparison of the different remote-sensing measurements of sea-level variations vs the tide (FES99 model) + IBE (computed from Penguin Point pressure gauge). All values are in m. The comparison is made with GPS data (a), ERS altimetry (b) and ERS SAR interferometry (c). m is the regression coefficient, r is the correlation and σ is the rms data dispersion.

Figure 2

Table 1. Dates and conditions of acquisition of the SAR data

Figure 3

Fig. 3. Interferogram corresponding to the May (158) ERS SAR scenes. The many parallel fringes visible on the tongue reveal the horizontal flexure of the tongue in between the two acquisitions. The grounding line is drawn in white.

Figure 4

Fig. 4. Modelling the horizontal elastic flexure of the tongue for a 1m s-1 lateral ocean current.

Figure 5

Fig. 5. Double-difference interferograms showing the horizontal flexure of the tongue. One fringe corresponds to 28 mm distance variation in the satellite line of sight. (a) Double difference Ii. Both vertical flexure at the grounding line and lateral flexure of the tongue are observed. (b) Double difference I2. Mainly vertical flexure near the grounding line is present.

Figure 6

Fig. 6. GPS measurements on the ice tongue. The measurement point is situated at the edge of the tongue (Fig. 1). (a) The horizontal position along the 4 days of observation. Abscissa is along-flow distance from the first measurement; ordinates are across flow distance (dark lines) and height (light lines).The right-side plots correspond to the same along-/across-flow position, but with the average velocity component removed. (b) The height level (light grey lines), along-flow velocity (black dots) and height change or time derivative of the height (triangles).Do Y, day of year. (c, d) The across-flow position vs the height level for ascending (c) and descending tides (d). +/− values correspond to eastward/westward respectively. Different dots shades correspond to different days of observation. (e, f) The along-flow speed vs the height level for ascending (e) and descending tides (f).