# Bumblebees minimize control challenges by combining active and passive modes in unsteady winds

Sridhar Ravi, Dmitry Kolomenskiy, Thomas Engels, Kai Schneider, Chun

Wang, Jörn Sesterhenn, Hao Liu

► To cite this version:

Sridhar Ravi, Dmitry Kolomenskiy, Thomas Engels, Kai Schneider, Chun Wang, et al.. Bumblebees minimize control challenges by combining active and passive modes in unsteady winds. Scientific Reports, 2016, 6 (35043), 10.1038/srep35043 . hal-01485154

HAL Id: hal-01485154

https://hal.science/hal-01485154

Submitted on 2 Jun 2024

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Bumblebees minimize control challenges by combining active and passive modes in unsteady winds

Sridhar Ravi, Dmitry Kolomenskiy, Thomas Engels, Kai Schneider, Chun Wang, Jörn Sesterhenn, Hao Liu

► To cite this version:

Sridhar Ravi, Dmitry Kolomenskiy, Thomas Engels, Kai Schneider, Chun Wang, et al., Bumblebees minimize control challenges by combining active and passive modes in unsteady winds. Scientific Reports, 2016, 6 (35043), 10.1038/srep35043. doi: 10.14851/54

HAL Id: hal-01485154

https://hal.science/hal-01485154v1

Submitted on 2 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Bumblebees minimize control challenges by combining active and passive modes in unsteady winds

Received: 04 April 2016

Accepted: 01 September 2016

Published: 18 October 2016

Sridhar Ravi<fnref n="1" />,²*, Dmitry Kolomenskiy<fnref n="1" />*, Thomas Engels³, Kai Schneider⁴, Chun Wang⁵, Jörn Sesterhenn³ & Hao Liu<fnref n="1" />,⁶

The natural wind environment that volant insects encounter is unsteady and highly complex, posing significant flight-control and stability challenges. It is critical to understand the strategies insects employ to safely navigate in natural environments. We combined experiments on free flying bumblebees with high-fidelity numerical simulations and lower-order modeling to identify the mechanics that mediate insect flight in unsteady winds. We trained bumblebees to fly upwind towards an artificial flower in a wind tunnel under steady wind and in a von Kármán street formed in the wake of a cylinder. Analysis revealed that at lower frequencies in both steady and unsteady winds the bees mediated lateral movement with body roll - typical casting motion. Numerical simulations of a bumblebee in similar conditions permitted the separation of the passive and active components of the flight trajectories. Consequently, we derived simple mathematical models that describe these two motion components. Comparison between the free-flying live and modeled bees revealed a novel mechanism that enables bees to passively ride out high-frequency perturbations while performing active maneuvers at lower frequencies. The capacity of maintaining stability by combining passive and active modes at different timescales provides a viable means for animals and machines to tackle the challenges posed by complex airflows.

Volant insects are able to maintain stable flight in natural environments, despite frequently encountering unpredictable, complex airflows. While laboratory studies of insect flight typically employ steady flow<fnref n="1" />, natural flows are often unsteady, containing de-stabilizing vortices and gusts over a wide range of spatial and temporal scales. These unsteady flows are generated by the interaction between wind and obstacles²,³, which results in chains of discrete vortices in the near wake of trees, branches and flowers, and fully developed, turbulent flow farther away.

The fact that unsteady airflow is ubiquitous in natural environments raises the fundamental question of whether insects need to actively respond to the full spectrum of unsteady flows that they experience in the natural world, or whether some of these disturbances can be handled with minimal energetic cost via passive mechanisms. Several recent studies have quantified the flight trajectories and kinematics of insects subjected to various flow perturbations, such as discrete gusts⁴, von Kármán vortex streets⁵,⁶, or steady vortices⁷. However, kinematics alone do not allow one to easily distinguish between passive body motions induced by a fluid disturbance and active, voluntary motions generated by the insect to maneuver or correct for a perturbation, particularly when the flow environment is unpredictable or unknown.

In this study, we focus on disentangling the active vs. passive components of a bumblebee's motion while flying through the unsteady, but structured, von Kármán vortex street generated behind a circular cylinder in flow. Bumblebees are ideal model organisms for studying insect flight in complex airflow, because they are relentless foragers that continue to fly in a wide range of weather conditions⁸,⁹. Prior work has shown that many flying animals, including bumblebees, are most sensitive to aerial disturbances along the lateral axis and rotations along the body's longitudinal axis (i.e. roll)⁴,⁶. Rolling motion and lateral translation are also critical components of

<footnote n="1">¹Graduate School of Engineering, Chiba University, Japan. ²School of Aerospace Mechanical and Manufacturing Engineering, RMIT University, Melbourne Australia. ³ISTA, Technische Universität Berlin, Müller-Breslau-Strasse 12, 10623 Berlin, Germany. ⁴I2M - UMR 7373 - CNRS, Centre de Mathématiques et d'Informatique, Aix-Marseille Université, France. ⁵School of Mechanical and Manufacturing Engineering, University of New South Wales, New South Wales, Australia. ⁶Shanghai-Jiao Tong University and Chiba University International Cooperative Research Centre (SJTU-CU ICRC), Shanghai, People's Republic of China. *These authors contributed equally to this work. Correspondence and requests for materials should be addressed to S.R. (email: dkolom@chiba-u.jp)</footnote>
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Bumblebees minimize control challenges by combining active and passive modes in unsteady winds

Received: 04 April 2016

Accepted: 01 September 2016

Published: 18 October 2016

Sridhar Ravi<fnref n="1" />,²*, Dmitry Kolomenskiy<fnref n="1" />*, Thomas Engels³, Kai Schneider⁴, Chun Wang⁵, Jörn Sesterhenn³ & Hao Liu<fnref n="1" />,⁶

The natural wind environment that volant insects encounter is unsteady and highly complex, posing significant flight-control and stability challenges. It is critical to understand the strategies insects employ to safely navigate in natural environments. We combined experiments on free flying bumblebees with high-fidelity numerical simulations and lower-order modeling to identify the mechanics that mediate insect flight in unsteady winds. We trained bumblebees to fly upwind towards an artificial flower in a wind tunnel under steady wind and in a von Kármán street formed in the wake of a cylinder. Analysis revealed that at lower frequencies in both steady and unsteady winds the bees mediated lateral movement with body roll - typical casting motion. Numerical simulations of a bumblebee in similar conditions permitted the separation of the passive and active components of the flight trajectories. Consequently, we derived simple mathematical models that describe these two motion components. Comparison between the free-flying live and modeled bees revealed a novel mechanism that enables bees to passively ride out high-frequency perturbations while performing active maneuvers at lower frequencies. The capacity of maintaining stability by combining passive and active modes at different timescales provides a viable means for animals and machines to tackle the challenges posed by complex airflows.

Volant insects are able to maintain stable flight in natural environments, despite frequently encountering unpredictable, complex airflows. While laboratory studies of insect flight typically employ steady flow<fnref n="1" />, natural flows are often unsteady, containing de-stabilizing vortices and gusts over a wide range of spatial and temporal scales. These unsteady flows are generated by the interaction between wind and obstacles²,³, which results in chains of discrete vortices in the near wake of trees, branches and flowers, and fully developed, turbulent flow farther away.

The fact that unsteady airflow is ubiquitous in natural environments raises the fundamental question of whether insects need to actively respond to the full spectrum of unsteady flows that they experience in the natural world, or whether some of these disturbances can be handled with minimal energetic cost via passive mechanisms. Several recent studies have quantified the flight trajectories and kinematics of insects subjected to various flow perturbations, such as discrete gusts⁴, von Kármán vortex streets⁵,⁶, or steady vortices⁷. However, kinematics alone do not allow one to easily distinguish between passive body motions induced by a fluid disturbance and active, voluntary motions generated by the insect to maneuver or correct for a perturbation, particularly when the flow environment is unpredictable or unknown.

In this study, we focus on disentangling the active vs. passive components of a bumblebee's motion while flying through the unsteady, but structured, von Kármán vortex street generated behind a circular cylinder in flow. Bumblebees are ideal model organisms for studying insect flight in complex airflow, because they are relentless foragers that continue to fly in a wide range of weather conditions⁸,⁹. Prior work has shown that many flying animals, including bumblebees, are most sensitive to aerial disturbances along the lateral axis and rotations along the body's longitudinal axis (i.e. roll)⁴,⁶. Rolling motion and lateral translation are also critical components of

<footnote n="1">¹Graduate School of Engineering, Chiba University, Japan. ²School of Aerospace Mechanical and Manufacturing Engineering, RMIT University, Melbourne Australia. ³ISTA, Technische Universität Berlin, Müller-Breslau-Strasse 12, 10623 Berlin, Germany. ⁴I2M - UMR 7373 - CNRS, Centre de Mathématiques et d'Informatique, Aix-Marseille Université, France. ⁵School of Mechanical and Manufacturing Engineering, University of New South Wales, New South Wales, Australia. ⁶Shanghai-Jiao Tong University and Chiba University International Cooperative Research Centre (SJTU-CU ICRC), Shanghai, People's Republic of China. *These authors contributed equally to this work. Correspondence and requests for materials should be addressed to S.R. (email: dkolom@chiba-u.jp)</footnote>
voluntary maneuvering, and some insects have been shown to employ roll as the preferred means to perform rapid flight maneuvers¹⁰,¹¹. We therefore chose to generate a von Kármán vortex street behind a vertical cylinder, which provides the greatest lateral and rolling perturbations⁶, and chose a cylinder size ($d = 25$ mm) corresponding to a typical bumblebee's wing span, which produces a chain of vortices shed at a rate of 23 Hz⁶. We compare the rolling and lateral flight dynamics of bumblebees flying in the vortex street to their flight in steady flow. The relative importance of active control mechanisms vs. passive dynamics would most readily be tested by varying the level of active control that an insect can employ as it flies through unsteady flow, and quantifying the resulting change in its body dynamics and flight trajectory. However, such manipulations are unfeasible in living specimens, but numerical simulations of flying insects offer a suitable alternative. Simulations that combine flapping wing aerodynamics and free-flight dynamics provide a powerful tool that allows for simultaneous manipulation of all kinematic parameters of the insect model, aerodynamic forces and torques, and the ambient flow field. These types of simulations have provided significant insight into both the fine scale flow structures that form over flapping wings and the overall dynamics of insect flight, see refs 12 and 13 for review, but such simulations have only been performed in still air hitherto. Here, we develop a model of a flapping bumblebee flying through a von Kármán vortex street that is nominally identical to the one utilized in experiments on real bees. The model accounts for flow properties at all physically relevant length scales in the oncoming flow and resolves the resulting lateral and roll flight dynamics of the insect due to its interaction with the wind. Such high fidelity simulations require extensive computational resources, and over 8192 cores and >2.2 billion grid points were used in this study, the highest number that has been used to resolve flapping flight to date.

