Wing and body kinematics of free cruising flight are described for 37 species of Panamanian butterflies ranging over two orders of magnitude in body mass. Butterflies exhibit considerable diversity in body and wing shape, but morphological design is, in general, isometric. Wing loading and mean body diameter show positive allometry. The cruising flight of butterflies is characterized by low wingbeat frequencies (here averaging 11 Hz), stroke amplitudes averaging 103°, and forward speeds in excess of 1ms−1. Body angles during flight are close to horizontal, and stroke plane angles are correspondingly high. Advance ratios are typically greater than 0.9, indicating that the forward and flapping velocity vectors are of comparable magnitude. Flight speed scales with morphological parameters in general accordance with predictions based on isometric design. Interspecifically, no consistent correlation exists between wing kinematics and absolute flight speed. However, maximum positional angle and stroke amplitude tend to increase while body angle decreases with increased relative flight speed.

The wide diversity of wing size and shape in butterflies (Lepidoptera: Papilionoidea and Hesperioidea) suggests considerable variation in the biomechanics and aerodynamics of flight. Although variation in wing kinematics and the characteristic flight styles of particular butterfly species have been commented upon anecdotally, there is little quantitative description of these phenomena. Unlike most other insects, many butterflies exhibit considerable vertical oscillation during horizontal flight, giving rise to the characteristic fluttering flight path. The biomechanical implications of these oscillating flight paths have never been investigated. Similarly, the basic aerodynamics of butterfly flight has not been systematically studied. The low aspect ratio wings characteristic of many butterfly species are generally associated with a large induced drag and correspondingly high power requirements of flight. However, only a detailed analysis of the aerodynamic mechanisms operating during flight will permit a functional evaluation of butterfly wing shape.

The biomechanics and flight patterns of butterflies are also of considerable ecological interest. It has been suggested that the frequently erratic flight of butterflies helps them to evade would-be predators (e.g. Humphries and Driver, 1970; Edmunds, 1974). It has also been demonstrated that specific flight patterns and wing and body shapes are associated with unpalatability in tropical mimetic butterflies (Chai, 1986).

Recently, surveys of flight morphology and kinematics in bats and birds have helped to interpret flight styles and functional correlates of wing and body form (e.g. Norberg, 1981, 1986; Norberg and Rayner, 1987; Rayner, 1987; Viscor and Fuster, 1987). Equivalent studies for the diverse array of flying insects are not yet available. It is the purpose of this study to present kinematic data and a biomechanical interpretation of wing and body morphology for the butterfly superfamily Papilionoidea (Ackery, 1984). This paper describes, for an assemblage of Panamanian butterflies, wing and body motions determined using video filming and three-dimensional image reconstruction. Kinematic analysis is confined to wing motions relative to the butterfly body and to general body orientation in space; correlations between flight kinematics and various aspects of wing and body morphology are investigated. Additional papers will evaluate, for the same butterflies, vertical oscillations of the body and the aerodynamics and mechanical power requirements of flight.

Study site and butterflies

All work was carried out on Barro Colorado Island, a research station located in Lake Gatun, Republic of Panama. The butterfly fauna of Barro Colorado Island is described by Huntington (1932). Butterflies chosen for the study belonged principally to the families Papilionidae, Pieridae and Nymphalidae (identifications following DeVries, 1987). However, one small riodinid butterfly and one skipper (Hesperioidea: Hesperiidae) were filmed. Butterflies were collected in either the laboratory clearing or the forest of Barro Colorado Island, and were generally filmed on the same day. Those butterflies filmed on the following day were kept overnight in small net cages and were fed a solution of sugar and water (1:1 v/v).

Filming procedures

Filming was carried out in a small screened insectary (3.4m×3m×2.ini). Butterflies were released in the insectary and allowed to fly freely. A video camera (Panasonic AG-160) equipped with a high-speed electronic shutter (1/1000 s) was used to film at 30 frames s−1; the high-speed shutter permitted sharp resolution of the moving wings. Butterflies generally cruised back and forth along one side of the insectary, opposite which the video camera was positioned horizontally using a spirit level. The camera was mounted on a tripod at a fixed height of about 1.5 m, and was generally 2 m from the insectary wall. Ambient temperatures in the insectary during filming ranged from 30 to 35 °C and relative humidities from 82 to 95%. Ambient air motion in the screened insectary was negligible. After approximately 20 min of free flight, the butterfly was captured and killed by gently squeezing the thorax, with little or no hemolymph loss.

Morphological parameters

