KINEMATICS AND EFFICIENCY OF STEADY SWIMMING IN ADULT AXOLOTLS (AMBYSTOMA MEXICANUM)

Transcription

1 The Journal of Experimental Biology 2, (1997) Printed in Great Britain The Company of Biologists Limited 1997 JEB KINEMATICS AND EFFICIENCY OF STEADY SWIMMING IN ADULT AXOLOTLS (AMBYSTOMA MEXICANUM) KRISTIAAN D AOÛT* AND PETER AERTS Department of Biology, University of Antwerp, Universiteitsplein 1, B-261 Antwerp, Belgium Accepted 18 April 1997 The kinematics of steady swimming at a wide range of velocities was analysed using high-speed video recordings (5 frames s 1 ) of eight individuals of Ambystoma mexicanum swimming through a tunnel containing stationary water. Animals in the observed size range ( m total body length) prefer to swim at similar absolute speeds, irrespective of their body size. The swimming mechanism is of the anguilliform type. The measured kinematic variables the speed, length, frequency and amplitude (along the entire body) of the propulsive wave are more similar to those of anguilliform swimming fish than to those of tadpoles, in spite of common morphological features with the latter, such as limbs, external gills and a tapering tail. The swimming speed for a given animal size correlates linearly with the tailbeat frequency (r 2 =.71), whereas the wavelength and tail-tip amplitude do not correlate with this variable. The shape of the amplitude profile along the body, however, is very variable between the different swimming bouts, even at Summary similar speeds. It is suggested that, for a given frequency, the amplitude profile along the body is adjusted in a variable way to yield the resulting swimming speed rather than maintaining a fixed-amplitude profile. The swimming efficiency was estimated by calculating two kinematic variables (the stride length and the propeller efficiency) and by applying two hydrodynamic theories, the elongatedbody theory and an extension of this theory accounting for the slope at the tail tip. The latter theory was found to be the most appropriate for the axolotl s swimming mode and yields a hydromechanical efficiency of.75±.4 (mean ± S.D.), indicating that Ambystoma mexicanum swims less efficiently than do anuran tadpoles and most fishes. This can be understood given its natural habitat in vegetation at the bottom of lakes, which would favour manoeuvrability and fast escape. Key words: Ambystoma mexicanum, axolotl, kinematics, efficiency, anguilliform swimming. Introduction Swimming is a complex biological phenomenon comprising interactions between neural stimulation, muscle function, fluid dynamics and body kinematics (Carling et al. 1994). All these factors influence each other, making animal swimming a challenging topic for scientific inquiry [for reviews, see Lighthill (1975) on hydromechanics, Videler (1993), Wardle et al. (1995) and Van Leeuwen (1995) on kinematics and muscle function, and Carling et al. (1994) on neural stimulation]. Vertebrate swimming has been studied most extensively in fishes. Less is known about swimming in other vertebrate classes. In amphibians, the only primary swimmers besides the fishes, swimming remains largely unstudied. For undulatory swimming fishes, a number of qualitatively distinct swimming modes have been defined (Breder, 1926). This classification is based on the number of waves present on the body during steady swimming and ranges from anguilliform (more than half a wave) to ostraciiform (only the tail blade oscillates). Many studies have investigated the kinematics of steady swimming, and specific differences have been linked to body shape, muscle stimulation pattern and the animal s ecological niche (Wardle and Videler, 1994). In contrast, there is a lack of basic information regarding amphibian swimming, with the exception of anuran tadpoles (Wassersug and Hoff, 1985; Wassersug, 1989). Recent work suggests that the swimming mode of anuran tadpoles differs significantly from that of fish and that it appears to be tailored to their larval development (Liu et al. 1996). The few available studies on salamanders (Caudata) suggest that they show an undulatory swimming mechanism with kinematics resembling those of anguilliform fishes (Blight, 1976; Hoff et al. 1989; Frolich and Biewener, 1992). The present study addresses steady swimming over a range of velocities in the adult axolotl Ambystoma mexicanum, an obligate aquatic, neotenic salamander. Its larval morphology differs in many aspects from that of typically anguilliform * kdaout@uia.ua.ac.be.