We compare the flight trajectories and body motions of live bumblebees flying in the wake of a cylinder to those of the model bee in similar flow conditions with no active control, and propose simple analytical models, to parsimoniously estimate the role of active control and identify the control strategies being implemented by flying bumblebees.

## Results

**Bumblebee Flight Dynamics.** During free flight in both steady and unsteady flows, the flight path of the bees typically consisted of large-amplitude motions at low frequencies (i.e., “casting”) that combined lateral translation and rolling (Fig. 1a)⁶. In unsteady flow, higher frequency translational and rolling motions were superimposed on top of these low-frequency casting motions (Fig. 1c, Supplementary Video 1), and earlier investigations revealed that these high-frequency motions occur primarily at the 23-Hz which is the vortex shedding frequency of the cylinder⁶. In contrast, the lateral casting movements typically occur at frequencies below 10 Hz, and this was used as the cut-off frequency for examining the low- vs. high-frequency components of motion. The overall results were found to be nominally insensitive to the choice of cut-off frequency, so long as it was smaller than the frequency of the von Kármán street.

The magnitude of lateral accelerations associated with low-frequency, casting motions was similar in both steady and unsteady flow ($p = 0.12$, Fig. 2a). High-frequency motions (>10 Hz) were not observed in steady flow. However in unsteady flow, the lateral accelerations associated with high-frequency motions were significantly higher than those associated with casting motions ($p < 10^{-7}$, Fig. 2a). The mean absolute roll angles of bees were also similar in both airflow conditions during casting ($p = 0.48$, Fig. 2b). The absolute roll angles measured at higher frequencies in unsteady flow were significantly smaller than those at lower frequencies ($p < 10^{-6}$, Fig. 2b).

During low-frequency casting, among the flight trails of all individuals, there was a strong, positive correlation between the instantaneous roll angle of bees and their lateral acceleration ($C = 0.89 \pm 0.16$ in steady flow, $C = 0.81 \pm 0.16$ in unsteady flow), and the strength of this correlation was similar in steady and unsteady flow ($p = 0.20$). However, for higher-frequency motions in unsteady flow, there was a strong, negative correlation between roll angle and lateral acceleration ($C = -0.81 \pm 0.13$), which was significantly different from the correlation for low-frequency motion in unsteady flow ($p < 10^{-17}$, Fig. 2c).

Due to the high cross correlations, to estimate the relationship between roll angle and lateral acceleration, we performed a linear regression between these variables for each flight trial of all individuals, for low-frequency casting motions in both flow conditions and for high-frequency motions in unsteady flow (Fig. 1b,d,e). We found the slope of the regression across all flight trials and compared mean values between flow conditions. For low-frequency casting motions, a statistically similar, positive regression slope was present in both wind conditions ($p = 0.866$, Fig. 2d). However, at higher frequencies in unsteady flow, the slope of this relationship was negative, and significantly different from the slope for low-frequency motions in unsteady flow ($p < 10^{-14}$, Fig. 2d).

**Assessing the Presence of Active Flight Control.** The high-fidelity numerical simulations were conducted to evaluate the passive flight dynamics of a flapping bumblebee with no active control, flying in a von Kármán vortex street with the same characteristics as the experimental unsteady flow condition. Flight was simulated over 0.25 sec (38 wing beats), during which time the bumblebee interacted with a number of vortices in the von Kármán street (Fig. 3, Supp. Video 2, Supp. Fig. 1). The main outcome of the numerical simulation is the time evolution of the aerodynamic forces and torques acting on the insect, as well as the resulting lateral and rolling motions of the ‘numerical’ bee. Because the simulations did not include voluntary low-frequency casting motions, the flight data of the model bee in the simulations were high-pass filtered and analyzed in the same way as the high-frequency motions of real bees flying in unsteady flow.

The high-pass filtered lateral and rotational positions and accelerations from the simulation showed good agreement with experimental measurements (Fig. 4a–d). The time-varying lateral and rolling motions were nominally within the standard deviations of the high-frequency components of motion measured on freely flying bees in unsteady wind, and mean values of absolute roll angle and lateral acceleration were not significantly different from measurements on real bees (roll angle: $p = 0.65$; lateral acceleration: $p = 0.37$).
voluntary maneuvering, and some insects have been shown to employ roll as the preferred means to perform rapid flight maneuvers¹⁰,¹¹. We therefore chose to generate a von Kármán vortex street behind a vertical cylinder, which provides the greatest lateral and rolling perturbations⁶, and chose a cylinder size ($d = 25$ mm) corresponding to a typical bumblebee's wing span, which produces a chain of vortices shed at a rate of 23 Hz⁶. We compare the rolling and lateral flight dynamics of bumblebees flying in the vortex street to their flight in steady flow. The relative importance of active control mechanisms vs. passive dynamics would most readily be tested by varying the level of active control that an insect can employ as it flies through unsteady flow, and quantifying the resulting change in its body dynamics and flight trajectory. However, such manipulations are unfeasible in living specimens, but numerical simulations of flying insects offer a suitable alternative. Simulations that combine flapping wing aerodynamics and free-flight dynamics provide a powerful tool that allows for simultaneous manipulation of all kinematic parameters of the insect model, aerodynamic forces and torques, and the ambient flow field. These types of simulations have provided significant insight into both the fine scale flow structures that form over flapping wings and the overall dynamics of insect flight, see refs 12 and 13 for review, but such simulations have only been performed in still air hitherto. Here, we develop a model of a flapping bumblebee flying through a von Kármán vortex street that is nominally identical to the one utilized in experiments on real bees. The model accounts for flow properties at all physically relevant length scales in the oncoming flow and resolves the resulting lateral and roll flight dynamics of the insect due to its interaction with the wind. Such high fidelity simulations require extensive computational resources, and over 8192 cores and >2.2 billion grid points were used in this study, the highest number that has been used to resolve flapping flight to date.

We compare the flight trajectories and body motions of live bumblebees flying in the wake of a cylinder to those of the model bee in similar flow conditions with no active control, and propose simple analytical models, to parsimoniously estimate the role of active control and identify the control strategies being implemented by flying bumblebees.

## Results

**Bumblebee Flight Dynamics.** During free flight in both steady and unsteady flows, the flight path of the bees typically consisted of large-amplitude motions at low frequencies (i.e., “casting”) that combined lateral translation and rolling (Fig. 1a)⁶. In unsteady flow, higher frequency translational and rolling motions were superimposed on top of these low-frequency casting motions (Fig. 1c, Supplementary Video 1), and earlier investigations revealed that these high-frequency motions occur primarily at the 23-Hz which is the vortex shedding frequency of the cylinder⁶. In contrast, the lateral casting movements typically occur at frequencies below 10 Hz, and this was used as the cut-off frequency for examining the low- vs. high-frequency components of motion. The overall results were found to be nominally insensitive to the choice of cut-off frequency, so long as it was smaller than the frequency of the von Kármán street.

The magnitude of lateral accelerations associated with low-frequency, casting motions was similar in both steady and unsteady flow ($p = 0.12$, Fig. 2a). High-frequency motions (>10 Hz) were not observed in steady flow. However in unsteady flow, the lateral accelerations associated with high-frequency motions were significantly higher than those associated with casting motions ($p < 10^{-7}$, Fig. 2a). The mean absolute roll angles of bees were also similar in both airflow conditions during casting ($p = 0.48$, Fig. 2b). The absolute roll angles measured at higher frequencies in unsteady flow were significantly smaller than those at lower frequencies ($p < 10^{-6}$, Fig. 2b).

During low-frequency casting, among the flight trails of all individuals, there was a strong, positive correlation between the instantaneous roll angle of bees and their lateral acceleration ($C = 0.89 \pm 0.16$ in steady flow, $C = 0.81 \pm 0.16$ in unsteady flow), and the strength of this correlation was similar in steady and unsteady flow ($p = 0.20$). However, for higher-frequency motions in unsteady flow, there was a strong, negative correlation between roll angle and lateral acceleration ($C = -0.81 \pm 0.13$), which was significantly different from the correlation for low-frequency motion in unsteady flow ($p < 10^{-17}$, Fig. 2c).

Due to the high cross correlations, to estimate the relationship between roll angle and lateral acceleration, we performed a linear regression between these variables for each flight trial of all individuals, for low-frequency casting motions in both flow conditions and for high-frequency motions in unsteady flow (Fig. 1b,d,e). We found the slope of the regression across all flight trials and compared mean values between flow conditions. For low-frequency casting motions, a statistically similar, positive regression slope was present in both wind conditions ($p = 0.866$, Fig. 2d). However, at higher frequencies in unsteady flow, the slope of this relationship was negative, and significantly different from the slope for low-frequency motions in unsteady flow ($p < 10^{-14}$, Fig. 2d).

**Assessing the Presence of Active Flight Control.** The high-fidelity numerical simulations were conducted to evaluate the passive flight dynamics of a flapping bumblebee with no active control, flying in a von Kármán vortex street with the same characteristics as the experimental unsteady flow condition. Flight was simulated over 0.25 sec (38 wing beats), during which time the bumblebee interacted with a number of vortices in the von Kármán street (Fig. 3, Supp. Video 2, Supp. Fig. 1). The main outcome of the numerical simulation is the time evolution of the aerodynamic forces and torques acting on the insect, as well as the resulting lateral and rolling motions of the ‘numerical’ bee. Because the simulations did not include voluntary low-frequency casting motions, the flight data of the model bee in the simulations were high-pass filtered and analyzed in the same way as the high-frequency motions of real bees flying in unsteady flow.