Morphological measurements on all butterflies were made within 10 min of filming; mass loss during this time was negligible. Basic morphological parameters determined for each insect included the following: body mass m (which includes the mass of both wing pairs), thoracic mass mt (expressed in both absolute terms and as a fraction of body mass,) wing length R, body length (expressed as a fraction of wing length) and total wing area S (the area of both wing pairs). Wing and body lengths were determined to ±0.1 mm using dial calipers. Wing base separation was also measured to ±0.1 mm, and was used to calculate the wingspan B (the distance between the extended wing tips) (see Norberg and Rayner, 1987). Wing area was determined with a leaf area meter (Delta-T Devices AMS), with the fore-and hindwings placed together in a position roughly equivalent to that observed in flapping flight (see Betts and Wootton, 1988). Areas of the fore- and hindwing were also determined separately (Sforc and Shind, respectively), and then-overlap during flapping flight Sover is given by (Sfore+Shind-S). In many butterflies, a portion of the hindwing folds over part or all of the abdomen and does not noticeably move when the wings flap. This overlap of the hindwing onto the body, Sbody, was determined by cutting and measuring this area separately from one of the two hindwings. Because this area does not function as part of the flapping wing surface, it was subtracted from the measurement of total wing area S. Wing loading pw (=mg/S, where g is gravitational acceleration) and aspect ratio were calculated for each insect.

While such basic morphological data as body mass and wing area are essential for a biomechanical analysis of flight, several non-dimensional parameters are also useful in evaluating the functional consequences of wing and body shape. In particular, the present morphological analysis makes considerable use of the moments about the wing base (and their non-dimensional radii) of wing area, wing mass and wing virtual mass, which is the mass of air accelerated by the moving wings. For greater details of the theory and methods involved, see Ellington (1984b).

Wing mass (the mass of both wing pairs) was measured to the nearest 0.1 mg and expressed as a percentage of the total body mass. The distribution of wing mass along the longitudinal wing axis was determined in the following manner. Beginning at the wing base, consecutive 5-mm strips were cut perpendicular to the wing length. Each strip was weighed to the nearest 0.1 mg. To correct for mass losses during the cutting and weighing procedure, the mass of each strip was multiplied by the ratio of the initial wing pair mass to the sum of masses of the individual strips. Non-dimensional radii of the first and second moments of wing mass [ and ] were determined using the strip masses. Using tracings of the wing pair outline, non-dimensional radii for the first three moments of wing area and of wing virtual mass were determined. Calculations of these area and virtual mass parameters were made using wing chords concentric to the wing base rather than chords perpendicular to the wing length, because of the low aspect ratio.of the wings (Betts and Wootton, 1988). A non-dimensional virtual wing mass (equal to the ratio of the wing virtual mass to that of a wing of equal R but with a constant wing chord equal to the mean wing chord) was also calculated for each wing pair. Virtual mass estimates are only approximate as these calculations assume a perpendicular wing chord, not the circular chord used here. A mean wing thickness ĥ, expressed as a percentage of the wing length, was also calculated.

The center of mass of the body was determined from masses of five body regions (head, two thoracic sections and two abdominal sections), assuming mass to be homogeneously distributed within each section (see R. Dudley and P. J. DeVries, in preparation). This procedure assumes that the center of mass is located on the longitudinal body axis. Because body length is large relative to mean body diameter (see Results), this is an appropriate assumption. The distance of the center of mass from the anterior tip of the body is termed , while is the distance of the center of mass to the wing base axis, and is the radius of gyration for the moment of inertia of the body about the wing base axis. These three parameters were expressed non-dimensionally as a fraction of body length L, as was the mean body diameter (Ellington, 1984b). Absolute values of the mean body diameter d were also calculated.

For comparison with results of previous workers, linear regressions were used to evaluate correlations between the body shape parameters and , the wing shape parameters and and the wing mass parameters and . Linear correlations between log-transformed values of the non-dimensional radii characterizing wing mass and shape were also performed to determine allometric power functions relating the different radii. Bivariate reduced major axis (RMA) correlations were determined between the log-transformed data of body mass and various wing and body dimensions to evaluate their allometry (see Rayner, 1985, for a discussion of appropriate statistical models in biomechanical analysis).

Flight kinematics

Video films were viewed using frame-by-frame playback on a video tape recorder (Hitachi VT 1570-A) connected to an AMDEK Video-300 VDU. For each butterfly, representative sequences of approximately horizontal, non-accelerating flight were identified through repeated viewing of the video tape. For each sequence, the following points on the butterfly image were manually digitized from each video frame: anterior and posterior tips of the body, wing base and both wing tips. The accuracy of the digitization process was checked by digitizing video images of variably sized circles. Typical distortion was approximately 4% of the true value for linear dimensions in the central regions of the screen where images were analyzed. Images near the screen edges were avoided. For all sequences, linear regressions of the coordinates of the wing base were used to determine the angle ξ between the flight path and horizontal. Sequences in which the absolute value of ξ was greater than 15° were rejected.