2 1864 K. D AOÛT AND P. AERTS fishes, being more similar to that of tadpoles in possessing a distinct head region, a tapering tail and limbs. Moreover, external gills are present. We analysed the basic kinematics and hydromechanical efficiency of adult axolotls over as wide a swimming speed range as possible, and a functional comparison was made between the neotenic salamander s swimming mechanism and that of both other amphibians and fishes. Materials and methods Eight adult individuals of the neotenic Mexican axolotl Ambystoma mexicanum Shaw, 1789 (Amphibia: Caudata: Ambystomatidae) were used in experiments. Their total body length (L) and mass ranged from.135 to.238 m and from 22.2 to 92.1 g, respectively (Table 1). During experiments, the animals were placed in an aquarium consisting of two open tanks (.5 m.5 m.4 m) connected by a closed 1 m long glass tunnel (.15 m wide.1 m high). A mirror was placed at 45 above the tunnel to permit filming of lateral and dorsal views simultaneously. The water temperature was 2±2 C. The animals were video-taped at 5 frames s 1 using an NAC HSV-1 high-speed video system, lit by three 65 W halogen spotlights (Trilite). They were stimulated to enter the tunnel section of the aquarium by gently touching the tail with a blunt stick. The animals then usually swam voluntarily through the tunnel to the opposite aquarium at the speed they preferred, with the limbs and gills fully adducted. The maximal burst speed might well be much higher than the speeds recorded here. For example, in Ambystoma californiense (total body length approximately.5 m), burst speeds of more than 15 L s 1 have been measured (Shaffer et al. 1991). For detailed analysis, we selected 28 sequences (out of approximately 6 recorded) using the following criteria: (1) a constant swimming speed (i.e. a virtually rectilinear time displacement plot), (2) a straight swimming direction (i.e. a virtually straight path of the least-oscillating body point), (3) positioning more than one body width away from the tunnel side walls and (4) showing at least three complete swimming cycles. The 28 sequences were also chosen to cover as wide a velocity range as possible. Our analyses were based on the dorsal views of the swimming sequences only. The lateral views were used to check the animal s vertical position in the water during the swimming sequence but yielded no further information since the animals always swam along the bottom of the tunnel, with their ventral body surface in close proximity to the tunnel floor. The position of the dorsal midline was digitised every four frames (resulting in an effective sampling frequency of 125 frames s 1 ) using an NAC-1 XY coordinator connected to a computer. The dorsal midline was always clearly identifiable owing to the presence of a dorsal finfold over almost the entire body and tail region. In addition, preliminary studies in our laboratory showed that the kinematic variables did not differ significantly when either the directly digitised midline or a midline calculated from the animal s outline was used. To obtain the desired variables from the raw data files, the following procedures were used. (1) The animal was described in each frame by 2 segments of equal length (.5 L), delimited by 21 equally spaced points numbered from head to tail and referred to below as body points. The slope at the tail tip was then calculated, i.e. the slope between the twentieth and twenty-first body point (Cheng and Blickhan, 1994). (2) The least-oscillating body point was then determined. The path of this body point, usually body point 5 (i.e. at.2 L), very closely approximated a straight line in all sequences analysed and defined the swimming direction (the x-axis, estimated using a linear regression through the path of the least-oscillating body point). (3) The positions of the body points in each frame were plotted with the x-axis vertical, i.e. with the animal swimming in the positive x-direction (upwards on the plots) and with each frame shifted a known distance to the right (Fig. 1). From these normalised plots, a number of measurements were taken to obtain the following kinematic variables (technique modified from Videler and Wardle, 1978; Fig. 1): swimming speed U (m s 1 ), wave speed V (m s 1 ), wave frequency f (Hz), wave amplitude A (m) and wavelength λ (m). U, V, A and λ were Table 1. Morphological variables and observed swimming speeds of Ambystoma mexicanum Ind. SVL Cloaca Mean U Range U Mean U sp Range U sp code L (m) (m) BP Mass (g) N (m s 1 ) (m s 1 ) (L s 1 ) (L s 1 ) Range Re A B C D E F G H Ind. code, individual code; L, total body length; SVL, snout vent length; cloaca BP, body point closest to cloaca (see Materials and methods for further details); N, number of observations used for the complete kinematic analysis; U, swimming speed; U sp, body-length-specific swimming speed; Re, Reynolds number (based on L).