The high-pass filtered lateral and rotational positions and accelerations from the simulation showed good agreement with experimental measurements (Fig. 4a–d). The time-varying lateral and rolling motions were nominally within the standard deviations of the high-frequency components of motion measured on freely flying bees in unsteady wind, and mean values of absolute roll angle and lateral acceleration were not significantly different from measurements on real bees (roll angle: $p = 0.65$; lateral acceleration: $p = 0.37$).
<figure><img src="image_3.png" /><figcaption>**Figure 1.** (a) Sample time series of lateral displacement and roll angle of a freely flying bumblebee in steady flow. (b) Instantaneous lateral acceleration vs. roll angle for lower-frequency (<10 Hz) motions of the flight shown in (a). (c) Sample time series of lateral displacement and roll angle for the same bee flying in the unsteady wake of a vertical cylinder. (d) Instantaneous lateral acceleration vs. roll angle for lower- (d) and higher-frequency (e) motions of the flight shown in (c). In (b,d,e) the linear regression between lateral acceleration and roll angle for the trial presented in (a,c) is shown by a dashed line.</figcaption></figure>

Discussion

Aerial locomotion at small scales is very changeling and insects utilize a variety of unsteady aerodynamic mech-
anisms to generate the necessary forces to not only stay aloft but also perform coordinated maneuvers. In the
natural environment volant organisms are also posed with wind unsteadiness. What are the added challenges
imposed by wind unsteadiness and what mechanisms do insects employ to mitigate them?

**Low-frequency motions.** In both steady and unsteady winds, bumblebees performed casting movements consisting of large, lateral displacements at low frequencies (Fig. 1a & Supp. Video 1). Such casting motions have also been observed in other insects, including wasps and honeybees¹⁴,¹⁵. While the significance of casting is not completely understood, some insects such as wasps perform casting flights to aid visual processing¹⁴. In the experiments conducted here, the high, positive cross-correlation between roll angle and lateral acceleration at the low frequencies associated with casting behavior was observed (Fig. 2c). We explain it using deterministic mathematical modelling and derive a causative relationship between the roll and the lateral motions, producing phase lock that leads to the positive correlation. In the following paragraph, we highlight the main results of the modelling. For more details, see Supplementary Material S2.1.

The “helicopter model” of flight control, in which the net aerodynamic force vector remains fixed relative to the body and body attitude is adjusted to alter the direction of motion, has been demonstrated in several species of insects and birds¹⁶,¹⁷. Based on this model, acceleration or deceleration along the longitudinal axis requires a change in body pitch, and lateral maneuvers require a body rotation around the roll axis, to reorient the aerodynamic force vector. Our data on low-frequency casting motions in bumblebees supports the helicopter model of flight control. Considering level translation only along the lateral axis, the relationship between body roll (ψ, in rad) (refer to Supp. Fig. 1 for coordinate system) and lateral acceleration based on the helicopter model is (g tan ψ) (refer S2.1 for derivation), where g is the acceleration due to gravity. The roll angles observed during casting were sufficiently small, and upon applying the small angle approximation, the relationship between the roll angle and lateral acceleration becomes linear with the proportionality coefficient nominally equal to the acceleration due to gravity (see S2.1 for derivation). In our experiments, since the bees maintained altitude while casting in both steady and unsteady flow, the slope of the regression line between body roll and lateral acceleration was close to the expected value of 9.8 m/s²/rad (steady airflow = 7.6 ± 1.2 m/s²/rad, unsteady airflow = 7.8 ± 1.5 m/s²/rad; Fig. 2d).

The high correlation and positive relationship between low-frequency rolling motions and lateral acceleration
in both steady and unsteady winds (Fig. 2c) suggests that the bees were performing the same type of controlled
maneuvers at low frequencies in both flow situations. Additionally, the presence of casting motions even in the
<figure><img src="image_3.png" /><figcaption>**Figure 1.** (a) Sample time series of lateral displacement and roll angle of a freely flying bumblebee in steady flow. (b) Instantaneous lateral acceleration vs. roll angle for lower-frequency (<10 Hz) motions of the flight shown in (a). (c) Sample time series of lateral displacement and roll angle for the same bee flying in the unsteady wake of a vertical cylinder. (d) Instantaneous lateral acceleration vs. roll angle for lower- (d) and higher-frequency (e) motions of the flight shown in (c). In (b,d,e) the linear regression between lateral acceleration and roll angle for the trial presented in (a,c) is shown by a dashed line.</figcaption></figure>

Discussion

Aerial locomotion at small scales is very changeling and insects utilize a variety of unsteady aerodynamic mech-
anisms to generate the necessary forces to not only stay aloft but also perform coordinated maneuvers. In the
natural environment volant organisms are also posed with wind unsteadiness. What are the added challenges
imposed by wind unsteadiness and what mechanisms do insects employ to mitigate them?

**Low-frequency motions.** In both steady and unsteady winds, bumblebees performed casting movements consisting of large, lateral displacements at low frequencies (Fig. 1a & Supp. Video 1). Such casting motions have also been observed in other insects, including wasps and honeybees¹⁴,¹⁵. While the significance of casting is not completely understood, some insects such as wasps perform casting flights to aid visual processing¹⁴. In the experiments conducted here, the high, positive cross-correlation between roll angle and lateral acceleration at the low frequencies associated with casting behavior was observed (Fig. 2c). We explain it using deterministic mathematical modelling and derive a causative relationship between the roll and the lateral motions, producing phase lock that leads to the positive correlation. In the following paragraph, we highlight the main results of the modelling. For more details, see Supplementary Material S2.1.

The "helicopter model" of flight control, in which the net aerodynamic force vector remains fixed relative to the body and body attitude is adjusted to alter the direction of motion, has been demonstrated in several species of insects and birds¹⁶,¹⁷. Based on this model, acceleration or deceleration along the longitudinal axis requires a change in body pitch, and lateral maneuvers require a body rotation around the roll axis, to reorient the aerodynamic force vector. Our data on low-frequency casting motions in bumblebees supports the helicopter model of flight control. Considering level translation only along the lateral axis, the relationship between body roll (ψ, in rad) (refer to Supp. Fig. 1 for coordinate system) and lateral acceleration based on the helicopter model is (g tan ψ) (refer S2.1 for derivation), where g is the acceleration due to gravity. The roll angles observed during casting were sufficiently small, and upon applying the small angle approximation, the relationship between the roll angle and lateral acceleration becomes linear with the proportionality coefficient nominally equal to the acceleration due to gravity (see S2.1 for derivation). In our experiments, since the bees maintained altitude while casting in both steady and unsteady flow, the slope of the regression line between body roll and lateral acceleration was close to the expected value of 9.8 m/s²/rad (steady airflow = 7.6 ± 1.2 m/s²/rad, unsteady airflow = 7.8 ± 1.5 m/s²/rad; Fig. 2d).

The high correlation and positive relationship between low-frequency rolling motions and lateral acceleration
in both steady and unsteady winds (Fig. 2c) suggests that the bees were performing the same type of controlled
maneuvers at low frequencies in both flow situations. Additionally, the presence of casting motions even in the
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Mean absolute (a) lateral acceleration and (b) roll angle of lower-frequency (<10 Hz) casting motions in steady and unsteady flow, and higher-frequency (>10 Hz) motions in unsteady flow. (c) Correlation coefficients (zero time-lag, normalized cross correlation) for the relationship between instantaneous lateral acceleration and roll angle for lower-frequency (<10 Hz) casting motions in steady and unsteady flow, and higher-frequency (>10 Hz) motions in unsteady flow. (d) The slope of the regression line between lateral acceleration and roll angle, in the same flow conditions as in (c). Asterisks indicate a significant difference at the p < 0.05 level, n = 14 for all conditions.</figcaption></figure>

absence of external flow perturbations (i.e., in steady flow) indicates that these motions were most likely voluntary (see Fig. 5a and Supp. Video 3 for schematic representation and animation of casting motion, respectively). Based on the regression between lateral acceleration and roll, bees are able to generate a lateral acceleration of around 1 m/s² (typical accelerations measured during casting) by rolling approximately 8 degrees towards the direction of intended motion.

**High-frequency motions.** In the unsteady flow condition, we observed high-frequency rolling motions and lateral accelerations superimposed on top of the slow casting motions, with a clear separation between these two types of oscillations in terms of their temporal scale (Fig. 1c and Supp. Video 1). These high-frequency oscillations in unsteady flow have previously been shown to occur primarily at the vortex-shedding frequency associated with the wake of the cylinder⁶. The most obvious difference between the low-frequency casting motions and the higher-frequency oscillations is in the phase relation between rolling and lateral motion. During casting in both flow conditions, we found a strong, positive correlation between roll and lateral acceleration (i.e., the bee rolls to the right and accelerates to the right), whereas at higher frequencies, there is an equally strong but negative correlation between these variables (Fig. 2c). The slope of the regression line between lateral acceleration and roll angle is negative and much larger in magnitude for the high-frequency motions (−80 ± 17 m/s²/rad; Fig. 2d); thus, mean lateral accelerations are significantly larger and mean roll angles significantly smaller for high- vs. low-frequency motions (Fig. 2). These findings suggest that the high-frequency oscillations are fundamentally different from the voluntary, low-frequency casting motions.

The high-frequency oscillations corresponding to the vortex shedding rate could be purely passive motions induced by the external airflow, or they could represent a combination of passive, externally-induced motions and active, corrective maneuvers. To distinguish between these hypotheses, we performed high-fidelity numerical simulations based on the assumption that if the high-frequency motions were purely passive, then a 'numerical' bumblebee with no flight control would exhibit similar motion dynamics as measured in the experiments.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Mean absolute (a) lateral acceleration and (b) roll angle of lower-frequency (<10 Hz) casting motions in steady and unsteady flow, and higher-frequency (>10 Hz) motions in unsteady flow. (c) Correlation coefficients (zero time-lag, normalized cross correlation) for the relationship between instantaneous lateral acceleration and roll angle for lower-frequency (<10 Hz) casting motions in steady and unsteady flow, and higher-frequency (>10 Hz) motions in unsteady flow. (d) The slope of the regression line between lateral acceleration and roll angle, in the same flow conditions as in (c). Asterisks indicate a significant difference at the p < 0.05 level, n = 14 for all conditions.</figcaption></figure>