Determination of wing motions and body orientation in three dimensions requires geometric reconstruction of the wing tip position using two-dimensional film images. Mathematical details of the reconstruction procedure used here can be found in Ellington (1984c), and will be briefly summarized. By considering the wing tip positions at both ends of a half-stroke, the three-dimensional motion of the wing tip relative to the body can be described. Points of interest on the insect’s wing and body are transformed from a coordinate system based on the film image to a coordinate system based on the insect’s body. Wing motions are assumed to be symmetrical with respect to the longitudinal body axis, which is the case in forward flight. A computer program which encoded the reconstruction procedure was written in MS-BASIC for use with a Macintosh Plus microcomputer.

For purposes of kinematic analysis, it was necessary to determine a reference length for each flight sequence. The low filming speed of the camera precluded use of the maximum projected wing length of a sequence as the reference length (cf. Ellington, 1984c). Instead, the projected body length was used. For all filming, the camera was nominally perpendicular to the flight path of the insect (see Results), and the projected body length was equivalent to the maximum possible projected length. An average projected body length was calculated for each flight sequence, and was used to estimate a maximum projected wing length by scaling appropriately from the morphological data. Occasionally, otherwise promising flight sequences could not be analyzed because the hindwing obscured the posterior tip of the body for a substantial number of frames.

The kinematic analysis assumes that wing motions are confined to a stroke plane which forms an angle β relative to horizontal. The maximum positional angle ϕmax and the minimum positional angle ϕmin refer to the wing tip position in the stroke plane at the top and bottom of the half-stroke, respectively. Stroke amplitude Φ, the angular extent of motion in the stroke plane, is given by ϕmax–ϕmin; the mean positional angle is equal to (ϕmaxmin)/2. Body angle χ is the angle relative to horizontal made by a fine segment connecting the anterior and posterior tips of the body, when the insect is viewed laterally. Measurements of body angle were made when the wings were at the top (χt) and at the bottom (χb) of the half-stroke; the change in body angle through the half-stroke, Δχ (=χtχb) and an average body angle were also calculated.

Proper identification of two video frames representing the top and bottom of the half-stroke is critical for a correct determination of wingbeat kinematics. Because, in general, these frames did not immediately follow one another (unless the wingbeat frequency was exactly 15Hz), it was necessary to analyze composite wingbeats formed from non-consecutive frames separated by intervening wingbeats. One frame was chosen which depicted the wing at the top of the half-stroke, and then the video tape was advanced frame by frame until the wing was shown at the bottom of the half-stroke. Video sequences were viewed repeatedly to guarantee correct identification of the frames containing the top and bottom of the half-stroke. Multiple composite wingbeats were analyzed in each video sequence to evaluate short-term variation in kinematic parameters (see Results).

For each flight sequence, the wingbeat frequency n was determined by counting the number of video frames required to complete a wingbeat, with occasional visual estimation of the wing positional angle. Because of the filming speed of the video camera (30 frames s−1), the accuracy of this procedure was low, approximately ±2 Hz. Five estimates of the wingbeat frequency were made and used to determine a mean value for each flight sequence. A similar frame counting procedure was also used for each sequence to determine the ratio of upstroke duration to downstroke duration (US:DS). For each butterfly, the number of frames in all analyzed sequences which depicted gliding flight (characterized by little or no wing motion from frame to frame) was determined. This quantity was divided by the total number of frames to determine the relative extent of gliding flight, G.

For all video filming, the optical axis of the camera was horizontal, so that the body angles and stroke plane angles determined by the kinematic analysis were expressed relative to true horizontal. Values of the roll η about the longitudinal body axis were calculated for each composite wingbeat (see Ellington, 1984c). Yaw a, which describes the rotation of the insect about the vertical axis, was also determined. Values of η are positive when the insect rolls away from the camera, and positive a indicates a yaw towards the camera. A yaw of 0° indicates that the longitudinal body axis is orthogonal to the camera. Non-zero yaw indicates that the actual projected body length deviates from the maximum possible projected length.

Relative flight speed was determined by calculating the number of body lengths travelled horizontally by the wing base per second. Horizontal airspeed V was then expressed in absolute terms (ms−1) by scaling with the measured body length. Fight speeds were related allometrically to wing and body morphology using RMA correlations on log-transformed data. RMA correlations were also used to evaluate the dependence of wing and body kinematics upon forward airspeed.