3 Kinematics of axolotl swimming 1865 Fig. 1. Technique for the determination of some of the kinematic variables measured (adapted from Videler and R1 Wardle, 1978). The first half-cycle of a wave is shown. Each curve shows the body of the animal divided into 2 α segments of equal length (connected by R the 21 body points) and represents one frame of the swimming sequence. The R1 first (uppermost) body point is the snout β tip. Subsequent frames (time interval B R2,2 8 ms) are shifted slightly to the right. R, β A line of reference (fixed position); R1, linear regression line through the leastoscillating body point, used to calculate R2,1 the swimming speed U; R1, R1 shifted vertically to enable the measurement of angle β; R2,1, R2,2, lines representing the successive wave fronts of the propulsive wave; β, angle used to calculate the wave speed V; A, segment used to calculate the wavelength λ, i.e. the distance between successive wave fronts; B, segment used to calculate the frequency f, i.e. the period between two successive wave fronts. For further details see text. scaled to the total body length, yielding specific speed U sp, specific wave speed V sp, specific amplitude A sp and specific wavelength λ sp. The reciprocal of λ sp gives the number of swimming waves present on the body at any instant during steady swimming. Although the use of the Reynolds number (Re) instead of body-length-specific values would allow direct comparison of kinematic variables under similar dynamic conditions (see Table 1), we chose instead to use the latter values. Since the swimming speed obtained at a given wave frequency is proportional to the body length (see Grillner and Kashin, 1976), dividing the swimming speed by the body length to give the specific swimming speed allows direct comparison of kinematic variables between differently sized individuals. We measured the wave amplitude A of every body point; the graph representing these amplitudes as a function of the position along the body is referred to below as the amplitude profile. A is defined as the maximal lateral deflection of a particular body point (or half the distance between the two extreme lateral positions of the body point). There was a large variation in A between the different sequences, but this could not be explained by individual variation (analysis of variance on body points 1, 6, 12 and 21, and on the intercept and slope of the line connecting body points 8 and 21, as an estimate of the steepness of the increase in amplitude down the body). The kinematic variables allowed for the estimation of the mechanical swimming efficiency in four different ways. First, the propeller efficiency, or slip, is given by the ratio U/V. It represents the amount of mechanical energy lost as slip, with the lowest losses as U/V approaches 1. Second, the specific stride length represents the number of body lengths travelled within one complete swimming cycle and is calculated as U/fL. Third, the hydrodynamic efficiency was calculated using Lighthill s (196) elongated-body theory with the efficiency µ EBT =(1+U/V)/2. This theory assumes that the animals use an undulatory mechanism with the wave amplitude increasing caudally, have an elongated (slender) body shape and zero slope at the tail tip. The first two of these assumptions are met by the axolotl. Fourth, we used Cheng and Blickhan s (1994) adapted elongated-body theory, which includes non-zero slopes at the tail tip, as were observed in the present study. The efficiency in this case, µ AEBT, is µ EBT decreased by an amount due to this tail-tip slope. According to Cheng and Blickhan (1994), we calculated the efficiency as: 1 1 µ AEBT = (1+β) α β 2 1 β with β being the propeller efficiency, U/V, and: 1 λ L A t α = 1, 2π L L A m with A t and A m being the amplitudes at the tail tip and at a distance L anteriorly, respectively. Note that equation 1 differs from Cheng and Blickhan s (1994) original equation 15 in the sign in the denominator of the final term. We assume this was a typographical error in the Cheng and Blickhan (1994) paper, as the application of equation 1 yields the correct results. In addition to a detailed kinematic analysis of 28 swimming sequences, we also measured the swimming speed in 15 other sequences (using simple distance/time measurements) to provide a larger data set for U and U sp. Results General kinematic description of axolotl steady swimming During steady swimming in Ambystoma mexicanum, waves of increasing amplitude travel down the body (see Fig. 2). At any instant, more than half a complete wave is visible on the body. According to Breder s (1926) terminology, the axolotl s swimming mechanism is thus of the anguilliform type. In 15 swimming sequences, swimming speed (U) ranged from.16 to.61 m s 1 and specific swimming speeds U sp, (1) (2)