absence of external flow perturbations (i.e., in steady flow) indicates that these motions were most likely voluntary (see Fig. 5a and Supp. Video 3 for schematic representation and animation of casting motion, respectively). Based on the regression between lateral acceleration and roll, bees are able to generate a lateral acceleration of around 1 m/s² (typical accelerations measured during casting) by rolling approximately 8 degrees towards the direction of intended motion.

**High-frequency motions.** In the unsteady flow condition, we observed high-frequency rolling motions and lateral accelerations superimposed on top of the slow casting motions, with a clear separation between these two types of oscillations in terms of their temporal scale (Fig. 1c and Supp. Video 1). These high-frequency oscillations in unsteady flow have previously been shown to occur primarily at the vortex-shedding frequency associated with the wake of the cylinder⁶. The most obvious difference between the low-frequency casting motions and the higher-frequency oscillations is in the phase relation between rolling and lateral motion. During casting in both flow conditions, we found a strong, positive correlation between roll and lateral acceleration (i.e., the bee rolls to the right and accelerates to the right), whereas at higher frequencies, there is an equally strong but negative correlation between these variables (Fig. 2c). The slope of the regression line between lateral acceleration and roll angle is negative and much larger in magnitude for the high-frequency motions (−80 ± 17 m/s²/rad; Fig. 2d); thus, mean lateral accelerations are significantly larger and mean roll angles significantly smaller for high- vs. low-frequency motions (Fig. 2). These findings suggest that the high-frequency oscillations are fundamentally different from the voluntary, low-frequency casting motions.

The high-frequency oscillations corresponding to the vortex shedding rate could be purely passive motions induced by the external airflow, or they could represent a combination of passive, externally-induced motions and active, corrective maneuvers. To distinguish between these hypotheses, we performed high-fidelity numerical simulations based on the assumption that if the high-frequency motions were purely passive, then a 'numerical' bumblebee with no flight control would exhibit similar motion dynamics as measured in the experiments.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Birds eye view of the flow field obtained with numerical simulation. Isosurfaces of dimensionless vorticity magnitude are shown as, light blue $|ω|=5$, red $|ω|=40$. Thin lines show the velocity in the horizontal plane that passes through the insect's center of gravity. The insect is colored in orange. Flow is shown at $t=246.1$ ms after the beginning of the simulation.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Time evolution of high pass filtered (hpf, >10 Hz) lateral accelerations (**a**), of roll accelerations (**b**), lateral displacements (**c**) and roll angles (**d**) from the numerical simulation and experiments on freely flying bees. Blue line shows results of the numerical simulation (for motions >10 Hz), and black line/gray shading show the mean and standard deviation of measurements from experiments in unsteady flow (motions >10 Hz). Dotted vertical lines represent wing beat periods.</figcaption></figure>

The close similarity in temporal patterns of lateral and rolling motions exhibited by real and simulated bumblebees over several cycles of flow disturbance (approximately 20 wing beats; Fig. 4a–d) lends support to the notion that the high-frequency motions are indeed predominantly due to passive interactions with vortices in the von Kármán street. This suggests that bees do not actively respond to these flow perturbations on a wingbeat-by-wingbeat basis.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Birds eye view of the flow field obtained with numerical simulation. Isosurfaces of dimensionless vorticity magnitude are shown as, light blue $|ω|=5$, red $|ω|=40$. Thin lines show the velocity in the horizontal plane that passes through the insect's center of gravity. The insect is colored in orange. Flow is shown at $t=246.1$ ms after the beginning of the simulation.</figcaption></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Time evolution of high pass filtered (hpf, >10 Hz) lateral accelerations (**a**), of roll accelerations (**b**), lateral displacements (**c**) and roll angles (**d**) from the numerical simulation and experiments on freely flying bees. Blue line shows results of the numerical simulation (for motions >10 Hz), and black line/gray shading show the mean and standard deviation of measurements from experiments in unsteady flow (motions >10 Hz). Dotted vertical lines represent wing beat periods.</figcaption></figure>

The close similarity in temporal patterns of lateral and rolling motions exhibited by real and simulated bumblebees over several cycles of flow disturbance (approximately 20 wing beats; Fig. 4a–d) lends support to the notion that the high-frequency motions are indeed predominantly due to passive interactions with vortices in the von Kármán street. This suggests that bees do not actively respond to these flow perturbations on a wingbeat-by-wingbeat basis.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Schematic representation of the relationship between lateral and roll accelerations and displacements during (a) low frequency casting maneuvers (helicopter model) and (b) while high frequency buffeting in unsteady winds (sailboat model). The red and green arrows in (a,b) respectively represent the self-initiated and wind-induced forces and torques respectively. Here $f_c$ (<10 Hz) and $f_v$ (≈23 Hz) are casting and von Kármán frequencies respectively and t is time. The schematic drawings exaggerate the amplitudes of all motions for the purpose of clearly showing the phase relationship between roll and lateral displacement. For in-scale visualization of the casting and buffeting motions, see Supplementary Videos S3 and S4, respectively. In low frequency casting (a), the lateral force is produced by tilt of the stroke plane, self-initiated by the bee. Therefore, lateral acceleration is of the same sign as the roll angle. In high frequency buffeting (b), the wind induces drag-based force and torque. Therefore, the lateral acceleration is of the same sign as the torque, i.e., opposite to the roll angle. See Supplementary Material S2 for mathematical derivation of the phase relationship between roll and lateral motions.</figcaption></figure>

To explain the dynamics of the passive, high-frequency motions that bees exhibit in unsteady flow, we propose an idealized “sailboat” model. It is based on the cause and effect relationship between forcing due to lateral wing, and motion of the bee, respectively, see Supplementary Material S2.2. The governing equations for the high-frequency motions exhibited by the insect are derived by considering a solid body representing a bumblebee subjected to a periodic lateral wind disturbance. This is justified since the von Kármán street induces similar lateral disturbances on a bee in flight. The relationship between the lateral acceleration and roll angle experienced would be as follows (see S2.2 for derivation):

$$ \Pi = \frac{F_y}{m\psi} = \frac{-4\pi^2 f_v^2 I_{xx}}{\bar{L}m} \quad (1) $$

Here, $\Pi$ is the coefficient of proportionality between lateral acceleration and roll angle ($\psi$, in rad), $m$ is the mass of the bee, $F_y$ is the lateral acceleration, $I_{xx}$ is the roll moment of inertia, $\bar{L}$ is the moment arm, or the distance between the center of area and center of mass of the bee, and $f_v$ is the frequency of the von Kármán street. The model asserts that the force induced by a lateral wind on the wings (which are located above the center of mass for much of the stroke cycle), and the misalignment between the center of pressure and center of mass of the body, leads to a lateral displacement and a corresponding roll, akin to a sailboat rolling when subjected to a side wind (see Fig. 5b and Supplementary Video 4). This model preserves the phase relationship between roll angle and lateral acceleration while buffeting in the experiments. By substituting parameters for the “mean” bumblebee from the experiments into Eq. 1 the proportionality between the lateral acceleration and the roll angle due to the wind disturbance is around -77 m/s²/rad (see S2.2 for elaboration), which is comparable to that measured in the free flying bees (-74 ± 13 m/s²/rad), see Fig. 2d. In Eq. 1 the distance between the center of aerodynamic pressure and center of mass critically influences the relationship, thus the dynamics of the high frequency motion may differ vastly between insects.
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Schematic representation of the relationship between lateral and roll accelerations and displacements during (a) low frequency casting maneuvers (helicopter model) and (b) while high frequency buffeting in unsteady winds (sailboat model). The red and green arrows in (a,b) respectively represent the self-initiated and wind-induced forces and torques respectively. Here $f_c$ (<10 Hz) and $f_v$ (≈23 Hz) are casting and von Kármán frequencies respectively and t is time. The schematic drawings exaggerate the amplitudes of all motions for the purpose of clearly showing the phase relationship between roll and lateral displacement. For in-scale visualization of the casting and buffeting motions, see Supplementary Videos S3 and S4, respectively. In low frequency casting (a), the lateral force is produced by tilt of the stroke plane, self-initiated by the bee. Therefore, lateral acceleration is of the same sign as the roll angle. In high frequency buffeting (b), the wind induces drag-based force and torque. Therefore, the lateral acceleration is of the same sign as the torque, i.e., opposite to the roll angle. See Supplementary Material S2 for mathematical derivation of the phase relationship between roll and lateral motions.</figcaption></figure>

To explain the dynamics of the passive, high-frequency motions that bees exhibit in unsteady flow, we propose an idealized “sailboat” model. It is based on the cause and effect relationship between forcing due to lateral wing, and motion of the bee, respectively, see Supplementary Material S2.2. The governing equations for the high-frequency motions exhibited by the insect are derived by considering a solid body representing a bumblebee subjected to a periodic lateral wind disturbance. This is justified since the von Kármán street induces similar lateral disturbances on a bee in flight. The relationship between the lateral acceleration and roll angle experienced would be as follows (see S2.2 for derivation):

$$ \Pi = \frac{F_y}{m\psi} = \frac{-4\pi^2 f_v^2 I_{xx}}{\bar{L}m} \quad (1) $$