Errors arising from perspective effects in the three-dimensional reconstruction of the video image were estimated from the formulae given by Ellington (1984c). Representative values for horizontal coverage of the film image and for the distance from the object plane to the focus of the lens were about 7 and 75 wing lengths, respectively. Maximum error arising from the neglect of perspective in reconstructed spatial angles of the wings was calculated to be less than I0, and the deviation from true horizontal of the horizontal plane based on the image projection was less than 2°.

Butterfly identification and morphometries

Morphological data were obtained for a total of 49 butterflies and one skipper, representing 37 species. Seven species were represented by two individuals (generally a male and a female), and three species were represented by three individuals. The taxonomic distribution of species was as follows: Papilionidae, 4; Pieridae, 5; Nymphalidae: Charaxinae, 2; Nymphalidae: Nymphalinae, 9; Nym-phalidae: Heliconiinae, 5; Nymphalidae: Melitaeinae, 1; Nymphalidae: Ithomii-nae, 2; Nymphalidae: Morphinae, 2; Nymphalidae: Brassolinae, 1; Nymphalidae: Satyrinae, 4; Riodinidae, 1; Hesperiidae, 1. Most butterflies were fairly large (average mass 267 mg) to ensure an adequate image size during filming, and there was little representation of butterflies of mass less than 90 mg.

Three separate flight sequences were analyzed for each of 15 individuals (see below), for which species identification, sex and general morphological parameters are given in Table 1. For the same butterflies, wing shape parameters are given in Table 2 and body parameters in Table 3. For comparative purposes, Tables 1–3 also include the minimum, mean and maximum values of all morphological parameters for the 50 butterflies examined. Because only one individual was generally used per species, intraspecific differences arising from geographical variation or sexual dimorphism were ignored. This variation is, however, likely to be small compared with the wide range of interspecific variation apparent in Tables 1–3.

Table 1.

Identification number (ID), species, sex, mass (m), wing length (R), aspect ratio (AR), wing loading (pw) and relative thoracic mass

(m^t)
of 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same morphological parameters for the full sample of 50 butterflies

Identification number (ID), species, sex, mass (m), wing length (R), aspect ratio (AR), wing loading (pw) and relative thoracic mass (m^t) of 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same morphological parameters for the full sample of 50 butterflies
Identification number (ID), species, sex, mass (m), wing length (R), aspect ratio (AR), wing loading (pw) and relative thoracic mass (m^t) of 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same morphological parameters for the full sample of 50 butterflies
Table 2.

Wing parameters of the 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same parameters for the full sample of 50 butterflies

Wing parameters of the 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same parameters for the full sample of 50 butterflies
Wing parameters of the 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same parameters for the full sample of 50 butterflies
Table 3.

Body parameters of the 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same parameters for the full sample of 50 butterflies

Body parameters of the 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same parameters for the full sample of 50 butterflies
Body parameters of the 15 butterflies for which three flight sequences were analyzed, and minimum, mean and maximum values of the same parameters for the full sample of 50 butterflies

Linear regressions were applied to the non-dimensional radii characterizing the wing shape and mass distribution, to compare the relationships for butterflies alone with those determined by Ellington (1984b) for a set of insects from various orders. Regressions relating to to and to were all significant at P<0.0001, and their slopes did not differ significantly from those previously published. The regression of on was also highly significant (P<0.0001), but 95% confidence limits of the slope fell outside that given previously. Least-squares regression was also performed on the non-dimensional radii of the first and second moments of body mass about the wing base axis, and . The resulting correlation was highly significant (P<0.001), and the slope of the regression did not differ significantly from that reported by Ellington (1984b).

In addition to linear correlations between morphological parameters, bivariate linear regressions of log-transformed data were used to derive allometric relationships between the non-dimensional parameters characterizing wing shape. Such allometric relationships between the following sets of variables were highly significant (P<0.001): and , and and and and and and . Three of these regressions were characterized by exponents significantly different from those determined by Ellington (1984b) for a diverse group of insect wings from six orders: on on and on No significant correlation was found between log-transformed values of and or between and .

For the full sample of 50 butterflies, the following morphological parameters increased isometrically with body mass (P<0.001): wingspan, wing length, body length, wing area, fore- and hindwing area, thoracic mass, wing mass, and fore- and hind wing mass. Wing loading (Fig. 1) and mean body diameter exhibited positive allometry. There was no significant correlation between aspect ratio and body mass (P>0.05).

Fig. 1.

Wing loading pw (in N m−2) vs body mass m (in kg) for 50 butterflies. The line shown is the reduced major axis correlation (r=0.74, P<0.0001); the equation is given by: pw=49.05m0.474. Wing loading exhibits positive allometry.

Fig. 1.