4 1866 K. D AOÛT AND P. AERTS Fig. 2. Dorsal-view images illustrating one complete steady swimming cycle in Ambystoma mexicanum. The grid at the bottom of the tunnel consists of.5 m.5 m squares. Frame interval, 2 ms. Note that the external gills and the limbs are folded against the body during steady swimming from.446 to 3.44 L s 1 (Fig. 3). Almost the entire speed range was covered in the 28 sequences that were selected for detailed analysis (see Tables 1, 2). At lower speeds than those we observed, axolotls do not use an undulatory swimming mechanism but instead use their limbs to crawl along the bottom. There is a strong bias towards slower swimming speeds (Fig. 3): more than half of the observations (82 out of 15) lie in the range L s 1. The mean relative swimming speeds, as well as the minimum relative speeds and maximum relative speeds, decrease as the body length increases: animals within the observed size range prefer to swim at equal absolute, rather than relative, swimming speeds (Fig. 4). This preferred swimming speed was also reflected in the 28 sequences that were selected for the detailed analysis (see Tables 1, 2; Figs 5 7). Analysis of the propulsive wave The main characteristics of the propulsive wave are its wavelength λ, frequency f, speed V (where V=λf) and amplitude A. The specific wavelength λ sp did not correlate significantly with specific swimming speed U sp and was.58±.11 L (mean ± S.D., N=28) for all individuals and swimming speeds (Fig. 5). Number of observations U sp (L s 1 ) Fig. 3. The distribution of observed specific swimming speeds (U sp) in 15 swimming bouts of eight axolotls. Note the bias towards specific swimming speeds of approximately 1. L s 1. Maximal burst speeds might well be much higher than those observed during steady swimming. Swimming speed, U (m s 1 ) Maximum 16 Mean 72 Minimum Body length, L (m) Fig. 4. Relationship between swimming speed U and body length L. Mean, minimum and maximum swimming speeds observed are shown for the four individuals with the most observations (i.e. more than 15; the number of observations is given above the mean value). Note that animals of very different sizes evidently prefer to swim at very similar speeds.

5 Kinematics of axolotl swimming 1867 Table 2. Kinematic variables, swimming speeds and calculated efficiencies of 28 sequences of Ambystoma mexicanum swimming Specific Ind. U sp V sp A(L L) A(L) λ sp stride Propeller Sequence code (L s 1 ) (L s 1 ) (m) (m) (L) length efficiency µ EBT µ AEBT 1 F C B D E F F F F F F F F F F F H G H G G F A A A A A A Mean S.D Ind. code, individual code referring to Table 1; U sp, specific swimming speed; V sp, specific (backward) wave speed; A(L L), amplitude at the 9 %L position; A(L), amplitude at the tail tip; λ sp, specific wavelength; µ EBT, efficiency according to the elongated-body theory (Lighthill, 196); µ AEBT, efficiency according to Cheng and Blickhan's (1994) adapted elongated-body theory; S.D., standard deviation. Wave frequency f correlated linearly with U sp : f=1.56u sp (r 2 =.71, P<.1; Fig. 6). Because λ sp did not change with U sp, a similar correlation was also found between the specific wave speed V and the wave frequency f (not shown). The amplitude of the propulsive wave is commonly described by the amplitude at the tail tip (Fig. 7). While this is relevant for animals, such as most fish, whose amplitude profile remains constant, for Ambystoma mexicanum we found high variability in the amplitude profile during different swimming bouts, even at similar velocities (Fig. 8). Although the absolute amplitudes varied to a large extent (for example, tail-tip amplitude varied by more than twofold), some characteristics were similar in all amplitude profiles. First, the least-oscillating point was situated at 15±6 % L (mean ± S.D. N=28). Second, the amplitude increased from this point towards the snout as a passive result of movements in the trunk and tail region due to the immobility of the head neck connection (Duellman and Trueb, 1986). Third, the amplitude increased greatly towards the tail tip (which has the greatest amplitude in all sequences at.97±.22 L, mean ± S.D., N=28). Swimming efficiency Several different variables can be used to estimate the mechanical efficiency of swimming animals. Simple estimates are the stride length (or distance travelled during one swimming cycle) and the propeller efficiency (also known as the slip factor), which is the ratio of the swimming speed to the propulsive wave speed, U/V (Webb, 1975). A longer stride length and a higher propeller efficiency are assumed to indicate a higher mechanical swimming efficiency. In Ambystoma mexicanum, the specific stride length is.35±.8 (mean ± S.D., N=28) and the propeller efficiency is.6±.9 (mean ± S.D., N=28; see Table 2; Fig. 9). These values are rather, but