Here, $\Pi$ is the coefficient of proportionality between lateral acceleration and roll angle ($\psi$, in rad), $m$ is the mass of the bee, $F_y$ is the lateral acceleration, $I_{xx}$ is the roll moment of inertia, $\bar{L}$ is the moment arm, or the distance between the center of area and center of mass of the bee, and $f_v$ is the frequency of the von Kármán street. The model asserts that the force induced by a lateral wind on the wings (which are located above the center of mass for much of the stroke cycle), and the misalignment between the center of pressure and center of mass of the body, leads to a lateral displacement and a corresponding roll, akin to a sailboat rolling when subjected to a side wind (see Fig. 5b and Supplementary Video 4). This model preserves the phase relationship between roll angle and lateral acceleration while buffeting in the experiments. By substituting parameters for the “mean” bumblebee from the experiments into Eq. 1 the proportionality between the lateral acceleration and the roll angle due to the wind disturbance is around -77 m/s²/rad (see S2.2 for elaboration), which is comparable to that measured in the free flying bees (-74 ± 13 m/s²/rad), see Fig. 2d. In Eq. 1 the distance between the center of aerodynamic pressure and center of mass critically influences the relationship, thus the dynamics of the high frequency motion may differ vastly between insects.
**Bimodal Flight Control.** While foraging in the outdoor environment, flying organisms are likely to be subjected to wind induced disturbances over a wide range of frequencies and scales, however the magnitude of fluctuations over different frequencies and scales is still unclear due to lack of outdoor measurements. Nevertheless, performing active corrections for disturbances at all time scales is not only infeasible due to sensorimotor delays but also suboptimal. In this respect, it is vital to optimize the flight control system such that behaviorally relevant tasks can continue to be performed with minimal trade off. In the context of insect flight through unsteady winds, we distinguish between two distinct modes, i.e., the active and passive. Passive interactions with high frequency disturbances may be a suitable strategy since they impart relatively small net displacements that are likely to average over shorter time scales – as in this study where the high frequency displacements were smaller and tended to average over a few centiseconds, see Figs 1c and 5. Do the bees perform any corrective maneuvers when flying in the von Kármán street? The unfiltered time histories of the numerical simulations reveal that though the high frequency motions were relatively small in magnitude, the unsteady wind also introduces instabilities at lower frequencies, see SF2, and if uncorrected, these instabilities tend to increase over time. Therefore to maintain stable flight, some level of flight control is necessary but at frequencies much lower than the frequency of the von Kármán street. The von Kármán street induced disturbances at 23 Hz or 43 ms, which is within the perceptive limits and realms of sensorimotor response times of insects⁴. Fruitflies subjected to mechanically induced rapid roll perturbations responded within 5 ms with active changes in wing kinematics¹⁸. The passive response to higher frequencies disturbances noted in our experiments could be attributed to the relatively smaller magnitude of perturbations presented as well as the vast difference in organism properties. Additionally, insect responses to aerial perturbations have been shown to be context-specific whereby escape responses are mediated by specialized neural controllers that can be triggered within much smaller time scales¹¹.

Combination of active and passive modes in unsteady fluid environments is not uncommon in the animal kingdom. Classical experiments on fish swimming in von Kármán streets reveal a characteristic locomotion pattern known as Kármán gaiting where fish passively comply with the vortex street resulting in a combined bending and swaying of the body¹⁹. This passive mechanism may likely reduce sensorimotor demands, however additional active stabilization is still necessary but at much longer timescales, see ref. 20 for review. The flight trajectories of the bees here do not resemble Kármán gaiting, but the mechanism may be considered analogous. The sensorimotor basis for the implementation of such mechanisms in insects is still an open question while we have a better understanding of aquatic locomotion.

Measurements made on Kármán gaiting fish reveal relatively limited additional cost of locomotion in unsteady fluid flows²⁰ but little is known about the energetic cost of flight in steady vs unsteady winds²¹. However unlike in aquatic locomotion, the variations in orientation angles while maneuvering and high frequency motions inflight would result in a reduction in altitude due to the reorientation of the aerodynamic force vector. Thus apart from implementing corrections to maintain body posture in unsteady wind, an overall increase in force production may be necessary to maintain altitude. Additionally performing active corrections may impose further metabolic cost but measurements on bumblebees in rarified medium and in load-lifting tests reveal that they can increase force production with relatively limited energetic expense²²,²³. Taken together we observe a novel flight control mechanism implemented by bumblebees to cope with complex and unpredictable winds, a mechanism that could be widespread among volant organisms and one that could be suitable for artificial miniature aerial vehicles flying in outdoor environments.

## Materials and Methods

**Study specimens.** Bumblebees (*Bombus impatiens*) from a commercial breeder (BioBest) were maintained in the lab and given continuous access to a foraging chamber, where they could feed freely from an artificial flower containing linalool-scented nectar. Nominally similar sized bumblebees (body length = 14 mm ± 0.5 mm, mass = 165 mg ± 10%) were selected for flight experiments. The bees were placed in a transparent chamber (0.4 × 0.4 × 0.4 m) without access to food, for approximately two hours prior to the experiment.

**Flight Tests.** All experiments were conducted in a 6 m-long, suction-type, open-return wind tunnel with a 0.9 × 0.5 × 0.5 m working section. Once sufficiently starved, each bee was placed in the wind tunnel (with no airflow) where it could feed from an artificial flower resembling the one in the foraging chamber. Once feeding commenced, the bee was allowed to feed for approximately 10 seconds, then separated from the nectar source and released at the downstream end of the wind tunnel. Upon release, if the bee did not fly towards the artificial flower, it was manually re-introduced to the nectar source and subsequently separated. This procedure was repeated until direct flights to the nectar source were observed. Once consistent behavior was established, wind was introduced and bees were filmed as they flew upstream. The wind-speed was set to ~2.55 m/s, which represents an intermediate cruising velocity for bumblebees⁹,²⁴.

Each bee was flown in two airflow conditions, steady and unsteady. In the wind tunnel with unimpeded (steady) flow, a uniform velocity profile was present across the interrogation volume (<2% variation in mean flow speed) and turbulence intensity (standard deviation/mean wind speed) was less than 1.2%. With a cylinder (d = 25 mm) positioned vertically at the inlet of the test section, a von Kármán street developed in the wake and vortex shedding occurred at 23 Hz, in agreement with the predicted vortex shedding Strouhal number of ≈0.2⁵,²⁶. The vertical cylinder induced lateral velocity fluctuations at the shedding rate over a range of ±1.2 m/s. Refer to⁶ for further elaboration of flow conditions.

We filmed bees and quantified airflow within a specific interrogation volume (a cube with side lengths of 100 mm, located 100 mm downstream from the cylinder). The downstream distance was chosen to avoid the recirculating region in the near wake of the cylinder and to allow the formation of a periodic von Kármán vortex street. A total of 14 bees were subjected to this assay, and paired trials in steady vs. unsteady flow were obtained for each individual.
**Bimodal Flight Control.** While foraging in the outdoor environment, flying organisms are likely to be subjected to wind induced disturbances over a wide range of frequencies and scales, however the magnitude of fluctuations over different frequencies and scales is still unclear due to lack of outdoor measurements. Nevertheless, performing active corrections for disturbances at all time scales is not only infeasible due to sensorimotor delays but also suboptimal. In this respect, it is vital to optimize the flight control system such that behaviorally relevant tasks can continue to be performed with minimal trade off. In the context of insect flight through unsteady winds, we distinguish between two distinct modes, i.e., the active and passive. Passive interactions with high frequency disturbances may be a suitable strategy since they impart relatively small net displacements that are likely to average over shorter time scales – as in this study where the high frequency displacements were smaller and tended to average over a few centiseconds, see Figs 1c and 5. Do the bees perform any corrective maneuvers when flying in the von Kármán street? The unfiltered time histories of the numerical simulations reveal that though the high frequency motions were relatively small in magnitude, the unsteady wind also introduces instabilities at lower frequencies, see SF2, and if uncorrected, these instabilities tend to increase over time. Therefore to maintain stable flight, some level of flight control is necessary but at frequencies much lower than the frequency of the von Kármán street. The von Kármán street induced disturbances at 23 Hz or 43 ms, which is within the perceptive limits and realms of sensorimotor response times of insects⁴. Fruitflies subjected to mechanically induced rapid roll perturbations responded within 5 ms with active changes in wing kinematics¹⁸. The passive response to higher frequencies disturbances noted in our experiments could be attributed to the relatively smaller magnitude of perturbations presented as well as the vast difference in organism properties. Additionally, insect responses to aerial perturbations have been shown to be context-specific whereby escape responses are mediated by specialized neural controllers that can be triggered within much smaller time scales¹¹.

Combination of active and passive modes in unsteady fluid environments is not uncommon in the animal kingdom. Classical experiments on fish swimming in von Kármán streets reveal a characteristic locomotion pattern known as Kármán gaiting where fish passively comply with the vortex street resulting in a combined bending and swaying of the body¹⁹. This passive mechanism may likely reduce sensorimotor demands, however additional active stabilization is still necessary but at much longer timescales, see ref. 20 for review. The flight trajectories of the bees here do not resemble Kármán gaiting, but the mechanism may be considered analogous. The sensorimotor basis for the implementation of such mechanisms in insects is still an open question while we have a better understanding of aquatic locomotion.

Measurements made on Kármán gaiting fish reveal relatively limited additional cost of locomotion in unsteady fluid flows²⁰ but little is known about the energetic cost of flight in steady vs unsteady winds²¹. However unlike in aquatic locomotion, the variations in orientation angles while maneuvering and high frequency motions inflight would result in a reduction in altitude due to the reorientation of the aerodynamic force vector. Thus apart from implementing corrections to maintain body posture in unsteady wind, an overall increase in force production may be necessary to maintain altitude. Additionally performing active corrections may impose further metabolic cost but measurements on bumblebees in rarified medium and in load-lifting tests reveal that they can increase force production with relatively limited energetic expense²²,²³. Taken together we observe a novel flight control mechanism implemented by bumblebees to cope with complex and unpredictable winds, a mechanism that could be widespread among volant organisms and one that could be suitable for artificial miniature aerial vehicles flying in outdoor environments.

## Materials and Methods

**Study specimens.** Bumblebees (*Bombus impatiens*) from a commercial breeder (BioBest) were maintained in the lab and given continuous access to a foraging chamber, where they could feed freely from an artificial flower containing linalool-scented nectar. Nominally similar sized bumblebees (body length = 14 mm ± 0.5 mm, mass = 165 mg ± 10%) were selected for flight experiments. The bees were placed in a transparent chamber (0.4 × 0.4 × 0.4 m) without access to food, for approximately two hours prior to the experiment.