Wing loading pw (in N m−2) vs body mass m (in kg) for 50 butterflies. The line shown is the reduced major axis correlation (r=0.74, P<0.0001); the equation is given by: pw=49.05m0.474. Wing loading exhibits positive allometry.

Wing and body kinematics

Most butterfly species examined flew parallel to and within 20 cm of the sunlit fide of the insectary. Typically, the butterfly cruised continuously back and forth or minutes on end, and flight sequences were taken from the center of the generally horizontal flight trajectory. For six individuals, only weak fluttering flights in the corners of the insectary were filmed, and these butterflies were therefore used for morphological data only. From the other butterflies, one flight sequence was analyzed for each of 15 individuals, two such sequences for 14 individuals, and three flight sequences were analyzed for each of the remaining 15 butterflies. The mean number of video frames per analyzed flight sequence was 8.4 (equivalent to about 0.3 s of flight). Mean values of kinematic parameters for each sequence were calculated using an average of 2.4 composite wingbeats. Whenever either two or three separate flight sequences of a particular butterfly were analyzed, kinematic results from these sequences were averaged to yield representative values of the flight speed and other kinematic parameters. In such cases, statistical analysis was also performed separately using the kinematic data obtained for sequences showing the lowest, intermediate (if it existed) and highest flight speed.

For the 15 butterflies with three analyzed sequences each, average values of wing and body kinematic parameters are given in Table 4, as are mean values of the advance ratio J (=V/2ΦnR). For each kinematic parameter, standard deviations were derived from the three sequences of each butterfly, and mean values of these are also given in Table 4, as are means of the kinematic data and the range for all butterflies (N=44) for which at least one flight sequence was analyzed. Flight sequences used for analysis generally displayed non-rolled flight perpendicular to the camera axis. For all analyzed wingbeats, the average absolute values of the roll η and yaw a were 8° and 9°, respectively.

Table 4.

Mean kinematic parameters of the 15 butterflies for which three flight sequences were analyzed, and average, minimum and maximum values of the same parameters for all analyzed sequences (44 butterflies)

Mean kinematic parameters of the 15 butterflies for which three flight sequences were analyzed, and average, minimum and maximum values of the same parameters for all analyzed sequences (44 butterflies)
Mean kinematic parameters of the 15 butterflies for which three flight sequences were analyzed, and average, minimum and maximum values of the same parameters for all analyzed sequences (44 butterflies)

Kinematics and morphology

For the sample of 15 butterflies with three analyzed flight sequences, forward airspeed was positively correlated with a variety of morphological parameters, including body mass, thoracic mass, wing length, wingspan, body length, mean body diameter, wing area and wing loading (Table 5). No significant correlation was found between aspect ratio and forward airspeed (P>0.05). When these correlations were performed using alternately the lowest, intermediate or highest airspeeds (and the corresponding kinematic results) determined for each butterfly, the results were unchanged. Table 5 includes the predicted scaling of flight speed given isometric design (see Norberg and Rayner, 1987). The allometric dependence of flight speed on body mass, thoracic mass, body length and wing loading differed significantly from predicted values.

Table 5.

Allometric relationships between forward airspeed V and morphological parameters for the sample of 15 butterflies each with three analyzed flight sequences

Allometric relationships between forward airspeed V and morphological parameters for the sample of 15 butterflies each with three analyzed flight sequences
Allometric relationships between forward airspeed V and morphological parameters for the sample of 15 butterflies each with three analyzed flight sequences

Considering all butterflies for which at least one flight was analyzed, correlations of airspeed with wing length, mean body diameter, wing area and aspect ratio were not significant (P>0.05). Correlations of airspeed with body mass, thoracic mass, body length, wingspan and wing loading were, however, significant (P<0.01), although characterized by correlation coefficients much lower than those for correlations using the sample of 15 butterflies with three flight sequences each.

Of the kinematic parameters in Table 4, only the maximum positional angle ϕmax was significantly correlated with mean forward airspeed V for the sample of 44 butterflies (P<0.02). For the 15 butterflies with mean airspeeds and kinematic parameters derived from three separate sequences, none of the kinematic parameters in Table 4 was positively correlated with forward airspeed. For both samples, mean body angle and relative speed were negatively correlated (P<0.05). For the sample of 15 butterflies, no other kinematic parameters were significantly correlated with , while for the sample of 44 butterflies only ϕmax and the stroke amplitude Φ showed a significant increase with V (P<0.05).

Mean wingbeat frequency decreased with body mass to the power -0.24 for the sample of 15 butterflies, and with mass to the power —0.26 for the sample of 44 butterflies (Fig. 2); both of these exponents were significantly less than 0 (P<0.05). For neither sample, however, was there a significant correlation between wingbeat frequency and wing loading (P>0.05).