6 1868 K. D AOÛT AND P. AERTS Specific body wavelength, λsp (L) Fig. 5. Specific body wavelength λ sp as a function of specific swimming speed U sp. Mean ± S.D. λ sp was.58±.11 L and was not correlated significantly with U sp. not extremely, low when compared with those of fishes (see Discussion). While stride length and propeller efficiency can be used for comparative purposes, hydrodynamic theories allow the estimation of the hydrodynamic propulsive efficiency or the ratio of the power resulting in propulsion to the mechanical power delivered by the animal. In the present study, we calculated the hydrodynamic efficiency µ using both Lighthill s (196) classical elongated-body theory (µ EBT ) and Cheng and Blickhan s (1994) expansion of this theory to include the slope at the tail tip (µ AEBT, see Table 2; Fig. 1). Other theories (Lighthill, 1971; Yates, 1983; Cheng et al. 1991) are less suited to the analysis of anguilliform swimming with large wave amplitudes. Mean µ AEBT is.75±.4 (mean ± S.D., N=28), which is 6 % lower than mean µ EBT (.8±.5; mean ± S.D., N=28). The influence of the slope at the tail tip thus cannot be neglected in Ambystoma mexicanum, and we consider the adapted elongated-body theory to be the most appropriate for assessing the axolotl s swimming efficiency (see Discussion for details). We found some variation in swimming efficiency between different swimming bouts (Table 2; Fig. 1), but this could not be explained by an individual effect (tested using ANOVA) and was not related to specific swimming velocity (see Fig. 1). Wave frequency, f (Hz) Fig. 6. Wave frequency f as a function of specific swimming speed U sp. A linear regression fitted through the data points (r 2 =.71; P<.1) indicates that, as in most fishes, swimming speed is directly related to wave frequency. Specific tail-tip amplitude, Asp (L) Fig. 7. Specific tail-tip amplitude A sp as a function of specific swimming speed U sp. Note the high variability in tail-tip amplitudes. See text for further details. Discussion The present study has shown that the adult axolotl is an anguilliform swimmer, as defined for fish by Breder (1926) and characterised by the occurrence of more than half a wave on the body at any time during steady swimming. The adult axolotl differs in several morphological characteristics from anguilliform fish such as the eel (Grillner and Kashin, 1976) and the lamprey (Williams et al. 1989). The most striking differences are the presence of limbs, which are used for crawling at low speeds (<.5 L s 1 ), and the large feather-like external gills. The limbs and the gills are folded against the body during swimming, probably under active control (Hoff et Specific tail-tip amplitude (L) Distance from snout (%) Fig. 8. Graph representing the envelope containing all 28 observed amplitude profiles (shaded area). Two sample amplitude profiles are highlighted (solid and dashed lines) to demonstrate the variability between different swimming bouts: both show approximately the same tail-tip amplitude, but have markedly different profiles. In one case (solid line), the amplitude is relatively low in the trunk and relatively high in the tail; in the other case (dashed line), the opposite is true.

7 Kinematics of axolotl swimming Swimming efficiency Swimming efficiency Fig. 9. Swimming efficiency as a function of specific swimming speed U sp. Efficiency is calculated as the specific stride length (filled circles) and the propeller efficiency or slip factor (open circles). al. 1989). The effect of adducted limbs on swimming kinematics and efficiency is likely to be small. Liu et al. (1996) have found that, in tadpoles before metamorphosis, the hindlimbs impose very little handicap on propulsion. The kinematic variables determined here for swimming axolotls can be compared quantitatively with those of fishes measured in previous studies. Videler (1993) pooled kinematic data from a variety of fish species and found that, in general, U=.71f (r 2 =.93; N=45). In Ambystoma mexicanum, we found that U sp.64f 1.91, indicating that the axolotl needs a markedly higher tailbeat frequency than do fishes in order to attain a particular swimming speed. The mean specific tail-tip amplitude determined here for Ambystoma mexicanum (.97±.22 L) agrees very well with Videler s (1993) general finding for fishes of.1 L. However, axolotl and fish swimming kinematics differ in the distribution of the wave amplitude along the body. While most kinematic parameters relate significantly to the specific swimming speed, the observed amplitudes clearly do not and, more puzzlingly, show a wide variability throughout the speed range investigated (Fig. 7). For example, we observed a more than twofold difference in tail-tip amplitude between two swimming bouts by the same animal at swimming speeds and frequencies that differed by only.7 % and 9.3 %, respectively. This unexpected and puzzling result can probably only be explained by the cumulative effects of several factors. First, small variations in swimming speed may have occurred, and these can have a marked effect on amplitude (Webb, 1975). Second, the noise in our amplitude data may account for some of the observed variability. Third, it is possible that increased drag due to subtle changes in gill or limb posture leads to an increased amplitude requirement in order to keep the swimming speed constant. However, we did not observe such changes during qualitative inspection of the high-speed video recordings. We believe that the three explanations given above are likely to explain only part of the amplitude variability, and we Fig. 1. Swimming efficiency as a function of specific swimming speed U sp. Efficiency is calculated using Lighthill s (196) elongatedbody theory (filled circles) and using Cheng and Blickhan s (1994) adapted elongated-body theory (open circles). Note that, in every case, the adapted elongated-body theory gives lower efficiency values because this theory takes into account the effect of the slope at the tail tip. therefore propose an alternative, additional hypothesis. In most studies, the amplitude profile along the body is considered to be constant, and therefore a description of the propulsive wave amplitude consists only of the amplitude at the tail tip. This is justifiable for most fish, in which the amplitude profile does not change drastically. However, it is clear from the present results that the axolotl has a highly variable amplitude profile and, therefore, that it is not sufficient to consider the tail-tip amplitude alone: one has to consider the entire amplitude profile. Indeed, for an anguilliform swimmer delivering thrust along almost the entire body length (Wardle and Videler, 1994; D Août et al. 1996), this makes good sense. The same amount of propulsion, needed to counterbalance the drag (a prerequisite for steady swimming), can then be delivered in different ways, using different amplitude profiles with different tail-tip amplitudes. It is unlikely that different amplitude profiles have equal efficiency, and one would expect evolution towards one optimal solution. So why has this not occurred in Ambystoma mexicanum? Ambystomatid salamanders live in the dense vegetation of shallow lakes (Smith and Smith, 1971; Scott, 1981). In such an environment, the ability to manoeuvre and to perform short strikes is important. These highly unsteady locomotor patterns impose totally different requirements on the locomotor system from those of steady swimming. It is likely then that adaptation towards an efficient steady swimming mechanism is not favoured in this environment. This hypothesis is further supported by the relatively low hydromechanical efficiency of steady swimming in Ambystoma mexicanum. Great care must be taken when calculating values for efficiencies. The classical parameters, propeller efficiency (slip) and stride length, are easy to calculate but, as they do not yield the true efficiency, they can only be used to compare the swimming efficiency of animals