**Flight Tests.** All experiments were conducted in a 6 m-long, suction-type, open-return wind tunnel with a 0.9 × 0.5 × 0.5 m working section. Once sufficiently starved, each bee was placed in the wind tunnel (with no airflow) where it could feed from an artificial flower resembling the one in the foraging chamber. Once feeding commenced, the bee was allowed to feed for approximately 10 seconds, then separated from the nectar source and released at the downstream end of the wind tunnel. Upon release, if the bee did not fly towards the artificial flower, it was manually re-introduced to the nectar source and subsequently separated. This procedure was repeated until direct flights to the nectar source were observed. Once consistent behavior was established, wind was introduced and bees were filmed as they flew upstream. The wind-speed was set to ~2.55 m/s, which represents an intermediate cruising velocity for bumblebees⁹,²⁴.

Each bee was flown in two airflow conditions, steady and unsteady. In the wind tunnel with unimpeded (steady) flow, a uniform velocity profile was present across the interrogation volume (<2% variation in mean flow speed) and turbulence intensity (standard deviation/mean wind speed) was less than 1.2%. With a cylinder (d = 25 mm) positioned vertically at the inlet of the test section, a von Kármán street developed in the wake and vortex shedding occurred at 23 Hz, in agreement with the predicted vortex shedding Strouhal number of ≈0.2⁵,²⁶. The vertical cylinder induced lateral velocity fluctuations at the shedding rate over a range of ±1.2 m/s. Refer to⁶ for further elaboration of flow conditions.

We filmed bees and quantified airflow within a specific interrogation volume (a cube with side lengths of 100 mm, located 100 mm downstream from the cylinder). The downstream distance was chosen to avoid the recirculating region in the near wake of the cylinder and to allow the formation of a periodic von Kármán vortex street. A total of 14 bees were subjected to this assay, and paired trials in steady vs. unsteady flow were obtained for each individual.
One trial was conducted for each bee in steady and unsteady winds and the order of airflow conditions presented was pseudo-randomized among the individuals. The experimental data analyzed here are the same as that used in ref. 6.

**Kinematic reconstruction and analysis.** Prior to experiments, a triangular marker was affixed to the dorsal surface of the thorax of each bee, to aid in quantifying its position and orientation during flight. The markers consisted of three black points representing the vertices of an isosceles triangle (measuring $2.7 \times 2.3$ mm) set upon a white background. During flight trials, bees were filmed as they flew through the interrogation volume using two Photron SA3 high-speed cameras sampling at 1000 Hz, placed above the wind tunnel at approximately 30° from the vertical. The recorded flight sequences were digitized using an open-source MATLAB-based routine, DLTdv5²⁷, utilizing the automated tracking feature to localize the three black points on the triangular marker throughout each sequence. The RMS of digitization error in localizing the centroids of marker points was less than 3 pixels in all trials, which is much smaller than the mean number of pixels separating the markers (~30). The digitized position data were initially passed through a 4^th^-order, Butterworth, low-pass filter to remove any higher frequency errors due to the digitization process, with a cutoff frequency of 30 Hz, which is lower than the Nyquist frequency (500 Hz). For analysis of displacements at higher frequencies, the displacements were filtered using a high pass filter with 10 Hz cut-off.

The Center of Mass (COM) of bumblebees was estimated in ref. 28 by taking mass and morphology measurements of individual sections of the bees and subsequently estimating COM based on the bee's posture inflight. The mass and body size of the bees used for experiments here were similar to that used in ref. 28 hence those measurements were used for COM estimation. This was considered reasonable and additionally the bees appeared to maintain posture in all flight trial conducted here. The COM during flight for bumblebees is generally located in the mid-transverse plane and immediately anterior to the petiole, connection between thorax and abdomen²⁸. In our experiments, for each flight trial among all individuals, in the first video frame of all cameras the COM was manually located based on the morphology, subsequently for that frame the relationship between the rigid triangular marker on the thorax and the COM was obtained. This was done by creating a local coordinate system and evaluating the rotation matrix relating the COM and the points on the rigid marker. For all other frames, upon digitization of the points on the triangular marker, the COM rotation matrix was used to calculate COM location at each instant in the flight trajectory.

Instantaneous velocity and accelerations of the bees were calculated through numeric differentiation of the COM, see Supplementary Figure 1 for coordinate system. The influence of flow conditions on the body orientation and rotation rates of bees was assessed by evaluating variation in roll, pitch and yaw angles of the triangular markers, using a rigid body assumption. The method detailed in ref. 6 was used to calculate the instantaneous orientation and rotation rates of the bees, and to evaluate the errors associated with data capture and subsequent analysis. Subsequently, during analysis the lateral and rotational velocities and accelerations were separated into low (<10 Hz) and high frequency components of motion by applying a 4^th^ order Butterworth 10-Hz low- or high-pass filter.

For each individual's flight trial we evaluated the instantaneous (zero-time-shifted) cross correlation between the roll angle and the lateral acceleration. This was performed for all individuals for flights in steady winds and unsteady winds at different frequency ranges. Subsequently the cross correlation coefficient for all individuals was pooled depending on airflow conditions and frequency ranges. To elucidate the relationship between the bee's roll and lateral translation we performed regression analysis whereby for each flight we calculated the slope of the linear line of best fit relating the instantaneous roll angle with the lateral acceleration for that flight.

**Numerical Modeling.** The numerical method employed in this study is described in greater detail in refs 29 and 30, and the bumblebee model has been introduced in ref. 31, based on flight data obtained in earlier experiments³². In this section we provide a short overview of the model, we refer to Supplementary Figure 1 for further elaboration. The bumblebee is approximated by three rigid elements: the body and two wings, which move with respect to each other. The characteristic size of the insect is the wing length, which is equal to R = 13.2 mm. The wings follow a prescribed periodic flapping motion with frequency f = 152 Hz, which is close to what bumblebees use in forward flight at 2.5 m/s³². Their positional angle varies as $\Phi = 24^\circ + 115^\circ/2 \sin(2\pi ft)$, the elevation with respect to the stroke plane is constant and equal to $\theta = 12.55^\circ$. The feathering angle $\alpha$ is equal to 70° during upstroke and -40° during downstroke, with sinusoidal variation at the reversals. The angle between the body longitudinal axis and the horizontal plane is equal to 24.5°, and the angle between the wing stroke plane and the horizontal plane is equal to 28°.

These wing kinematics ensure trimmed flight at $u_\infty = 2.5$ m/s, assuming that the body mass is 175 mg. The position of the insect varies in time according to the solid-body dynamics model²⁹. In this study, only two degrees of freedom are taken into account: lateral displacement and roll rotation about the longitudinal axis of the body. The roll moment of inertia is equal to $I_{xx} = 10.92 \times 10^{-10}$ kg m². As with the experimental data gathered from bees in unsteady flow, the calculated translational and rotational positions and accelerations of the simulated bee were initially filtered using 4^th^ order Butterworth low-pass filter with cut-off frequency 30 Hz. Subsequently to isolate the higher frequency motions, the data was high-pass filtered using a 4^th^ order Butterworth 10 Hz high-pass filter.

To generate unsteady flow in the simulation, a vertical cylinder (diameter = 25 mm, Re = 4000) is placed in front of the insect, at approximately 80 mm from the insect's center of mass (Fig. 3, Supplementary Fig. 1) from the beginning of the simulation. A salient detail on one side of the cylinder breaks the symmetry of the flow, that starts impulsively and forms a vortex street after about 50 ms. The flow is governed by the incompressible Navier-Stokes equations, and the surrounding air has density $\rho = 1.177$ kg/m³ and kinematic viscosity $\nu = 1.57 \times 10^{-5}$ m²/s. The flow domain is a rectangular channel of length 132 mm, having a $105.6 \times 52.8$ mm cross-section. The no-slip boundary condition at the surfaces of the insect and cylinder as well as the flow outlet condition, are imposed using the volume penalization method (see refs 29 and 30 for details). The penalized
One trial was conducted for each bee in steady and unsteady winds and the order of airflow conditions presented was pseudo-randomized among the individuals. The experimental data analyzed here are the same as that used in ref. 6.

**Kinematic reconstruction and analysis.** Prior to experiments, a triangular marker was affixed to the dorsal surface of the thorax of each bee, to aid in quantifying its position and orientation during flight. The markers consisted of three black points representing the vertices of an isosceles triangle (measuring $2.7 \times 2.3$ mm) set upon a white background. During flight trials, bees were filmed as they flew through the interrogation volume using two Photron SA3 high-speed cameras sampling at 1000 Hz, placed above the wind tunnel at approximately 30° from the vertical. The recorded flight sequences were digitized using an open-source MATLAB-based routine, DLTdv5²⁷, utilizing the automated tracking feature to localize the three black points on the triangular marker throughout each sequence. The RMS of digitization error in localizing the centroids of marker points was less than 3 pixels in all trials, which is much smaller than the mean number of pixels separating the markers (~30). The digitized position data were initially passed through a 4th-order, Butterworth, low-pass filter to remove any higher frequency errors due to the digitization process, with a cutoff frequency of 30 Hz, which is lower than the Nyquist frequency (500 Hz). For analysis of displacements at higher frequencies, the displacements were filtered using a high pass filter with 10 Hz cut-off.

The Center of Mass (COM) of bumblebees was estimated in ref. 28 by taking mass and morphology measurements of individual sections of the bees and subsequently estimating COM based on the bee's posture inflight. The mass and body size of the bees used for experiments here were similar to that used in ref. 28 hence those measurements were used for COM estimation. This was considered reasonable and additionally the bees appeared to maintain posture in all flight trial conducted here. The COM during flight for bumblebees is generally located in the mid-transverse plane and immediately anterior to the petiole, connection between thorax and abdomen²⁸. In our experiments, for each flight trial among all individuals, in the first video frame of all cameras the COM was manually located based on the morphology, subsequently for that frame the relationship between the rigid triangular marker on the thorax and the COM was obtained. This was done by creating a local coordinate system and evaluating the rotation matrix relating the COM and the points on the rigid marker. For all other frames, upon digitization of the points on the triangular marker, the COM rotation matrix was used to calculate COM location at each instant in the flight trajectory.

