ABSTRACT
The flight pattern of many fly species consists of straight flight segments interspersed with rapid turns called body saccades, a strategy that is thought to minimize motion blur. We analyzed the body saccades of fruit flies (Drosophila hydei), using high-speed 3D videography to track body and wing kinematics and a dynamically scaled robot to study the production of aerodynamic forces and moments. Although the size, degree and speed of the saccades vary, the dynamics of the maneuver are remarkably stereotypic. In executing a body saccade, flies perform a quick roll and counter-roll, combined with a slower unidirectional rotation around their yaw axis. Flies regulate the size of the turn by adjusting the magnitude of torque that they produce about these control axes, while maintaining the orientation of the rotational axes in the body frame constant. In this way, body saccades are different from escape responses in the same species, in which the roll and pitch component of banking is varied to adjust turn angle. Our analysis of the wing kinematics and aerodynamics showed that flies control aerodynamic torques during the saccade primarily by adjusting the timing and amount of span-wise wing rotation.
INTRODUCTION
When exploring a local environment, many species of fly exhibit a flight pattern consisting of straight segments interspersed with rapid turns called body saccades (Collett and Land, 1975; Tammero and Dickinson, 2002; Wagner, 1986; Wehrhahn et al., 1982). This distinct flight pattern might serve many functions, but one likely advantage is that the quick turns allow flies to restrict the time periods during which their visual system is severely compromised by motion blur (Collett and Land, 1975; Hateren and Schilstra, 1999; Land, 1999; Schilstra and Hateren, 1999). Whatever the function, it appears to be the preferred means of changing direction in some species. For example, a recent analysis of Drosophila melanogaster found that more than 80% of all changes in heading occur via body saccades (van Breugel and Dickinson, 2012).
There is some controversy regarding the neural mechanisms that trigger saccades in Drosophila (Dickinson, 2014). Several studies of both free and tethered flight behavior suggest that most saccades are triggered by visual expansion and thus represent collision-avoidance reflexes that protect flies from flying into large obstacles or avoiding clutter (Censi et al., 2013; Reiser and Dickinson, 2013; Stewart et al., 2010; Tammero and Dickinson, 2002). Other studies, however, suggest that some saccades are triggered internally by a deliberately stochastic process that functions to optimize the animal's search efficiency (Maye et al., 2007; Reynolds and Frye, 2007). In addition, flies exhibit rapid turns in other contexts, such as when they lose contact with an odor plume (van Breugel and Dickinson, 2014), or in response to a rapid visual expansion as might be created by an approaching predator (Muijres et al., 2014). It is not known, however, whether rapid turns elicited by different stimuli or internal triggers operate via a single common motor program.
The aerodynamic basis of saccades, which is the main subject of this paper, has been investigated previously by Fry and co-workers using the fruit fly Drosophila melanogaster (Fry et al., 2003). Based on a relatively small number of high-speed video sequences, these authors proposed a model in which flies change course primarily by creating torque around their yaw axis (defined in that study as perpendicular to the longitudinal body axis), which they accomplish by simultaneously changing stroke amplitude and deviations out of the stroke plane. In addition, these authors reported that the changes in wing kinematics were biphasic, which they interpreted as indicating that flies generate first torque and then counter-torque during each maneuver. Finally, they presented a simple model in which the dynamics about the yaw axis were dominated by inertia during the brief maneuvers, consistent with the production of counter-torque (Fry et al., 2003). Subsequent authors, however, challenged some conclusions of this simple model. Hesselberg and Lehmann (2007) noted that because of the reciprocating pattern of wing motion, the damping about the yaw axis is quite large and should quickly dominate dynamics during turns – a calculation that was supported by subsequent models and measurements (Cheng et al., 2010; Dickson et al., 2010). Hedrick and co-workers (2009) went so far as to suggest that animals over a very large size range need only produce a small amount of counter-torque during saccades, and can rely primarily on passive damping to coast to a stop after initiating a turn. Further, studies of corrective maneuvers in Drosophila suggested that flies generate yaw torque by regulating the angle of attack of the wing during the upstroke and downstroke, and not by altering either stroke amplitude or stroke deviation (Bergou et al., 2010).
In this paper, we employed 3D high-speed videography to capture the wing and body motion of the fruit fly Drosophila hydei Sturtevant 1921 during free-flight body saccades. Although this study repeats the basic paradigm developed by Fry and co-workers over a decade ago (Fry et al., 2003), the improvements in high-speed cameras as well as the utility of an automated machine vision-tracking system allowed us to revisit saccade dynamics with greater resolution and statistical rigor. The results demonstrate that body saccades are a remarkably stereotyped behavior combining a brief banked turn (requiring rotation and counter-rotation in roll and pitch) with a unidirectional rotation about the yaw axis. By measuring the changes in wing motion during saccades and using a dynamically scaled robot, we were able to determine the relative importance of different features of wing motion in generating forces and torques. The results help to resolve some of the recent controversies regarding the dynamics of saccadic turns in Drosophila and other insects.
- A
amplitude
- a={ax, ay, az}
acceleration vector in the world reference frame
- ahor
horizontal acceleration
- an
Fourier series coefficient
- aN
horizontal acceleration normal to the flight path
- aT
horizontal acceleration tangential to the flight path
- az
vertical acceleration
- bn
Fourier series coefficient
- C
normalized aerodynamic damping coefficient matrix
- dt
time step between measurements (inverse of camera frame rate)
- dU
flight speed relative to the start of the maneuver
- f
wingbeat frequency
- F
normalized aerodynamic force vector in the world reference frame
- g
gravitational acceleration scalar
- g={0,0,g}
gravitational acceleration vector
- I
normalized inertia coefficient matrix
- l
wing length
- m
body mass
- mg
body weight
- Q
process noise covariance matrix of Kalman filter
- R
measurement noise covariance matrix of Kalman filter
- t
time relative to the start of the saccade
- T
normalized aerodynamic torque vector in the body reference frame
- T′
an arbitrarily defined aerodynamic torque vector
- t*
normalized time within a wingbeat relative to the start of the downstroke
- U={u,v,w}
velocity vector in the world reference frame
- X={x,y,z}
position vector in the world reference frame
- XB={xB,yB,zB}
position vector in the body reference frame
- α
rotation angle of the wing around its long axis
- β
pitch angle of the aerodynamic force vector in the body reference frame
- γ
deviation angle of the wing out of the stroke plane
- Δt
duration of a saccadic turn
- Δα
shift of the mean wing rotation angle of the left wing relative to the right
- Δκ
shift of the mean wing kinematics angle of the left wing relative to the right
- Δσ
turn angle of the saccade
- Δτ
phase shift of the left wing movement relative to the right wing movement
- η
body roll angle derived from roll rate ωx
- ηE
body roll Euler angle in the world reference frame
- θ
body pitch angle derived from pitch rate ωy
- θE
body pitch Euler angle in the world reference frame
- κ
wing kinematics angle
- μ
torque axis angle within the stroke plane relative to the roll torque axis
- μ1
primary torque axis angle (mean torque angle during the initial phase of the banked turn)
- μ2
counter-torque axis angle (mean torque angle during the counter-torque phase of the banked turn)
- μ1⊥
angle within the stroke plane and orthogonal to µ1
- μ2⊥
angle within the stroke plane and orthogonal to µ2
- ξ
roll angle of the aerodynamic force vector in the body reference frame
- σ
heading (the direction of the horizontal body velocity component)
- σa
direction of the horizontal body acceleration component
- τ
phase within a wingbeat relative to the start of the downstroke
- φ
stroke angle of the wing within the stroke plane
- ψ
body yaw angle derived from yaw rate ωz
- ψE
body yaw Euler angle in the world reference frame
- Ω={ωx, ωy, ωz}
normalized rotation rate vector of the body in the body reference frame
- ωx
normalized body roll rate in the body reference frame
- ωy
normalized body pitch rate in the body reference frame
- ωz
normalized body yaw rate in the body reference frame
RESULTS
We tracked a total of 44 flight sequences (Fig. 1, see Materials and methods), each consisting of a straight flight segment followed by a single body saccade (see supplementary material Movies 1 and 2). Assuming that there is no difference between left- and right-hand turns, we mirrored all left-hand turns, and then aligned all sequences based on time and heading with respect to the start of the saccade (Fig. 1E,F; see also supplementary material Movies 3–6). Heading is defined as the angular direction of the flight path, not the body orientation. The saccadic turn angles (Δσ), which quantify the total angular change in heading, varied substantially among saccades from approximately 20 deg to almost 180 deg (Fig. 1E,F and Fig. 2A,B), with an average of 93±27 deg (mean±s.d., N=44). During the maneuvers, flight speed tended to dip briefly and then gradually rise (Fig. 2C). The duration of turn (Δt=tstop−tstart) was 49±18 ms (N=44), or approximately nine wingbeats, although this is an underestimate of the entire maneuver because tstart and tstop were defined using finite thresholds (see Materials and methods). These basic saccade metrics are similar to those previously reported on free-flying D. melanogaster (Tammero and Dickinson, 2002; Fry et al., 2003; van Breugel and Dickinson, 2012).
To examine how flies alter heading during saccades, we measured the magnitude and direction of horizontal and vertical accelerations throughout the maneuver (Fig. 2D–F). Immediately at the start of the saccade, flies generate a horizontal force that results in a sideways acceleration. The magnitude (ahor/g) of this sideways acceleration first increases and then decreases (Fig. 2E), while its orientation (σa) remains relatively constant (Fig. 2D). Vertical acceleration (az/g) remains near zero throughout the entire maneuver (Fig. 2F).
The time course of roll rate, pitch rate and yaw rate during the maneuvers (Fig. 1C and Fig. 2G–I) along with their integrals and derivatives (supplementary material Fig. S1) show that flies rotate about all principal body axes, and thus body saccades constitute a banked turn. Flies also increase force production so the vertical component remains roughly equal to body weight (Fig. 2F,J). The orientation of F in the body reference frame remains constant (Fig. 2K,L), consistent with the so-called ‘helicopter model’ of insect flight (David, 1978; Götz and Wandel, 1984). Thus, as with more rapid escape maneuvers (Muijres et al., 2014), a fly generates sideways accelerations during saccades by rotating its body rather than by adjusting the orientation of the force vector in the body frame.
To determine how the flies control roll, pitch and yaw throughout the saccade, we estimated torque about these axes as the sum of torque required to overcome inertia (inertia torque) and torque required to overcome damping (damping torque) (Fig. 3, see Eqn 2 in Materials and methods). Torque about the yaw axis consists mostly of damping torque, confirming that yaw rotations during a saccade are highly damped (Hedrick et al., 2009; Hesselberg and Lehmann, 2007). In contrast, torque about the roll axis consists about equally of inertial torque and damping torque, whereas damping torque about the pitch axis is negligible compared with inertial torque. The fact that rotations about the pitch axis are poorly damped during saccades is supported by a recent study of forward flight dynamics in fruit flies (Elzinga et al., 2014).
A simple means of implementing a banked turn would be to rotate the body about a fixed axis in the stroke plane and then to counter-rotate to continue level flight, while at the same time generating a yaw rotation to align body orientation with the new flight heading. The magnitude of the turn could then be adjusted by regulating the amount of torque produced, and not its direction. Evidence that flies might implement such a simple control scheme is shown in Fig. 4. The torque vectors for the primary rotation phase and the counter-rotation phase for all 44 sequences (in blue and orange, respectively) aligned remarkably well (Fig. 4A,B). The average torque vector axes for the initial rotation and counter-rotation are defined as the primary torque axis µ1 and the counter-torque axis µ2, which are oriented 36 and 8 deg from the longitudinal body axis, respectively. The time history of the direction (Fig. 4A) and magnitude (Fig. 4C,D) of the torque component in the stroke plane exhibits a biphasic shape indicative of rotation and counter-rotation about the µ1 and µ2 axis, respectively (Fig. 4C). In contrast, little torque is generated about the orthogonal axes (μ1⊥ and μ2⊥, respectively, Fig. 4D; see also supplementary material Fig. S2). Thus, animals execute the banked turn by generating torque and counter-torque about two control axes (µ1 and µ2) whose orientation remains constant from saccade to saccade (Fig. 4A,B). The magnitude of the torque produced about the control axes, however, does correlate with the turn angle Δσ (Fig. 4F), which suggests the mechanism by which flies regulate the size of the heading change.
Simultaneously with the rotations about the µ1 and µ2 axes, but within a longer time scale, flies generate a unidirectional yaw torque in the same direction as the change in heading (Fig. 4E). Also, the magnitude of mean yaw torque during the turn is positively correlated with turn angle (Fig. 4H). Although not necessary for changing the direction of the flight path, the yaw rotation is required to align the longitudinal body axis with the new heading. This correction to minimize sideslip is not fully completed in most sequences because of limitations of our visualization volume (Fig. 4E). Note that by rotating about the yaw axis while its body is banked, the fly will generate a head-down movement within the world frame of reference. This could explain why the initial body rotation axis µ1 and counter-rotation axis µ2 are not aligned and µ2 includes a smaller pitch (down) component.
The sequence of 21 averaged wingbeats (N=44 trials) show that all kinematics parameters (wingbeat frequency, stroke angle, deviation angle and wing rotation angle, Fig. 1C) change during a saccade, and that these modifications are all very subtle (Fig. 5A–G). Wingbeat frequency increases by only a few hertz and all modifications in wing angles are less than 5 deg. Nevertheless, replaying the averaged kinematics sequence on the robotic fly generated normalized forces (|F|) and torques (Troll, Tpitch and Tyaw) that were similar to those estimated from body dynamics (using Eqns 1 and 2 in Materials and methods, respectively; Fig. 5H–K). For example, the time history for the torque about the roll axis measured on the robotic fly exhibits the torque and counter-torque that is predicted from body dynamics (Fig. 5I). Thus, although the measured changes in wing motion are subtle, they appear sufficient to capture the requisite changes in forces and moments reasonably well. Because the robotic fly is fixed and cannot translate and rotate in response to the forces and moments it generates, we did not expect (nor obtain) a perfect match between the measured forces and moments and those calculated from body dynamics.
To better understand how flies modulate wingbeat kinematics to control torque around each individual control axis during a saccade, we extracted the wing kinematics associated with peak torque about the µ1, µ2 and yaw axes from the entire dataset and replayed these kinematics on the robotic fly (Fig. 6, see Materials and methods). The kinematics correlated with peak torque production exhibit distortion of all three wing angles (Fig. 6A–C) and, when played through the robot, produced the expected torques (Fig. 6D–F). The fly's wing motion is able to create positive yaw torque throughout almost the entire wingbeat, with the exception of brief periods during stroke transitions (Fig. 6D). In contrast, torque production about the µ1 and µ2 axes is more complicated in that the time history includes both positive and negative excursions and the magnitude of the transient peaks is quite large relative to the average value (Fig. 6E). Most of the variations in torque production relative to the torque generated by the steady flight wingbeat (i.e. the symmetric wingbeat that produced weight support and no net torque, Fig. 1D), occur at the start of the upstroke and downstroke, just after stroke reversal (Fig. 6E). The wingbeat patterns that generate peak torque about the µ1 and µ2 axes create very little mean torque about the orthogonal axes (μ1⊥ and μ2⊥), even though the magnitude of the oscillations through the stroke is quite large (Fig. 6F).
Next, we constructed a set of systematically distorted stroke patterns ranging from steady flight conditions to kinematics that generate peak torque about the µ1, µ2 and yaw axes, and then replayed these on the robot (see Materials and methods). For yaw, torque measured using the robot matched the values derived from body motion (Eqn 1) throughout the entire range of distorted kinematics (average difference ∼2%, Fig. 7A). For the µ1 and µ2 axes, the torques derived from body motion are about 66% and 59% (respectively) that measured with the robotic fly, suggesting that the dynamic model for roll and pitch (Eqns 2–4) may be oversimplified or that the torque measured using a robot in a fixed reference frame does not accurately model the free-flight case. We also measured the cross-talk (i.e. the torque generated around orthogonal axes) produced by the kinematics associated with production of peak torque around the different torque axes. For the yaw case, the cross-talk about the roll and pitch axes was −24% and 27% of |Tyaw|, respectively (Fig. 7A). For the µ1 case, the cross-talk about the μ1⊥ and yaw axes was 3% and 45% of |Tµ1|, respectively, and for the µ2 case, the cross-talk around the μ2⊥ and yaw axes was 23% and 47% of |Tµ2|, respectively (Fig. 7B). The cross-talk between the yaw torque and torque about the axes in the stroke plane (Fig. 7A,B) shows that our method did not enable us to completely separate the effect of torque and force production about the different principal axes. This is most likely due to the fact that during saccadic maneuvers, flies increase force production and produce torque about the different axes in synchrony. So, wingbeats that produced large torque about the axes in the stroke plane tended to also produce high yaw torque and increased aerodynamic forces (Fig. 5).
Next, we varied the kinematics for each of the three wing angles in isolation to determine their relative contribution to torque (Fig. 7C,D). As was also the case for evasive maneuvers (Muijres et al., 2014), the sum of torques resulting from modulating the different components of wing motion separately matches the torque generated by modulating all components simultaneously, indicating a remarkable degree of linearity. For all torque axes, wing rotation angle had the strongest effect on changes in total torque, whereas changes in stroke amplitude and stroke deviation contribute modestly to the control of torque around the µ1 and µ2 axes and make almost no contribution to the control of yaw torque (Fig. 7C,D).
As suggested by a previous study of D. melanogaster, a change in the mean offset of the rotation angle will create yaw torque by increasing the angle of attack during one half-stroke and decreasing it on the other (Bergou et al., 2010). However, another means by which wing rotation might influence torque is through changes in the timing of wing rotation relative to stroke reversal, which can change forces via unsteady mechanisms (Dickinson et al., 1999). To explore the relative importance of these two potential mechanisms, we estimated both the relative phase shift and mean offset of the time course of the wing kinematics angles for the left and right wings during the strokes that generated peak torque (Fig. 8; supplementary material Fig. S3). Fig. 8C–F shows that there is a phase shift of ∼5 deg in the wing rotation angle of the left wing relative to the right wing during wing strokes that produce peak torque around the µ1 and µ2 axes, but there is no evidence for a change in mean offset. In the case of yaw torque, a phase shift of ∼3 deg is accompanied by an offset of ∼4 deg. Thus, although torque about the yaw, µ1 and µ2 axes is primarily controlled by changes in the time course of wing rotation, our results suggest that the relative mechanisms are different. Flies regulate the torque about the µ1 and µ2 axes by modulating unsteady rotational lift mechanisms during stroke reversal (Dickinson et al., 1999), whereas yaw torque is controlled by a combination of unsteady effects at stroke reversal and differences in drag during the translational phase of the two half-strokes (Bergou et al., 2010). This interpretation is consistent with the time history of the changes in torque throughout the wingbeat for the kinematics that produce peak torques (Fig. 6). Yaw torque is produced mostly during the translation phase of the wingbeat (Fig. 6D), which is indicative of an asymmetry in drag production between the two wings. In contrast, changes in torque about the µ1 and µ2 axes (relative to a steady wing stroke) occur mostly at the start of each wingstroke (Fig. 6E), suggesting modulations in rotational effects.
DISCUSSION
Our analysis showed that fruit flies perform body saccades by executing stereotyped banked turns (Fig. 2). The axes of the initial rotation (µ1) and subsequent counter-rotation (µ2) are aligned 36 and 8 deg from the roll axis of the fly, respectively (Fig. 4B). Flies control the size of the turn by regulating the magnitude of torque around these rotation axes (Fig. 4F), and not by adjusting their orientation as they do during more rapid escape maneuvers (Muijres et al., 2014). Flies also rotate unidirectionally around the yaw axis during saccades to correct for the misalignment between body orientation and heading (i.e. sideslip) that accumulates as a result of the banked turn.
Using a simple dynamic model, we were able to estimate the relative contribution of inertia and damping during a saccade (Fig. 3). We found that yaw dynamics are dominated by passive damping that results from the reciprocal flapping pattern, as suggested by Hesselberg and Lehmann (2007). Pitch dynamics, in contrast, are dominated by inertia, and the contribution of damping and inertia in roll dynamics is roughly equal. Thus, saccade dynamics are quite complex and accurate models must include both inertial and damping terms (Bergou et al., 2010; Cheng et al., 2010; Dickson et al., 2010; Fry et al., 2003; Hedrick et al., 2009; Hesselberg and Lehmann, 2007). The model for rapid turns proposed by Hedrick and co-workers (2009) in which insects rely primarily on passive damping to generate counter-torque may not be as general as proposed, assuming that other insects also employ banked turns and do not simply rotate about the yaw axis.
The maneuver that we have measured in fruit flies using high-speed videography is similar to the body saccades of blowflies, measured elegantly by Schilstra and Hateren (1999) using tiny inductive coils. Thus, at least two species of fly, encompassing a rather large range in scale, execute banked turns to produce rapid changes in flight direction. Given the high degree of similarity within their nervous systems, it is likely that the body saccades of Drosophila and Calliphora are generated by homologous circuits.
Visually elicited escape maneuvers (Muijres et al., 2014) and voluntary body saccades in fruit flies both consist of a banked turn, but the dynamics of the two maneuvers are different enough to suggest that they are produced by distinct motor programs. Although body saccades are fast, the changes in heading during an evasive maneuver are faster, consistent with a more pronounced rotation of the body that reorients the mean force vector to produce a larger horizontal component. As a consequence, flies do not maintain weight support during the initial stages of an escape maneuver, whereas they do during body saccades (Fig. 2F,J). Similarly, yaw is poorly controlled during the initial stages of an evasive maneuver, resulting in large sideslip angles that are corrected long after the fly changes heading. During a body saccade, the yaw correction is better coordinated with the banked turn, so that sideslip is minimized throughout the maneuver. Perhaps the greatest difference between the two maneuvers relates to the manner by which the magnitude of the change in heading is controlled. Flies regulate the turn angle of a saccade by varying torque magnitude about two highly stereotypic axes (µ1 and µ2), whereas during evasive maneuvers the turn angle is controlled by adjusting the direction of the body rotation axis within the stroke plane. One interpretation that unifies all these differences is that evasive maneuvers may be optimized to alter flight heading as quickly as possible at the expense of flight control and motion blur, whereas body saccades are optimized to restrict retinal slip to a brief period.
The transition from rotation to counter-rotation during a body saccade is quite fast and such a pattern might be generated in a feed-forward manner by a central motor program or, alternatively, the initial rotation might trigger the counter-rotation via a sensory-mediated reflex. The halteres are a likely source of such feedback, as they mediate compensatory reactions to imposed rotations (Dickinson, 1999; Nalbach, 1994; Ristroph et al., 2010; Sherman and Dickinson, 2003). Previous studies of saccades using a magnetic tether in which the animal is free to rotate about its yaw axis suggest that haltere feedback, but not visual feedback, plays a role in terminating saccades (Bender and Dickinson, 2006). Interpretation of these prior experiments using magnetic tethers is complicated, however, by the new free-flight data, which show that Drosophila bank to change direction at the start of the saccade. Given their time course, the slow unidirectional saccade-like rotations that flies exhibit on magnetic tethers most likely represent the slower yaw phase of a free-flight saccade. Thus, the manipulation experiments performed by Bender and Dickinson (2006) might indicate that haltere feedback is involved in regulating the duration of the slower yaw correction phase of a saccade but do not directly address the question of whether feedback triggers the faster counter-rotation phase of the initial banked turn.
By exploiting the high-throughput capabilities of our tracking system, we were able to make accurate measurements of the changes in wing motion during saccades (Fig. 6). All three wing angles (φ, γ, α) exhibit a biphasic modulation during the time course of the saccade, as expected from the production of torque and counter-torque (Fig. 5B–G). These data are consistent with the previous observations of Fry and co-workers (2003), although these authors misinterpreted this biphasic pattern as indicating active breaking around the yaw axis, when it is more likely that they observed evidence for the counter-rotation about the µ2 axis. Replaying the pattern of wing motion on a stationary robot generated a time history of forces and moments that matched those derived from free-flight body dynamics reasonably well, but not perfectly (Fig. 5H–K). The match for total flight force and yaw was quite good, whereas roll torque based on the wing kinematics and the robot measurements was larger than roll torque derived from body dynamics. This mismatch is perhaps not too surprising, given that the saccade basically consists of a rapid roll and counter-roll (i.e. the µ1 and µ2 axes are not too far off the roll axis). Thus, errors that derive from the fact that the kinematics from a rotating fly were replayed on a stationary robot would be particularly large for this degree of freedom. Other sources of error include the possible inaccuracies of our dynamics model, which did not include cross-terms and used damping coefficients based on steady-state approximations.
By mining the entire database, we were able to determine the pattern of wing motion that correlated with peak torque production around the µ1, µ2 and yaw axes (Fig. 6), and then to determine the relative contribution of the three wing angles to the moments by playing systematically distorted wing patterns through the robotic fly (Fig. 7). The results indicate that changes in the time history of the wing rotation angle, which strongly influences the angle of attack, are by far the most important for regulating torque about the µ1, µ2 and yaw axes (Fig. 7C,D). However, the torque modulations about the yaw axis and about the two axes in the stroke plane (µ1 and µ2) appear to occur via two distinct aerodynamic mechanisms (Fig. 8). As suggested by Bergou and co-workers (2010), changes in the mean offset of the wing rotation angle magnitude lower the angle of attack during one half-stroke while raising it on the other, thus producing an upstroke-to-downstroke imbalance in drag and thus net torque around the axis normal to the stroke plane. Although our results support this mechanism (Fig. 6A,D and Fig. 7A,B), we also found that the flies adjust the relative phase of wing rotation as well, thus creating additional yaw torque via unsteady mechanisms at stroke reversal (Dickinson et al., 1999; Sane and Dickinson, 2002). Flies also create torque about the µ1 and µ2 axes via changes in wing rotation angle, but in these cases the effect appears to be mediated almost entirely by changes in phase, and thus through unsteady effects at stroke reversal (Fig. 6E and Fig. 7C–F). One possible explanation for this difference is that changes in angle of attack during the translational portion of the stroke would be undesirable for regulating roll because they would generate large cross-talk in yaw torque.
A quite surprising result of our analysis was the relatively small importance of stroke amplitude and stroke deviation in the control of torque during saccades (Fig. 7C,D). This was particularly true for yaw torque, an observation that complicates interpretation of many tethered flight studies, which collectively show that flies generate large changes in stroke amplitude in response to both visual and mechanosensory rotations about the yaw axis (e.g. Sherman and Dickinson, 2003), as well as transient spontaneous changes that have been ubiquitously interpreted as fictive saccades. These changes in stroke amplitude were quite large and clearly correlated with yaw torque (Tammero, 2004). Why do tethered flies generate such large changes in stroke amplitude that do not seem necessary to generate yaw torque in free flight? One possibility is that the stroke amplitude signal measured during fictive saccades is indicative of the roll and pitch required for a banking maneuver, and not for yaw production per se. If true, this has immediate implications for the underlying circuitry as it is noteworthy that the stroke amplitude changes associated with fictive saccades are unidirectional, i.e. there is no evidence of an attempt at a programmed counter-rotation.
Recently, Schnell and co-workers (2014) suggested that the basic optomotor circuit in Drosophila includes an integral feedback term that might be mediated by Ca2+ dynamics in the terminals of the interneurons that encode horizontal rotation. Because the putative integrator winds up with prolonged stimulation, flies generate extremely large motor responses during typical open-loop optomotor experiments. Although this hypothesis explains why the optomotor responses in tethered flight are so large relative to the kinematics changes exhibited during free-flight maneuvers, it does not explain why flies generate such enormous changes in stroke amplitude when presented with stimuli that should elicit the production of yaw torque. Perhaps flies respond to visual rotation about the yaw axis by producing not just yaw torque but also roll torque, because they are attempting to turn via banking, as has been suggested by previous authors (Mronz and Lehmann, 2008). This hypothesis could be tested in the future by measuring free-flight responses to horizontal motion.
MATERIALS AND METHODS
Animals and experimental setup
The methods used in this study were identical to those described in detail elsewhere (Muijres et al., 2014) and are only briefly outlined here. Experiments were performed on 1–5 day old male and female D. hydei, from a laboratory stock reared on a 14 h:10 h light:dark cycle. Each day, approximately 50, 1 day old flies were released in the experimental chamber 4 h before their subjective dawn, after which experiments ran for 8 h. The experimental chamber consisted of a transparent cylindrical enclosure, surrounded by a panoramic array of green LED panels (Reiser and Dickinson, 2008) that provided a uniform illumination of 70 lx. Flies were filmed using three synchronized high-speed cameras (Photron SA5 with AF Nikor 60 mm lenses, lens aperture=f/22), which viewed the central portion of the arena from above and from two orthogonal side positions (Fig. 1A). The cameras operated at 7500 frames s−1 with an image resolution of 1024×1000 pixels, exposure time of 1/30,000 s, and image depth of 12 bits. Because of limitations in optics, the region of interest was restricted to a cube ∼40 mm on each side. Each camera view was backlit using high-intensity infrared light panels, which were strobed in synchrony with every exposure. To maintain the inside temperature at ∼25°C, we passed refrigerated air around the outside of the flight chamber. At the start and end of each 5 day recording period, we calibrated the camera system using direct linear transformation (DLT) (Christoph Reinschmidt; http://isbweb.org/software/movanal.html).
Throughout each experimental session, the cameras sampled continuously. When a fly flew directly through the center of the region of interest, it tripped an infrared laser cross-beam, which automatically triggered the storage of 372 ms of data both before and after the trigger event. We captured more than 300 flight sequences, of which the majority consisted of a straight flight path, but occasionally a fly would perform a saccadic maneuver (supplementary material Movie 1). Although fruit flies are largely insensitive to infrared light, we were concerned that the high intensity of the infrared trigger lasers might elicit behavioral responses. For this reason, we only analyzed the saccades that flies initiated before passing through the trigger point. Given our sampling methods, we have no way of knowing whether any given saccade was elicited by visual expansion or via some internal stochastic event.
Measuring body and wingbeat kinematics
We manually selected 44 sequences from the entire data set for detailed analysis, which we subjectively classified as body saccades based on visual inspection of the raw video sequences. To extract kinematics throughout these saccades, we used an automatic machine-vision system as described in a recent analysis of escape maneuvers (Fig. 1B) (Muijres et al., 2014). The tracking routine provided us with Kalman-filtered estimates of the body and wing kinematics throughout each flight sequence (Fig. 1C; supplementary material Movie S2) (Muijres et al., 2014).
For all Kalman filters, the cross-product values in the error covariance matrices (R and Q) were set to zero, and measurement noise covariance matrix R was set to identity. Thus, adjusting the parameters in the process noise covariance matrix Q controlled Kalman smoothing. For positional data, we used a linear Kalman filter and the smoothing parameters of matrix Q were scaled according to a Taylor series. Q values corresponding to position were set to dt2, velocity values were set to 1, and acceleration values were set to 1/dt2, where dt is the time step between two measurements (inverse of the frame rate).
To filter body and wing orientation data, we used an extended Kalman filter with two filtering steps (Yun and Bachmann, 2006), because quaternion update is non-linear. For body orientation filtering, we used the following Q parameters. For the first iteration, Q parameters associated with angular velocities and angular accelerations were set to zero, and parameters associated with quaternions were set to 0.003. For the second iteration, angular velocity parameters were set to 0.0001, angular acceleration parameters were set to 0.0001/dt2, and quaternion parameters were set to 1. For wing orientation filtering, we used a similar strategy but less smoothing was desired. For the first iteration, angular velocity parameters and angular acceleration parameters were set to zero, and quaternion parameters in Q were set to 1. For the second iteration, angular velocity parameters in Q were set to 1/dt2, angular acceleration parameters were set to 0 and quaternion parameters were set to 1.
Body position data thus consist of Kalman-filtered estimates of position, X(t), velocity, U(t), and linear acceleration, a(t). Based on U(t) and a(t), we determined the start and end of each saccade by estimating horizontal accelerations normal and tangential to the flight path [aN(t) and aT(t), respectively]. Using an expectation maximization-based clustering analysis on aN(t) and aT(t) for all measured flight sequences (Muijres et al., 2014), we divided flight sequences into steady and maneuvering segments, such that steady segments satisfied aN<0.19 g and −0.17 g <aT<0.14 g, where g is gravitational acceleration. The times at which a saccade starts (tstart) and ends (tend) were defined as the moments when body acceleration passed these threshold values. Note that these points in time do not encompass the entire length of the maneuver, but the method enabled us to systematically and objectively define and align the saccades.
Based on the velocity data, we determined flight heading, σ(t), which we used to determine the turn angle, Δσ, for each saccade. U(t) was also used to determine changes in flight speed, dU(t), throughout the maneuver relative to the start of each sequence. Linear acceleration, a(t), was used to estimate the direction, σa, and magnitude, ahor/g, of normalized acceleration within the horizontal plane, as well as the magnitude of the vertical acceleration, az/g. Body orientation was expressed in the body Euler angles (yaw, ψE; pitch, θE; and roll, ηE) in the world reference frame and normalized angular velocity of the body, Ω={ωx, ωy, ωz}, in the body reference frame (Fig. 1C). All rotation rates were normalized using fsteady, where fsteady is the wingbeat frequency during steady flight for D. hydei (∼189 Hz; from Muijres et al., 2014).
Each tracked wingbeat was expressed by the wingbeat frequency, f, and three Euler angles within the body reference frame: wing stroke angle φ, stroke deviation angle γ and wing rotation angle α (Fig. 1C). Note that all body and wing orientation variables are defined relative to the stroke plane, which is defined as horizontal during steady flight (Fig. 1C). For D. hydei, this is at an inclination angle of 47.5 deg relative to the long axis of the body (Muijres et al., 2014).
Estimating aerodynamic forces and torques based on body dynamics
where F is the aerodynamic force vector normalized with body weight, mg, and g={0,0,g} is the gravitational acceleration vector. The orientation of this force vector in the body frame was defined as the roll (ξ) and pitch (β) angle of F relative to the stroke plane normal (Fig. 1C), and was calculated based on the body Euler angles and F.
Thus, our simplified model assumes that all interaction coefficients in both C and I are negligible.
The damping coefficients in C were based on damping estimates for D. melanogaster reported in the literature (Cheng et al., 2009; Dickson et al., 2010). The yaw-damping coefficient for D. hydei (Cyaw) was estimated by scaling values measured for D. melanogaster using a robotic model (Dickson et al., 2010). Roll and pitch damping for D. hydei were based on computational estimates of damping coefficients in D. melanogaster (Cheng et al., 2009). For consistency among the damping coefficients for the three degrees of freedom, we first linearly scaled all damping coefficients estimated for D. melanogaster by Cheng et al. (2009) equally, such that yaw damping was equal to that reported by Dickson et al. (2010), and then scaled these values to the slightly larger species, D. hydei. This resulted in the following normalized roll, pitch and yaw damping coefficients: Croll=0.22, Cpitch=0.08, Cyaw=0.41. Note that the damping coefficient for yaw is 5 times greater than that for pitch and twice as large as that for roll.
Inertia coefficients within the stroke plane reference frame were estimated based on a cylindrical body model with body mass m and pitch angle of 47.5 deg and a wing model consisting of a horizontal disk divided into 100 concentric rings (Fig. 1C). Each ring has a homogeneously distributed mass equal to the mass of the local spanwise wing section plus its added fluid mass. Following Ellington (1984), total wing mass was estimated as 5% of body mass, and added mass was equal to a cylindrical fluid mass circumventing each wing section. This resulted in normalized roll inertia Iroll=0.64, normalized pitch inertia Ipitch=1.07 and normalized yaw inertia Iyaw=0.57.
Correlating wingbeat kinematics with torque production
where mod(κ,T′)i is the wingbeat modification variable for kinematics angle κ (representing φ, γ or α) of the ith wingbeat in the dataset of |T′| producing wingbeats. κsteady is the equivalent kinematics angle distribution of the average steady wingbeat of D. hydei, based on 1603 wingbeats measured by Muijres et al. (2014) (Fig. 1D).
Measuring aerodynamic forces and torques using a dynamically scaled robot
Apart from estimating aerodynamic forces and torques from body dynamics (Eqns 1 and 2, respectively), we also estimated forces and torques from wing kinematics, using a dynamically scaled robot in a fixed body reference frame (Dickinson et al., 1999). This technique enabled us to study the aerodynamics of saccadic maneuvers in a systematic and detailed manner. We replayed the wingbeat kinematics in a fixed body reference frame because aerodynamic and inertial effects of body rotations were modeled using Eqn 2. Note that aerodynamics effects due to body translations were ignored.
To assess the accuracy of our methods, we compared torque measurements based on the wingbeat kinematics with the torque estimates based on body dynamics (Eqn 2). This analysis was performed on two sets of data. The first dataset consisted of the mean wing kinematics and body dynamics throughout the average saccade, estimated by aligning the wingbeats from all maneuvers. This enabled us to test qualitatively whether the measured changes in wingbeat kinematics captured the forces and torques produced throughout the saccade. The second dataset was based on the systematic analysis described above, in which we correlated wingbeat kinematics with torque production. Using Eqn 7, we constructed a set of kinematics patterns that should produce a systematically increasing amount of body torque about a specific body axis T′. The range of body torques was chosen such that it captured the complete behavioral envelope of measured torque production. We defined a body torque distribution spanning a range from zero body torque (steady flight) to a torque equal to approximately three times the standard deviation of the |T′| distribution for all measured wingbeats, and parsed it into eight values. For each, we created the matching deformed wingbeat kinematics using Eqn 7. We then replayed the set of systematically distorted wing kinematics on the robotic fly and measured the resulting stroke-averaged forces and torques about all three orthogonal body axes, and compared these values with the values of body torques. This approach also enabled us to determine cross-talk between torque modulations about the different orthogonal axes. The distributions of |T′| and the corresponding wing kinematics constructed using Eqn 7 were also used to study the effects of stroke, deviation and wing rotation angle on torque production. Using a method similar to that described by Muijres et al. (2014), we systematically modulated one wing kinematics variable, while maintaining steady kinematics for the other angles, and replayed these on the robot.
Cases in which body torque was equal to three standard deviations of the |T′| distribution, which we considered as an estimate for peak torque production, were analyzed in more detail using two complementary methods. First, we replayed the wingbeat kinematics that corresponded to this peak torque (estimated using Eqn 7) on the robot and measured force and torque throughout the wingbeat. Second, we analyzed the changes in wing kinematics that result in peak torque production. For each kinematics angle, the difference in kinematics between the left and right wing might result from temporal phase shift, from a shift in the mean value, or from higher order modulations in the time history (Fig. 1D). For each kinematics variable, we estimated the temporal phase shift and offset in mean value between the left and right wing by systematically translating the left wing data along the time (τ) and ordinate (κ) axes. For each combination of Δτ and Δκ, we determined the root mean square error (RMSE) between the right wing and shifted left wing values. The combination of Δτ and Δκ at which RMSE was minimum defines the phase shift and angular offset between left and right wing for that wing kinematics angle and |T′| axis combination. The corresponding magnitude of RMSE quantifies how well the wing kinematics modulations are described by Δτ and Δκ, relative to any higher order modulations.
Acknowledgements
The authors wish to thank Johan Melis, Steve Safarik and Darren Howell for help with this project.
Author contributions
F.T.M. and M.H.D. planned experiments and wrote the paper. F.T.M. and N.A.I. collected data in flight arena, ran the automated machine vision tracking software and identified sequences with saccadic turns. M.J.E. collected data using the robotic fly. F.T.M. analyzed the data, with the help of M.J.E. and M.H.D.
Funding
This research was supported by grants (to M.H.D.) from the Air Force Office of Scientific Research (FA9550-10-1-0368), the Paul G. Allen Family Foundation, and the US Army Research Laboratory (DAAD 19-03-D-0004) and (to F.T.M.) the Swedish Research Council (Vetenskapsrådet).
References
Competing interests
The authors declare no competing or financial interests.