8 187 K. D AOÛT AND P. AERTS using a similar swimming style. Their use is based on the assumption that the lower the slip and the higher the distance travelled per cycle, the higher the mechanical swimming efficiency will be. Hydrodynamic theories, in contrast, aim to provide values of the true hydrodynamic efficiency (the ratio of energy resulting in propulsion to mechanical energy delivered by the animal), but one must be aware that the available theories each have specific drawbacks. The widely used elongated-body theory (Lighthill, 196) is the most appealing because of its simplicity, the efficiency µ EBT being given by µ EBT =(1+U/V)/2. However, its use is confined to small-amplitude motions, a point addressed by Lighthill s (1971) large-amplitude elongated-body theory. While the latter theory has proved very useful for carangiform swimmers with a distinct trailing edge, for which the theory was developed, it is of questionable use for anguilliform swimmers with a tapering tail and has therefore not been used here. Both theories are based on a balance between (1) the moments imparted to the water by the animal and shed in the post trailing-edge wake and (2) the propulsion. Since anguilliform swimmers are thought to deliver propulsion along most of the body length (D Août et al. 1996; Wardle and Videler, 1993; Wardle et al. 1995), the more recent three-dimensional waving plate theory (Cheng et al. 1991) seems more attractive for the study of anguilliforms because it is based on the vortex ring around the body (represented by a sequence of articulating plates of different heights). However, this theory requires more complex input from the experiments and is much more elaborate to calculate, which might be the reason why Lighthill s (196, 1971) theories remain standard in most studies. Moreover, a recent computational fluid-dynamic approach to tadpole swimming (Liu et al. 1996) yielded very similar efficiencies to Lighthill s (196) elongated body theory. Cheng and Blickhan (1994) adapted Lighthill s (196) elongated-body theory to remove the assumption of zero slope at the tail tip. Indeed, in the axolotl (and also in most fish), the tail tip does not move parallel to the sagittal plane during most of the swimming cycle, but has a certain slope (see Fig. 2). This slope will result in a decreased efficiency relative to that calculated according to the elongated-body theory. Cheng and Blickhan s (1994) results demonstrate that the adapted elongated-body theory approaches the results from the three-dimensional waving plate theory and remains simple to apply. We consider this theory the most appropriate for estimating the axolotl s hydromechanic swimming efficiency. Only Cheng and Blickhan (1994) have used the adapted elongated-body theory to date, and the axolotl s efficiency can therefore best be compared with that of other swimmers using published values for propeller efficiency and stride length and using elongated-body theory. Propeller efficiency and stride length Hoff et al. (1989), using small Ambystoma larvae, reported much higher relative swimming speeds than those in the present study, which agrees with our observation in adult axolotls that the smaller the individual, the higher the specific swimming speed. The larvae swam with a propeller efficiency of.65 at swimming speeds greater than 6 L s 1. At lower velocities, the propeller efficiency was slightly lower, making our results (mean propeller efficiency.6) consistent with their findings. Anuran tadpoles (of Bufo, Rana and Xenopus) perform markedly better, with propeller efficiencies of (Wassersug and Hoff, 1985; Hoff and Wassersug, 1986). Much better propeller efficiencies are achieved by several fishes, e.g..8 for Salmo (Webb, 1975). A similar pattern emerges for the stride length data. Larval Ambystoma have a stride length of.5 L (Hoff et al. 1989), Bufo and Rana tadpoles.6 L (Wassersug and Hoff, 1985; Hoff and Wassersug, 1986) and subcarangiform fishes.7 L (Bainbridge, 1958). Our result for Ambystoma mexicanum,.35 L, is markedly lower. The knifefish, Notopterus notopterus, has an even shorter stride length of.3 L (Blake, 1983). Hydrodynamic theories Wassersug and Hoff (1985) calculated the swimming efficiency of anuran tadpoles to be.83, slightly higher than that found here for Ambystoma mexicanum (.8). They also found that efficiency tended to be lower at lower swimming speeds. Liu et al. (1996) calculated the hydrodynamic efficiency of anuran tadpoles using a computational fluiddynamic (CFD) approach and found an efficiency of.82, in very close agreement with the elongated-body theory values for these tadpoles (.83). Fishes and some mammals swim much more efficiently and elongated-body theory efficiencies greater than.9 are not uncommon. The kinematics and the estimates of mechanical swimming efficiency reported in the present study show that Ambystoma mexicanum uses an anguilliform swimming mechanism with a highly variable wave amplitude profile and with a low efficiency. This is consistent with their ecological niche and the likely constraints imposed by this environment. We thank Dr Frits De Vree for his expertise and support during the research in his laboratory, and Dr R. Van Damme and two anonymous reviewers for their input and constructive comments on an earlier version of the manuscript. This research was supported by grant from the Flemish Institute for the Encouragement of Scientific Research in Science and Technology (I.W.T.) to K.D. P.A. is a research associate from the Fund for Scientific Research Flanders (F.W.O.). References BAINBRIDGE, R. (1958). The speed of swimming of fish as related to size and to the frequency and amplitude of the tail beat. J. exp. Biol. 35, BLAKE, R. W. (1983). Median and paired fin propulsion. In Fish Biomechanics (ed. P. W. Webb and D. Weihs), pp New York: Praeger Scientific. BLIGHT, A. R. (1976). Undulatory swimming with and without waves of contraction. Nature 264,