Fig. 2.

Wingbeat frequency n (in Hz) vs body mass m (in kg) for 44 butterflies. The line shown is the reduced major axis correlation (r=0.52, P<0.001); the equation is given by: n=1.194m−0.258.

Fig. 2.

Wingbeat frequency n (in Hz) vs body mass m (in kg) for 44 butterflies. The line shown is the reduced major axis correlation (r=0.52, P<0.001); the equation is given by: n=1.194m−0.258.

Many butterflies glided frequently and, of the species examined here, 60% glided to some extent (Table 4). Betts and Wootton (1988) noted that, in papilionid butterflies, the fore- and hindwings are generally separated during gliding. Because of the video camera orientation used here, it was not possible to evaluate in detail the separation between the fore- and hindwings in gliding utterflies. In most species, however, this separation did not appear to be substantial and, particularly in the Pieridae and in the satyrine nymphalid Pierella luna, fore- and hindwings appeared to overlap during gliding flight. Overlap of the hindwings onto the body may improve aerodynamic flow and reduce body drag in either flapping or gliding flight. However, no significant correlation was found between the relative overlap of the wings on the body (=sbody/S) and either forward airspeed or the extent of gliding (P>0.05).

Because of the low filming frequency of the video camera, detailed analysis of wing orientation and deformation at the ends of half-strokes was not possible. Wing rotation at the ends of half-strokes was not, however, particularly pronounced in most flight sequences examined. Particularly for those butterflies with a significant overlap of the hindwing onto the body, the wing appeared to be maintained at an angle relative to horizontal close to that of the mean body angle throughout the wingbeat. Wing deformation during each half-stroke appeared to be minimal. At the beginning of the wingbeat, a pronounced clap and subsequent peeling apart of the two wing pairs was apparent in many butterflies (see Ellington, 1984a). This was most evident in the Pieridae (e.g. Aphrissa boisduvalii).

Butterflies are far from morphologically uniform, as even a cursory examination of Tables 1-3 will reveal. Masses of the butterflies examined here range over two orders of magnitude, and there is considerable variation in wing and body shape parameters. Nonetheless, wing loading and aspect ratio are, in general, lower in butterflies than in other insects. Thus, wing areas relative to body mass are higher and wings are generally broader. Correspondingly, relative wing mass is substantially higher in butterflies than in other insect orders. Wing lengths typically exceed body lengths by 54 %, mean body diameters of butterflies average 12 % of the body length, and the center of mass of the body is located in the generally elongated abdomen. With the exception of mean body diameter and wing loading, morphological parameters increase isometrically with body mass. Wing loading in butterflies exhibits positive allometry, which has also been found by Norberg and Rayner (1987) for bats over a wide mass range.

In the sample of 50 butterflies analyzed for morphology, allometric relationships between the non-dimensional radii of wing area, wing mass and wing virtual mass are, in general, the same as those found by Ellington (1984b). However, there are exceptions, particularly the correlation between and ,which also differs from a linear regression derived from the area radii of Betts and Wootton (1988) for 11 butterflies. For the butterflies examined here, wing area is typically much more skewed distally [higher ; see Table 2] than in the insects examined by Ellington (1984b) and Betts and Wootton (1988). The butterflies examined by Betts and Wootton (1988) were primarily papilionids, whose broad wing bases tend to reduce the radius of the centroid of wing area. However, the correlation relating to does not differ significantly from that of Ellington (1984b) or from that of Betts and Wootton (1988). Also, no significant correlations were found relating and to , and such correlations are not significant for the butterfly mass and area radii given in Betts and Wootton (1988). These results suggest that, while wing mass distributions may be conservative across the Insecta, wing area distributions are not necessarily so.

Paralleling the morphological diversity of butterflies is substantial diversity in wing and body kinematics (Table 4). Several general trends are nonetheless apparent. Body angles during free flight are close to horizontal (averaging 16°), and stroke plane angles are correspondingly high (averaging 67°). Through the course of a wingbeat, body angles tend to vary, at times substantially, and on average decrease by 14° from top to bottom of the half-stroke. Stroke amplitudes are not particularly high, typically being close to 100°, although maximum positional angles are large (average 69°). Wingbeat frequencies of butterflies are low, averaging here around 11 Hz. In the insectary, butterflies flew at airspeeds generally in excess of 1ms−1, and the advance ratio of the wings averaged 0.93. The magnitude of the forward velocity vector was thus comparable to that of the flapping velocity vector. In contrast, the advance ratios for two bumblebee workers flying at 1 ms−1 were 0.13 and 0.16, and the effects of flapping dominated flow over the wings (Dudley and Ellington, 1990).