Instantaneous velocity and accelerations of the bees were calculated through numeric differentiation of the COM, see Supplementary Figure 1 for coordinate system. The influence of flow conditions on the body orientation and rotation rates of bees was assessed by evaluating variation in roll, pitch and yaw angles of the triangular markers, using a rigid body assumption. The method detailed in ref. 6 was used to calculate the instantaneous orientation and rotation rates of the bees, and to evaluate the errors associated with data capture and subsequent analysis. Subsequently, during analysis the lateral and rotational velocities and accelerations were separated into low (<10 Hz) and high frequency components of motion by applying a 4th order Butterworth 10-Hz low- or high-pass filter.

For each individual's flight trial we evaluated the instantaneous (zero-time-shifted) cross correlation between the roll angle and the lateral acceleration. This was performed for all individuals for flights in steady winds and unsteady winds at different frequency ranges. Subsequently the cross correlation coefficient for all individuals was pooled depending on airflow conditions and frequency ranges. To elucidate the relationship between the bee's roll and lateral translation we performed regression analysis whereby for each flight we calculated the slope of the linear line of best fit relating the instantaneous roll angle with the lateral acceleration for that flight.

**Numerical Modeling.** The numerical method employed in this study is described in greater detail in refs 29 and 30, and the bumblebee model has been introduced in ref. 31, based on flight data obtained in earlier experiments³². In this section we provide a short overview of the model, we refer to Supplementary Figure 1 for further elaboration. The bumblebee is approximated by three rigid elements: the body and two wings, which move with respect to each other. The characteristic size of the insect is the wing length, which is equal to $R = 13.2$ mm. The wings follow a prescribed periodic flapping motion with frequency $f = 152$ Hz, which is close to what bumblebees use in forward flight at $2.5 \text{ m/s}³²$. Their positional angle varies as $\Phi = 24° + 115°/2 \sin(2\pi ft)$, the elevation with respect to the stroke plane is constant and equal to $\theta = 12.55°$. The feathering angle $\alpha$ is equal to $70°$ during upstroke and $-40°$ during downstroke, with sinusoidal variation at the reversals. The angle between the body longitudinal axis and the horizontal plane is equal to $24.5°$, and the angle between the wing stroke plane and the horizontal plane is equal to $28°$.

These wing kinematics ensure trimmed flight at $u_\infty = 2.5 \text{ m/s}$, assuming that the body mass is 175 mg. The position of the insect varies in time according to the solid-body dynamics model²⁹. In this study, only two degrees of freedom are taken into account: lateral displacement and roll rotation about the longitudinal axis of the body. The roll moment of inertia is equal to $I_{xx} = 10.92 \times 10^{-10} \text{ kg m}^2$. As with the experimental data gathered from bees in unsteady flow, the calculated translational and rotational positions and accelerations of the simulated bee were initially filtered using 4th order Butterworth low-pass filter with cut-off frequency 30 Hz. Subsequently to isolate the higher frequency motions, the data was high-pass filtered using a 4th order Butterworth 10 Hz high-pass filter.

To generate unsteady flow in the simulation, a vertical cylinder (diameter = 25 mm, $Re = 4000$) is placed in front of the insect, at approximately 80 mm from the insect's center of mass (Fig. 3, Supplementary Fig. 1) from the beginning of the simulation. A salient detail on one side of the cylinder breaks the symmetry of the flow, that starts impulsively and forms a vortex street after about 50 ms. The flow is governed by the incompressible Navier-Stokes equations, and the surrounding air has density $\rho = 1.177 \text{ kg/m}^3$ and kinematic viscosity $\nu = 1.57 \times 10^{-5} \text{ m}^2/\text{s}$. The flow domain is a rectangular channel of length 132 mm, having a $105.6 \times 52.8 \text{ mm}$ cross-section. The no-slip boundary condition at the surfaces of the insect and cylinder as well as the flow outlet condition, are imposed using the volume penalization method (see refs 29 and 30 for details). The penalized
equations are solved numerically using a Fourier pseudo-spectral method. The boundary condition on the side-walls of the channel is periodic. The flow domain is discretized using a uniform Cartesian grid consisting of 1920 × 1536 × 768 points (>>2.2 billion grid points). The volume penalization parameter is equal to η = 2.5 × 10⁻⁴. The resulting simulated flow environment accurately reproduced the unsteady flow environment used in experiments, with a vortex shedding rate and transient velocities matching those measured in the experiments.

**Statistical Analysis.** Statistical significance of experimental results was analyzed by performing paired t-tests (n = 14 individuals in all cases) in MATLAB, between the low-frequency components of flight in steady vs. unsteady flow, and between the low- and high-frequency components of flight in unsteady flow. The mean lateral and rotational displacement and accelerations obtained from the simulation were compared to those measured on freely flying bees in the experiments by performing a one-sample t-test.

**References**

1. Shyy, W., Kang, C.-K., Chirarattananon, P., Ravi, S. & Liu, H. Aerodynamics, sensing and control of insect-scale flapping-wing flight. *Proc. R. Soc. A Math. Phys. Eng. Sci.* **472** (2016).

2. Stull, R. B. *An Introduction to Boundary Layer Meteorology*. (Springer Netherlands, 1988). doi: 10.1007/978-94-009-3027-8.

3. Watkins, S., Milbank, J. & Loxton, B. J. & Melbourne, W. H. Atmospheric Winds and Their Implications for Microair Vehicles. *AIAA J.* **44**, 2591–2600 (2006).

4. Vance, J. T., Faruque, I. & Humbert, J. S. Kinematic strategies for mitigating gust perturbations in insects. *Bioinspir. Biomim.* **8**, 016004 (2013).

5. Ortega-Jimenez, V. M., Greeter, J. S. M., Mittal, R. & Hedrick, T. L. Hawkmoth flight stability in turbulent vortex streets. *J. Exp. Biol.* **216**, 4567–4579 (2013).

6. Ravi, S., Crall, J. D., Fisher, A. & Combes, S. A. Rolling with the flow: bumblebees flying in unsteady wakes. *J. Exp. Biol.* **216**, 4299–4309 (2013).

7. Ortega-Jimenez, V. M., Mittal, R. & Hedrick, T. L. Hawkmoth flight performance in tornado-like whirlwind vortices. *Bioinspir. Biomim.* **9**, 025003 (2014).

8. Crall, J. & Combes, S. Blown in the wind: Bumblebee temporal foraging patterns in naturally varying wind conditions. In *Integrative and Comparative Biology* E270 (Oxford Univ. Press Inc. 2013).

9. Riley, J. R. et al. Compensation for wind drift by bumble-bees. *Nature* **400**, 126–126 (1999).

10. Dudley, R. Mechanisms and implications of animal flight maneuverability. *Integr. Comp. Biol.* **42**, 135–140 (2002).

11. Muijres, F. T., Elzinga, M. J., Melis, J. M. & Dickinson, M. H. Flies evade looming targets by executing rapid visually directed banked turns. *Science* **344**, 172–177 (2014).

12. Shyy, W., Aono, H., Kang, C.-K. & Liu, H. *An Introduction to Flapping Wing Aerodynamics* (Cambridge University Press, 2013).

13. Sun, M. Insect flight dynamics: Stability and control. *Rev. Mod. Phys.* **86**, 615–646 (2014).

14. Boeddeker, N., Dittmar, L., Stürzl, W. & Egelhaaf, M. The fine structure of honeybee head and body yaw movements in a homing task. *Proc. Biol. Sci.* **277**, 1899–1906 (2010).

15. Braun, E., Dittmar, L., Boeddeker, N. & Egelhaaf, M. Prototypical components of honeybee homing flight behavior depend on the visual appearance of objects surrounding the goal. *Front. Behav. Neurosci.* **6**, 1 (2012).

16. Thomas, A. L. R. & Taylor, G. K. Animal Flight Dynamics II. Longitudinal Stability in Flapping Flight. *J. Theor. Biol.* **214**, 351–370 (2001).

17. Ros, I. G., Bassman, L. C., Badger, M. A., Pierson, A. N. & Biewener, A. A. Pigeons steer like helicopters and generate down- and upstroke lift during low speed turns. *Proc. Natl. Acad. Sci. USA* **108**, 19990–19995 (2011).

18. Beatus, T., Guckenheimer, J. M. & Cohen, I. Controlling roll perturbations in fruit flies. *J. R. Soc. Interface* **12**, 20150075–20150075 (2015).

19. Liao, J. C., Beal, D. N., Lauder, G. V. & Triantafyllou, M. S. Fish exploiting vortices decrease muscle activity. *Science* **302**, 1566–1569 (2003).

20. Liao, J. C. A review of fish swimming mechanics and behaviour in altered flows. *Philos. Trans. R. Soc. Lond. B. Biol. Sci.* **362**, 1973–1993 (2007).

21. Ortega-Jimenez, V. M., Sapir, N., Wolf, M., Variano, E. A. & Dudley, R. Into turbulent air: size-dependent effects of von Kármán vortex streets on hummingbird flight kinematics and energetics. *Proc. Biol. Sci.* **281**, 20140180 (2014).

22. Dillon, M. E. & Dudley, R. Surpassing Mt. Everest: extreme flight performance of alpine bumble-bees. *Biol. Lett.* **10**, 20130922 (2014).

23. Hedenstrom, A., Ellington, C. P. & Wolf, T. J. Wing wear, aerodynamics and flight energetics in bumblebees (Bombus terrestris): an experimental study. *Funct. Ecol.* **15**, 417–422 (2001).

24. Ellington, C. P. Limitations on animal flight performance. *J. Exp. Biol.* **160**, 71–91 (1991).

25. Roshko, A. Experiments on the flow past a circular cylinder at very high Reynolds number. *J. Fluid Mech.* **10**, 345 (1961).

26. Vickery, B. J. Fluctuating lift and drag on a long cylinder of square cross-section in a smooth and in a turbulent stream. *J. Fluid Mech.* **25**, 481 (1966).

27. Hedrick, T. L. Software techniques for two- and three-dimensional kinematic measurements of biological and biomimetic systems. *Bioinspir. Biomim.* **3**, 034001 (2008).

28. Mountcastle, A. M., Ravi, S. & Combes, S. A. Nectar vs. pollen loading affects the tradeoff between flight stability and maneuverability in bumblebees. *Proc. Natl. Acad. Sci.* **112**, 10527–10532 (2015).