9 Kinematics of axolotl swimming 1871 BREDER, C. M. (1926). The locomotion of fishes. Zoologica 4, CARLING, J. C., BOWTELL, G. AND WILLIAMS, T. L. (1994). Swimming in the lamprey: modelling the neural pattern generation, the body dynamics and the fluid mechanics. In Mechanics and Physiology of Animal Swimming (ed. L. Maddock, Q. Bone and J. M. V. Rayner), pp Cambridge: Cambridge University Press. CHENG, J.-Y. AND BLICKHAN, R. (1994). Note on the calculation of propeller efficiency using elongated body theory. J. exp. Biol. 192, CHENG, J.-Y., ZHUANG, L.-X. AND TONG, B.-G. (1991). Analysis of swimming three-dimensional waving plates. J. Fluid Mech. 232, D AOÛT, K., AERTS, P. AND DE VREE, F. (1996). The timing of muscle strain and activation during steady swimming in a salamander, Ambystoma mexicanum. Neth. J. Zool. 46, DUELLMAN, W. E. AND TRUEB, L. (1986). Biology of Amphibians. New York: McGraw-Hill. FROLICH, L. M. AND BIEWENER, A. A. (1992). Kinematic and electromyographic analysis of the functional role of the body axis during terrestrial and aquatic locomotion in the salamander Ambystoma tigrinum. J. exp. Biol. 162, GRILLNER, S. AND KASHIN, S. (1976). On the generation and performance of swimming in fish. In Neural Control of Locomotion (ed. R. M. Herman, S. Grillner, P. S. G. Stein and D. G. Stuart), pp New York: Plenum. HOFF, K., HUQ, N., KING, V. A. AND WASSERSUG, R. J. (1989). The kinematics of larval salamander swimming. Can. J. Zool. 67, HOFF, K. AND WASSERSUG, R. J. (1986). The kinematics of swimming in larvae of the clawed frog, Xenopus laevis. J. exp. Biol. 122, LIGHTHILL, M. J. (196). Note on the swimming of slender fish. J. Fluid Mech. 9, LIGHTHILL, M. J. (1971). Large-amplitude elongated-body theory of fish locomotion. Proc. R. Soc. Lond. B 179, LIGHTHILL, M. J. (1975). Mathematical Biofluiddynamics. Philadelphia: Society for Industrial and Applied Mathematics, 281pp. LIU, H., WASSERSUG, R. J. AND KAWACHI, K. (1996). A computational fluid dynamics study of tadpole swimming. J. exp. Biol. 199, SCOTT, P. W. (1981). Axolotls. Neptune City, NJ: T. F. H. Publications. SHAFFER, H. B., AUSTIN, C. C. AND HUEY, R. B. (1991). The consequences of metamorphosis on salamander (Ambystoma) locomotor performance. Physiol. Zool. 64, SMITH, H. M. AND SMITH, R. B. (1971). Synopsis of the Herpetofauna of Mexico, vol. 1, Analysis of the Literature on the Mexican Axolotl. Augusta: Lundberg Press. VAN LEEUWEN, J. L. (1995). Review article: The action of muscles in swimming fish. Exp. Physiol. 8, VIDELER, J. J. (1993). Fish Swimming. London: Chapman & Hall. VIDELER, J. J. AND WARDLE, C. S. (1978). New kinematic data from high speed cine film recordings of swimming cod (Gadus morhua). Neth. J. Zool. 28, WARDLE, C. S. AND VIDELER, J. J. (1993). The timing of the electromyogram in the lateral myotomes of mackerel and saithe at different swimming speeds. J. Fish Biol. 42, WARDLE, C. S. AND VIDELER, J. J. (1994). The timing of lateral muscle strain and EMG activity in different species of steadily swimming fish. In Mechanics and Physiology of Animal Swimming (ed. L. Maddock, Q. Bone and J. M. V. Rayner), pp Cambridge: Cambridge University Press. WARDLE, C. S., VIDELER, J. J. AND ALTRINGHAM, J. D. (1995). Tuning in to fish swimming waves: body form, swimming mode and muscle function. J. exp. Biol. 198, WASSERSUG, R. J. (1989). Locomotion in amphibian larvae (or Why aren t tadpoles built like fishes?). Am. Zool. 29, WASSERSUG, R. J. AND HOFF, K. (1985). The kinematics of swimming in anuran larvae. J. exp. Biol. 119, 1 3. WEBB, P. W. (1975). Hydrodynamics and energetics of fish propulsion. Bull. Fish. Res. Bd Can. 19, WILLIAMS, T. L., GRILLNER, S., SMOLJANINOV, V. V., WALLÉN, P., KASHIN, S. AND ROSSIGNOL, S. (1989). Locomotion in lamprey and trout: the relative timing of activation and movement. J. exp. Biol. 143, YATES, G. T. (1983). Hydromechanics of body and caudal fin propulsion. In Fish Biomechanics (ed. P. W. Webb and D. Weihs), pp New York: Praeger Scientific.