Allometric relationships between forward airspeed and various morphological parameters generally corresponded to expectations based on isometric scaling (Table 5). Correlations were strongest when multiple flight sequences were averaged for each individual. For the sample of 15 butterflies each with three analyzed flight sequences, several correlations deviated significantly from the expected relationship based on isometry, most notably the variation of airspeed with body mass and wing loading. For flying animals in general, higher airspeeds are correlated with increased wing loadings. Chai and Srygley (1990) also found for 52 butterfly species in Costa Rica that airspeed was positively correlated with wing loading. No significant correlation was found here between flight speed and aspect ratio. In bats, Norberg and Rayner (1987) found that field airspeeds increased with aspect ratio, but suggested that this was an indirect effect arising from a positive correlation of aspect ratio with wing loading.

No systematic variation in wing kinematics with absolute flight speed is apparent in butterflies, although in some cases maximum positional angles and stroke amplitudes do show a significant increase with relative flight airspeed. Mean body angles decrease significantly with relative, but not with absolute, flight speeds. Betts and Wootton (1988) found for individual papilionid butterflies that slower flight was characterized by an increase in the stroke amplitude, stroke plane angle and upstroke: downstroke ratio, and by a decrease in the wingbeat frequency. Chai and Srygley (1990) suggested that fast-flying butterflies have deeper wingbeats, but no correlation was found here between stroke amplitude or minimum positional angle and flight speed. The highest relative flight speed observed here was attained by a riodinid (Thisbe irenea), which flew at 140 Ls−1. Its mean body angle was 7° and the stroke plane angle 65°; stroke amplitude was fairly high at 123°. The fastest (in non-dimensional terms) butterfly studied by Betts and Wootton (1988) was Precis iphita (Nymphalidae), which flew at a low body angle and high stroke plane angle, and exhibited the highest stroke amplitude and upstroke: downstroke ratio of those butterflies studięd. In absolute terms, the fastest butterfly studied here was Archaeoprepona demophon (Nymphalidae: Charaxinae), which in one sequence flew at 3.2 ms−1. At this speed, the body angle was 8° and the stroke plane angle was 85°. Stroke amplitude (111°) was, however, only slightly greater than the sample average (103°). Interestingly, wing aspect ratio in A. demophon was below the sample average. Betts and Wootton (1988) suggested that low aspect ratio wings would be associated with slow flight but, as mentioned above, no such relationship was found here.

Filming in an enclosed insectary raises the possibility that flight behavior is not normal or is modified in some way unknown to the experimenter. This proposition is difficult to evaluate, given the absence of knowledge concerning flight behavior under natural conditions. There is little evidence to suggest, however, that the insectary dimensions seriously constrained flight activity or speeds. Butterflies cruised continuously in the insectary for minutes on end. The highest airspeeds measured in the insectary are substantial; as mentioned previously, the fastest flier (A. demophon) attained a speed of 3.2ms−1. The butterfly Caligo illioneus (Nymphalidae: Brassolinae) is of comparable mass (>800mg) to A. demophon, and in the insectary was on average slower by 54 % (N=2 for each species). Field observations confirm that cruising flight in this species is much slower than in A. demophon. In the same subfamily as A. demophon is the mimetic butterfly Consul fabius, which unlike other charaxines has a fairly slow cruising flight (DeVries, 1987). The one filmed individual of this species flew in the insectary much more slowly than did A. demophon (average of 1.6ms−1 vs 2.9ms−1). While these comparisons are suggestive, only direct measurements of airspeeds for butterflies in natural free flight will evaluate the validity of studying flight behavior in an enclosed space.

Individual butterflies possess a considerable diversity of flight modes (sensuBetts and Wootton, 1988), ranging from hovering (at flowers or during mate selection and oviposition) to fast forward flight. Hovering in butterflies appears to be much more poorly controlled than in many dipterans and hymenopterans, with considerable lateral and vertical oscillation of the body. Slow deliberate flight characterizes butterflies in search of oviposition or nectar sites, while startled or migrating butterflies are probably much faster (see also Nikolaev, 1974). Here, the cruising flight of butterflies in the insectary corresponds to the fast forward flight mode described by Betts and Wootton (1988). Advance ratios of this flight mode are high, averaging 0.9 (Table 4). Because butterflies are capable of such variation in flight modes and performance, it is correspondingly difficult to determine a ‘typical’ flight speed for a given species. Additionally, intraspecific differences may exist in flight speed between the sexes, and in wing and body morphology depending upon nutritional state and age. By confining analysis of different butterflies to typical cruising flight in the same insectary and under approximately equivalent microclimatic conditions, flight behavior is operationally standardized.