29. Engels, T., Kolomenskiy, D., Schneider, K. & Sesterhenn, J. FluSI: A novel parallel simulation tool for flapping insect flight using a Fourier method with volume penalization. *SIAM J. Sci. Comput.* (2016).

30. Kolomenskiy, D. & Schneider, K. A Fourier spectral method for the Navier–Stokes equations with volume penalization for moving solid obstacles. *J. Comput. Phys.* **228**, 5687–5709 (2009).

31. Engels, T., Kolomenskiy, D., Schneider, K., Lehmann, F.-O. & Sesterhenn, J. Bumblebee flight in heavy turbulence. *Phys. Rev. Lett.* **116**, 028103 (2016).

32. Dudley, R. & Ellington, C. P. Mechanics of forward flight in bumblebees: I. kinematics and morphology. *J. Exp. Biol.* **148**, 19 (1990).

**Acknowledgements**

This work was granted access to the HPC resources of IDRIS (Institut du Développement et des Ressources en Informatique Scientifique) under project number i20152a1664. SR and DK gratefully acknowledges the financial support from the JSPS (Japan Society for the Promotion of Science) Postdoctoral Fellowship, JSPS KAKENHI Grant Number 15F15718 and 15F15061 respectively. TE, KS, and JS thank the French-German University for financial support. The authors would also like to thank Sacey Combes, James Crall and Nick Gravish for their valuable comments.
equations are solved numerically using a Fourier pseudo-spectral method. The boundary condition on the side-walls of the channel is periodic. The flow domain is discretized using a uniform Cartesian grid consisting of 1920 × 1536 × 768 points (>>2.2 billion grid points). The volume penalization parameter is equal to η = 2.5 × 10⁻⁴. The resulting simulated flow environment accurately reproduced the unsteady flow environment used in experiments, with a vortex shedding rate and transient velocities matching those measured in the experiments.

**Statistical Analysis.** Statistical significance of experimental results was analyzed by performing paired t-tests (n = 14 individuals in all cases) in MATLAB, between the low-frequency components of flight in steady vs. unsteady flow, and between the low- and high-frequency components of flight in unsteady flow. The mean lateral and rotational displacement and accelerations obtained from the simulation were compared to those measured on freely flying bees in the experiments by performing a one-sample t-test.

**References**

1. Shyy, W., Kang, C.-K., Chirarattananon, P., Ravi, S. & Liu, H. Aerodynamics, sensing and control of insect-scale flapping-wing flight. *Proc. R. Soc. A Math. Phys. Eng. Sci.* **472** (2016).

2. Stull, R. B. *An Introduction to Boundary Layer Meteorology*. (Springer Netherlands, 1988). doi: 10.1007/978-94-009-3027-8.

3. Watkins, S., Milbank, J. & Loxton, B. J. & Melbourne, W. H. Atmospheric Winds and Their Implications for Microair Vehicles. *AIAA J.* **44**, 2591–2600 (2006).

4. Vance, J. T., Faruque, I. & Humbert, J. S. Kinematic strategies for mitigating gust perturbations in insects. *Bioinspir. Biomim.* **8**, 016004 (2013).

5. Ortega-Jimenez, V. M., Greeter, J. S. M., Mittal, R. & Hedrick, T. L. Hawkmoth flight stability in turbulent vortex streets. *J. Exp. Biol.* **216**, 4567–4579 (2013).

6. Ravi, S., Crall, J. D., Fisher, A. & Combes, S. A. Rolling with the flow: bumblebees flying in unsteady wakes. *J. Exp. Biol.* **216**, 4299–4309 (2013).

7. Ortega-Jimenez, V. M., Mittal, R. & Hedrick, T. L. Hawkmoth flight performance in tornado-like whirlwind vortices. *Bioinspir. Biomim.* **9**, 025003 (2014).

8. Crall, J. & Combes, S. Blown in the wind: Bumblebee temporal foraging patterns in naturally varying wind conditions. In *Integrative and Comparative Biology* E270 (Oxford Univ. Press Inc. 2013).

9. Riley, J. R. et al. Compensation for wind drift by bumble-bees. *Nature* **400**, 126–126 (1999).

10. Dudley, R. Mechanisms and implications of animal flight maneuverability. *Integr. Comp. Biol.* **42**, 135–140 (2002).

11. Muijres, F. T., Elzinga, M. J., Melis, J. M. & Dickinson, M. H. Flies evade looming targets by executing rapid visually directed banked turns. *Science* **344**, 172–177 (2014).

12. Shyy, W., Aono, H., Kang, C.-K. & Liu, H. *An Introduction to Flapping Wing Aerodynamics* (Cambridge University Press, 2013).

13. Sun, M. Insect flight dynamics: Stability and control. *Rev. Mod. Phys.* **86**, 615–646 (2014).

14. Boeddeker, N., Dittmar, L., Stürzl, W. & Egelhaaf, M. The fine structure of honeybee head and body yaw movements in a homing task. *Proc. Biol. Sci.* **277**, 1899–1906 (2010).

15. Braun, E., Dittmar, L., Boeddeker, N. & Egelhaaf, M. Prototypical components of honeybee homing flight behavior depend on the visual appearance of objects surrounding the goal. *Front. Behav. Neurosci.* **6**, 1 (2012).

16. Thomas, A. L. R. & Taylor, G. K. Animal Flight Dynamics II. Longitudinal Stability in Flapping Flight. *J. Theor. Biol.* **214**, 351–370 (2001).

17. Ros, I. G., Bassman, L. C., Badger, M. A., Pierson, A. N. & Biewener, A. A. Pigeons steer like helicopters and generate down- and upstroke lift during low speed turns. *Proc. Natl. Acad. Sci. USA* **108**, 19990–19995 (2011).

18. Beatus, T., Guckenheimer, J. M. & Cohen, I. Controlling roll perturbations in fruit flies. *J. R. Soc. Interface* **12**, 20150075–20150075 (2015).

19. Liao, J. C., Beal, D. N., Lauder, G. V. & Triantafyllou, M. S. Fish exploiting vortices decrease muscle activity. *Science* **302**, 1566–1569 (2003).

20. Liao, J. C. A review of fish swimming mechanics and behaviour in altered flows. *Philos. Trans. R. Soc. Lond. B. Biol. Sci.* **362**, 1973–1993 (2007).

21. Ortega-Jimenez, V. M., Sapir, N., Wolf, M., Variano, E. A. & Dudley, R. Into turbulent air: size-dependent effects of von Kármán vortex streets on hummingbird flight kinematics and energetics. *Proc. Biol. Sci.* **281**, 20140180 (2014).

22. Dillon, M. E. & Dudley, R. Surpassing Mt. Everest: extreme flight performance of alpine bumble-bees. *Biol. Lett.* **10**, 20130922 (2014).

23. Hedenstrom, A., Ellington, C. P. & Wolf, T. J. Wing wear, aerodynamics and flight energetics in bumblebees (Bombus terrestris): an experimental study. *Funct. Ecol.* **15**, 417–422 (2001).

24. Ellington, C. P. Limitations on animal flight performance. *J. Exp. Biol.* **160**, 71–91 (1991).

25. Roshko, A. Experiments on the flow past a circular cylinder at very high Reynolds number. *J. Fluid Mech.* **10**, 345 (1961).

26. Vickery, B. J. Fluctuating lift and drag on a long cylinder of square cross-section in a smooth and in a turbulent stream. *J. Fluid Mech.* **25**, 481 (1966).

27. Hedrick, T. L. Software techniques for two- and three-dimensional kinematic measurements of biological and biomimetic systems. *Bioinspir. Biomim.* **3**, 034001 (2008).

28. Mountcastle, A. M., Ravi, S. & Combes, S. A. Nectar vs. pollen loading affects the tradeoff between flight stability and maneuverability in bumblebees. *Proc. Natl. Acad. Sci.* **112**, 10527–10532 (2015).

29. Engels, T., Kolomenskiy, D., Schneider, K. & Sesterhenn, J. FluSI: A novel parallel simulation tool for flapping insect flight using a Fourier method with volume penalization. *SIAM J. Sci. Comput.* (2016).

30. Kolomenskiy, D. & Schneider, K. A Fourier spectral method for the Navier–Stokes equations with volume penalization for moving solid obstacles. *J. Comput. Phys.* **228**, 5687–5709 (2009).

31. Engels, T., Kolomenskiy, D., Schneider, K., Lehmann, F.-O. & Sesterhenn, J. Bumblebee flight in heavy turbulence. *Phys. Rev. Lett.* **116**, 028103 (2016).

32. Dudley, R. & Ellington, C. P. Mechanics of forward flight in bumblebees: I. kinematics and morphology. *J. Exp. Biol.* **148**, 19 (1990).

**Acknowledgements**

This work was granted access to the HPC resources of IDRIS (Institut du Développement et des Ressources en Informatique Scientifique) under project number i20152a1664. SR and DK gratefully acknowledges the financial support from the JSPS (Japan Society for the Promotion of Science) Postdoctoral Fellowship, JSPS KAKENHI Grant Number 15F15718 and 15F15061 respectively. TE, KS, and JS thank the French-German University for financial support. The authors would also like to thank Sacey Combes, James Crall and Nick Gravish for their valuable comments.
**Author Contributions**

S.R. and D.K. conceptualized the experiments and simulations. S.R. performed the experiments while, D.K. and T.E. performed the simulations. S.R., D.K. and T.E. analyzed the results. All authors contributed to writing the manuscript.

**Additional Information**

**Supplementary information** accompanies this paper at http://www.nature.com/srep

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Ravi, S. *et al*. Bumblebees minimize control challenges by combining active and passive modes in unsteady winds. *Sci. Rep.* **6**, 35043; doi: 10.1038/srep35043 (2016).

<figure><img src="image_8.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2016
## Author Contributions

S.R. and D.K. conceptualized the experiments and simulations. S.R. performed the experiments while, D.K. and T.E. performed the simulations. S.R., D.K. and T.E. analyzed the results. All authors contributed to writing the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at http://www.nature.com/srep

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Ravi, S. *et al*. Bumblebees minimize control challenges by combining active and passive modes in unsteady winds. *Sci. Rep.* **6**, 35043; doi: 10.1038/srep35043 (2016).

<figure><img src="image_11.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2016