Even then, analysis of multiple flight sequences for any one individual is appropriate. The strong correlations between flight speed and morphology which emerge here suggest that at least some degree of realism and standardization in flight behavior has been attained.

Phil DeVries and Bob Srygley provided valuable advice during the course of this study. I thank the Smithsonian Tropical Research Institute for use of superb field, sites and laboratories on Barro Colorado Island. This work was supported by a Smithsonian Institution Postdoctoral Fellowship.

Ackery
,
P. R.
(
1984
).
Systematic and faunistic studies on butterflies
.
Symp. R. entomol. Soc., Land
.
11
,
9
21
.
Betts
,
C. R.
and
Wootton
,
R. J.
(
1988
).
Wing shape and flight behaviour in butterflies (Lepidoptera: Papilionoidea and Hesperioidea): a preliminary analysis
.
J. exp. Biol
.
138
,
271
288
.
Chai
,
P.
(
1986
).
Field observations and feeding experiments on the responses of rufous-tailed jacamars (Galbula ruficauda) to free-flying butterflies in a tropical rainforest
.
Biol. J. Linn. Soc
.
29
,
161
189
.
Chai
,
P.
and
Srygley
,
R. B.
(
1990
).
Predation and flight, morphology, and temperature of Neotropical rainforest butterflies
.
Am. Nat.
(in press)
.
Devries
,
P. J.
(
1987
).
The Butterflies of Costa Rica
.
Princeton
:
Princeton University Press
.
Dudley
,
R.
and
Ellington
,
C. P.
(
1990
).
Mechanics of forward flight in bumblebees. I. Kinematics and morphology
.
J. exp. Biol
.
148
,
19
52
.
Edmunds
,
M.
(
1974
).
Defence in Animals
.
Harlow
:
Longman Group Limited
.
Ellington
,
C. P.
(
1984a
).
The aerodynamics of hovering insect flight. I. The quasi-steady analysis
.
Phil. Trans. R. Soc. Ser. B
305
,
1
15
.
Ellington
,
C. P.
(
1984b
).
The aerodynamics of hovering insect flight. II. Morphological parameters
.
Phil. Trans. R. Soc. Ser. B
305
,
17
40
.
Ellington
,
C. P.
(
1984C
).
The aerodynamics of hovering insect flight. HI. Kinematics
.
Phil. Trans. R. Soc. Ser. B
305
,
41
78
.
Humphries
,
D. A.
and
Driver
,
P. M.
(
1970
).
Protean defence by prey animals
.
Oecologia
5
,
285
302
.
Huntington
,
E. I.
(
1932
).
A list of the Rhopalocera of Barro Colorado Island, Canal Zone, Panama
.
Bull. Am. Mus. nat. Hist
.
63
,
191
230
.
Nikolaev
,
N. A.
(
1974
).
Speeds and types of flights of some cabbage butterflies/Pieridae and Nymphalidae (Lepidoptera: Pieridae, Nymphalidae)
.
Bull. Moscow Univ
.
4
,
17
20
.
(in Russian)
.
Norberg
,
U. M.
(
1981
).
Allometry of bat wings and legs and comparison with bird wings
.
Phil. Trans. R. Soc. Ser. B
292
,
359
398
.
Norberg
,
U. M.
(
1986
).
Wing form and flight mode in bats
.
In Recent Advances in the Study of Bats
(ed.
M. B.
Fenton
,
P.
Racey
and
J. M. V.
Rayner
), pp.
43
56
.
Cambridge
:
Cambridge University Press
.
Norberg
,
U. M.
and
Rayner
,
J. M. V.
(
1987
).
Ecological morphology and flight in bats (Mammalia: Chiroptera): wing adaptations, flight performance, foraging strategy and echolocation
.
Phil. Trans. R. Soc. Ser. B
316
,
335
427
.
Pennycuick
,
C. J.
(
1975
).
Mechanics of flight
.
In Avian Biology
(ed.
D. S.
Farner
and
J. R.
King
), pp.
1
75
.
London
:
Academic Press
.
Rayner
,
J. M. V.
(
1985
).
Linear relations in biomechanics: the statistics of scaling functions
.
J. Zool., Lond
.
206
,
415
439
.
Rayner
,
J. M. V.
(
1987
).
Form and function in avian flight
.
In Current Ornithology
, vol.
5
(ed.
R. F.
Johnston
), pp.
1
65
.
New York
:
Plenum Press
.
Viscor
,
G.
and
Fuster
,
J. F.
(
1987
).
Relationships between morphological parameters in birds with different flying habits
.
Comp. Biochem. Physiol
.
87A
,
231
249
.