Abstract
Birds morph their wing shape to accomplish extraordinary manoeuvres1,2,3,4, which are governed by avian-specific equations of motion. Solving these equations requires information about a birdâs aerodynamic and inertial characteristics5. Avian flight research to date has focused on resolving aerodynamic features, whereas inertial properties including centre of gravity and moment of inertia are seldom addressed. Here we use an analytical method to determine the inertial characteristics of 22 species across the full range of elbow and wrist flexion and extension. We find that wing morphing allows birds to substantially change their roll and yaw inertia but has a minimal effect on the position of the centre of gravity. With the addition of inertial characteristics, we derived a novel metric of pitch agility and estimated the static pitch stability, revealing that the agility and static margin ranges are reduced as body mass increases. These results provide quantitative evidence that evolution selects for both stable and unstable flight, in contrast to the prevailing narrative that birds are evolving away from stability6. This comprehensive analysis of avian inertial characteristics provides the key features required to establish a theoretical model of avian manoeuvrability.
Similar content being viewed by others
Main
There is currently no theory that provides hypotheses to guide studies of avian manoeuvrability. This is not owing to a lack of physical understanding; manoeuvrability can be broadly defined as a birdâs ability to change the magnitude and direction of its velocity vector7,8. Similar to comparable uncrewed aerial vehicles, a birdâs flight dynamics and thus its manoeuvrability are dictated by its governing equations of motion. For example, aircraft dynamics depend on a minimum of six equations, three translational and three rotational, that can be derived from Newtonâs second law and its rotational counterpart5,9:
Where v is velocity vector and Ï is the angular velocity vector. These equations can be combined to solve for a flyerâs acceleration (translationally: \(\frac{{\rm{d}}{\bf{v}}}{{\rm{d}}t}\) and rotationally: \(\frac{{\rm{d}}{\boldsymbol{\omega }}}{{\rm{d}}t}\)), but this requires knowledge of both the aerodynamically informed external forces (F) and moments (M) as well as the inertial characteristics, including the mass (m) and moment of inertia tensor(I). However, avian inertial characteristics are not currently available with sufficient breadth or resolution.
Therefore, avian flight manoeuvrability is often evaluated experimentally by tracking individuals to measure accelerations during observed manoeuvres1,3,4. However, tracking data do not provide a birdâs maximal manoeuvring capabilities or allow extrapolation to unobserved behaviours. Determining these attributes requires a robust and general framework for manoeuvrability, equivalent to the manoeuvrability equations for aircraft8,10. Obtaining generalizable data is further complicated because aerodynamic and inertial characteristics vary substantially within and among species, and even dynamically for an individual bird1,11. For example, birds can initiate manoeuvres by morphingâthat is, changing the orientation and shape of their wings, body and tail7,12,13. To progress towards a theoretical formulation of avian manoeuvrability, there has been a marked and justifiable focus on resolving the aerodynamic characteristics of a bird in flight14,15,16. However, studies often overlook the equally essential inertial properties (Fig. 1a) or use static morphology approximations for individual species13,17,18,19,20. Here we fill this gap by investigating the variable inertial characteristics of flying birds to provide the necessary next step towards establishing a general framework of avian manoeuvrability.
Another challenge to solving a flying birdâs equations of motion is how to properly formulate the equations. For example, the equations can be simplified by defining the origin at the centre of gravity (Fig. 1a), which is equivalent to the centre of mass in a constant gravitational field9. If the centre of gravity moves substantially relative to the body, additional terms in the equations are required to properly capture flight dynamics10. Physically shifting a birdâs morphology shifts the centre of gravity, but it is not known how much the centre of gravity moves as a bird morphs. In addition, the rotational inertiaâquantified by the mass moment of inertia tensor (I) about the originâis also affected by morphing (Fig. 1a, b). This symmetric matrix describes the body mass distribution, where diagonal elements quantify the distribution relative to the major axes (Ixx, roll; Iyy, pitch; and Izz, yaw) and off-diagonal elements quantify distribution within the three major geometric planes9 (only Ixz is non-zero for symmetric configurations; Fig. 1a).
We calculated a birdâs centre of gravity and I to evaluate avian manoeuvrability through the lens of agility and static stability. Agility encompasses a birdâs ability to perform linear accelerations (axial agility) and angular accelerations (torsional agility)7, and depends on both the centre of gravity5 and I. In contrast, static stability refers to the initial tendency to return towards an equilibrium after a disturbance14. We quantified static pitch stability with the static margin, which is the distance between the centre of gravity and neutral point5,16 (Fig. 1a). If the neutral point is behind the centre of gravity, the static margin will be positive and thus stable. Often, stability is inversely related to agility because larger manoeuvring forces and moments are sometimes necessary to overcome stabilizing forces and moments14.
To determine how inertial characteristics vary during wing morphing, we developed a general analytical method to quantify any flying birdâs centre of gravity and I, and used a comparative analysis to investigate 22 species spanning the phylogeny defined by Prum et al.21, except for Palaeognathae as this clade contains largely flightless birds. First, we measured geometric and mass properties of cadavers and used motion tracking on cadaveric wings to extract the range of extension and flexion for the elbow and wrist (Fig. 1g). We limited our study to solely investigate the role of wing morphing due to elbow and wrist flexion and extension because previous studies have shown that this range of motion (ROM) enables a substantial shift in the neutral point14,16. The investigated ROM defines a birdâs physical capability to adjust its inertial characteristics and includes wing configurations outside of those probably used in flight. In addition, we assumed that the shoulder was set to allow a comparable wing orientation (see Methods) and that the tail is furled, but these degrees of freedom have an important role in avian flight control22 and warrant future morphing studies. Finally, we developed an open-source R package (AvInertia) that models birds as a composite structure of simple geometric objects and uses morphological data to calculate the centre of gravity and I for any bird using any wing configuration (Fig. 1câf, Methods). We validated this methodology with previous static wing measurements (Fig. 1h, Methods).
Centre of gravity is relatively constant
With our validated results, we first tested the effect of the elbow and wrist ROM on a birdâs centre of gravity when its wings are held symmetrically. We found that the ROM had a minimal effect on the position of the centre of gravity (Fig. 2b, opaque polygons). The maximum shifts along the x-axis and z-axis (\({x}_{\bar{{\rm{C}}}{\rm{G}}}\) and \({z}_{\bar{{\rm{C}}}{\rm{G}}}\); normalized by the full birdâs lengthâthe subscript CG refers to the centre of gravity) were 3% (great blue heron (Ardea herodias), 2.0 cm) and 2% (barn owl (Tyto alba), 0.7cm), respectively (Fig. 2b). Despite the small magnitude, wrist extension consistently shifted \({x}_{\bar{{\rm{C}}}{\rm{G}}}\) forwards (Pâ<â0.002) and the wrist angle explained a high amount of variance in the data leading to a high effect size, quantified by partial eta-squared (η2)23, 24. We found that partial η2âwas greater thanâ0.34 for all species (Fig. 2e). Similarly, elbow extension tended to shift \({x}_{\bar{{\rm{C}}}{\rm{G}}}\) forwards, but its effect size varied across species. Both elbow and wrist extension predominately shifted \({z}_{\bar{{\rm{C}}}{\rm{G}}}\) dorsally, but the magnitude and effect size varied. We could not differentiate the log-transformed mean \({x}_{\bar{{\rm{C}}}{\rm{G}}}\) or \({z}_{\bar{{\rm{C}}}{\rm{G}}}\) position from those expected if birds were simply scaled by preserving all length scales (that is, isometry) (Fig. 2f, Extended Data Table 1).
The small effect of the elbow and wrist on the location of the centre of gravity led us to question whether this would carry over to shoulder joint motion as well. To obtain a conservative estimate, we assumed that wings could rotate about the humeral head by 90° forwards, backwards, up and down (Fig. 2b, transparent squares). This revealed that the maximum \(\triangle {x}_{\bar{{\rm{C}}}{\rm{G}}}\) and \(\triangle {z}_{\bar{{\rm{C}}}{\rm{G}}}\) shifts were 18% (10.9âcm) for the great blue heron, approximately sixfold greater than that achieved with elbow and wrist morphing alone. Such a large shift in the centre of gravity probably cannot be neglected when formulating the equations of motion. At the other extreme, the Lady Amherstâs pheasant (Chrysolophus amherstiae) had a negligible shift of 1% (1.4âcm) with shoulder joint motion. Across the full range of taxa, we found a significant positive relationship between \({x}_{\bar{{\rm{C}}}{\rm{G}}}\) due to shoulder motion and the ratio of maximum wingspan to body length (Fig. 2d, Extended Data Table 1). This trend suggests that proper modelling of flight dynamics for birds with wings substantially longer than their body length will require an estimation of the expected centre of gravity shift to verify whether a fixed centre of gravity is an appropriate assumption.
Although the full birdâs centre of gravity defines its symmetric flight dynamics, the wing-only parameters can give insight into asymmetric configurations. We found that the elbow and wrist ROM caused the centre of gravity of the wing to shift along the y axis (\(\triangle {y}_{{\bar{{\rm{C}}}{\rm{G}}}_{{\rm{w}}{\rm{i}}{\rm{n}}{\rm{g}}}}\), normalized by the maximum half span) by between 10% (black swift (Cypseloides niger)) and 27% (American white pelican (Pelecanus erythrorhynchos)) (Fig. 2c), where the most distal \({y}_{{\bar{{\rm{C}}}{\rm{G}}}_{{\rm{w}}{\rm{i}}{\rm{n}}{\rm{g}}}}\) was 28% (western grebe (Aechmophorus occidentalis)). Additionally, \(\triangle {y}_{{\bar{{\rm{C}}}{\rm{G}}}_{{\rm{w}}{\rm{i}}{\rm{n}}{\rm{g}}}}\) was positively associated with the arm-to-hand wing ratio (Extended Data Table 1), such that birds with longer hand wings than arm wings (like the swift) would have a reduced capacity to shift the wingâs centre of gravity. The centre of gravity shift was largely driven by elbow extension (Pâ<â0.001, partial η2â>â0.51; Fig. 2e) whereas the effect of the wrist varied across species. These results highlight a well-conserved proximal location of the wing centre of gravity across species. Contrary to a previous study25, we did not find that the log-transformed mean \({y}_{{\bar{{\rm{C}}}{\rm{G}}}_{{\rm{w}}{\rm{i}}{\rm{n}}{\rm{g}}}}\) differed from isometric expectations (Fig. 2g, Extended Data Table 1).
Morphing affects lateral inertia
The centre of gravity is crucial to formulating the governing equations, but their solution depends on a birdâs rotational inertia. Like the centre of gravity, we found that a birdâs rotational inertia (log-transformed mean diagonal components of I) scaled isometrically with body mass (Fig. 3a, Extended Data Table 1). However, we found that elbow and wrist extension provided a more than 11-fold Ixx increase (heron) and a 3-fold Izz increase (heron and owl) (Fig. 3c). This capability was largely driven by elbow extension (Fig. 3b), which had a significant effect on both Ixx (Pâ<â0.001, partial η2â>â0.23; except for Leachâs storm petrel (Hydrobates leucorhous)) and Izz (Pâ<â0.009, partial η2â>â0.45). The absolute values of Iyy and Ixz were minimally affected by joint extension and the effect size varied substantially across species (Fig. 3b). We next computed the contribution of each major body part to the overall rotational inertia for birds with wings at maximum elbow and wrist extension (Fig. 3dâf). Because the wings were extended along the y axis, this captures approximately the lowest wing contribution to Iyy but the highest wing contribution to Ixx. The percentage contribution of each body part varied substantially across the species, but as expected the wings were responsible for the majority of Ixx. These results indicate that elbow and wrist ROM provides substantial inertial control over the roll and yaw axes (Ixx, Izz), but less so for the pitch axis (Iyy), although species-specific differences were also apparent in our results. Incorporating the shoulder joint ROM would increase the wingâs contribution to inertial pitch control.
Inertia informs the pitch agility metric
We next tested whether inertial characteristics could be used to estimate a birdâs pitch agility. However, because both inertia and aerodynamics are fundamental to flight dynamics, we first used aerodynamic theory and data from a rigid gull wing16 to obtain an estimate for the neutral point, and thus the static margin for each configuration (Methods and Supplementary Methods). Using these results, we derived a novel pitch agility metric that is proportional to the angular acceleration about the y axis due to a change in the angle of attack (a form of torsional agility; Fig. 1a, Methods). Note that agility in a stable configuration indicates that the developed acceleration would tend to return the bird towards an equilibrium position. We found that the pitch agility range decreases as body mass increases, which was expected because flight speed and body size scale positively with mass26 (Fig. 4a, Extended Data Table 1). These results are further driven by the static margin whose range also decreases as mass increases (Fig. 4b, Extended Data Table 1). Incorporating the shoulder joint ROM would broaden the static margin range because the resultant neutral point shift is probably larger than the centre of gravity shift as evidenced by morphing uncrewed aerial vehicles with shoulder-inspired joints27,28.
Evolutionary pressures on stability
Next, we looked for evidence of selective evolutionary pressures on avian pitch agility and stability. We investigated the static margin specifically because it is both a component of the pitch agility metric and dictates the static stability of a flying bird. We identified the configurations with the maximum and minimum static margin for each individual (Extended Data Fig. 1) and then calculated the mean of each trait for each species (Fig. 4b). We found that four species were entirely stable, one species was entirely unstable, and 17 species had the capacity to shift between stable and unstable flight (Fig. 4b, e). Using these data, we found that an Ornstein Uhlenbeck model was significantly favoured over a Brownian motion model for both the maximum (ÎAICc = â8.24; Extended Data Fig. 2b) and minimum static margin (ÎAICc = â5.01; Extended Data Fig. 2c), where AICc is the Akaike information criterion with correction for smaller sample sizes. Further, we found that the optimal static margin phenotype (θsm) was stable for the maximum static margin (26% of the maximum root chord, strength of selection (αOU) = 0.53, variance (Ï2) = 14.2 Ã 10â3), whereas the optimal phenotype for the minimum static margin was unstable (â15% of the maximum root chord, αOU = 0.06, Ï2 = 2.7 Ã 10â3) (Fig. 4b). This suggests that evolutionary pressures act to maintain birdsâ ability to transition between stable and unstable flight. The strength of selection (αOU) was relatively low, but our results were robust to measurement errors (Extended Data Figs. 3, 4) and to a preliminary estimation of a neutral point shift due to the tail (Supplementary Methods). Further, an Ornstein Uhlenbeck model was a good fit for the mean \({x}_{\bar{{\rm{C}}}{\rm{G}}}\)such that the phenotypic optimum (θCG) was 10% of the body length behind the humeral head (ÎAICcâ=ââ8.23, αOU = 0.11, Ï2 = 0.1 Ã 10â3; Extended Data Fig. 2a). The stability of this centre of gravity position depends on the location of the neutral point (Fig. 4c, d).
Although studies have suggested that modern birds may be capable of stabilized flight13,14,16, it is widely believed that birds have evolved to be unstable in pitch to enhance manoeuvrability6. Our results offer a new perspective on the evolution of avian flight: evolutionary pressures may be maintaining the ability to shift between stable and unstable configurations. Elbow and wrist flexion and extension alone offer birds the capacity to shift between these pitch stability modes but, if and when a flying bird does shift between these modes remains to be seen. As highlighted by Thomas and Taylor13, dynamically switching between stable and unstable modes probably requires substantially different control algorithms, and thus switching between these modes would necessitate a complex flight control system. Further, our findings offer insight on how birds perform slow glides with positive tail lift29. By maintaining the capacity to relocate the wingâbody neutral point in front of the centre of gravity, birds may achieve an equilibriumâalbeit unstableâflight condition.
It is important to highlight that further work is required to incorporate the inter- and intra-specific aerodynamic capabilities, shoulder and tail ROM, and in vivo configurations to definitively confirm the optimal phenotype(s) for static pitch stability. We expect that the shoulder joint will enhance the available pitch control and the ability to shift between modes owing to an increased static margin range; the extent of this enhancement will depend on each speciesâ shoulder ROM5,27,28. Future work is also required to extend this analysis to the roll and yaw axes to discuss lateral agility and stability, which will need to account for aerodynamic and inertial coupling5. Finally, 23% of the species in our study were unable to shift between stable and unstable modes with the elbow and wrist alone, and thus there are many combinations of stability characteristics in modern birds.
Conclusions
In summary, our results reveal that elbow and wrist ROM have a small relative effect on the centre of gravity location and pitch inertia, but have a substantial effect on the roll and yaw inertia. Although inter- and intra-specific variation is apparent, we found that the measured range of wrist and elbow motion alone is sufficient to enable switching between stable and unstable flight in 17 out of 22 bird species. Further, an evolutionary analysis shows that the phenotypic optimum maximum and minimum static margin supports the ability to transition between stable and unstable flight, suggesting the need for a complex flight control system. Collectively, investigating the inertial characteristics of flying birds throughout elbow and wrist ROM brings us a step closer to establishing a fundamental theory to quantify and evaluate avian manoeuvrability.
Methods
Collection of morphological data
We obtained morphological data for 36 adult specimens representing 22 species (Fig. 2a) from frozen cadavers acquired from the Cowan Tetrapod Collection at the Beaty Biodiversity Museum (University of British Columbia, Vancouver, Canada). Sample size was a function of the availability and quality of specimens from the museum as we could only rely on fully-intact, well-preserved specimens. The cadavers were inspected to ensure adequate condition and completeness, after which we measured the full body mass, wingspan, and body length. Next, we disarticulated the wing at the shoulder joint, taking care to ensure that each wingâs skin, propatagial elements, and feathers remained intact. One wing from each cadaver was used to determine wing ROM and corresponding wing shape change (see âDetermination of the elbow and wrist ROMâ). The cadaver was further dissected to obtain length and mass measurements for the head, neck, torso, wing components, legs, and tail (refer to Supplementary Methods for details on each measurement). We obtained the centre of gravity coordinates for the torso (body without head, neck, tail, wings) by manually balancing the torso and measuring the distance from the clavicle reference point to the balanced position. Note that because of the preservation of the storm petrel specimens, we estimated the mass on the basis of humerus bone length and the torso centre of gravity as being proportional to that of the gull. Finally, we individually weighed and photographed each flight feather, enabling geometric parameters to be extracted using ImageJ software30. Refer to the publicly available data for details on all assumptions used for extracting the morphological measurements. Note that this study consisted of a single experimental group and thus randomization and blinding was not necessary.
Determination of the elbow and wrist ROM
To determine the wing ROM and corresponding shape change, we actuated the cadaver wings throughout the full range of extension and flexion of the elbow and wrist joints by hand (following methods established by Baliga et al.11, Fig. 1g). We tracked the location of 10 reflective markers each 4Â mm in diameter (grey and white points in Fig. 1câf, refer to Supplementary Methods for details) with automated 3D data capture at 30 frames per second using a 4- or 5-camera tracking system (OptiTrack, NaturalPoint). Using tools from NaturalPoint, each recording was calibrated to have less than 0.5âmm overall mean reprojection error. Joint angles were calculated as the interior angle defined by three key points: points 1, 2 (vertex) and 3 for the elbow, and points 2, 3 (vertex) and 4 for the wrist (Supplementary Methods).
Developing AvInertia
We developed an open source R package (AvInertia) to calculate the centre of gravity and moment of inertia tensor (I) for any flying bird (Fig. 1a) in RStudio31 (version 1.3.1093) running R32 (version 4.0.3). A high-level overview of the code methodology follows in this section. Further details are provided in the Supplementary Methods, as each individual component of the avian models required specific procedures and approximations.
To allow a generalized approach, we used a common methodology from mechanics to estimate the centre of gravity and inertia components using simple geometric shapes9. We elected to use as many elements as possible to allow the best resolution. For each species, we first modelled the birdâs body without the wings as a composite of five components: head, neck, torso, legs and tail. To determine the inertial properties of the wings, we aligned each wing configuration extracted from the ROM measurements so that the wrist joint was in line with the shoulder joint along the y and z axes and so that the wrist joint was aligned with the first secondary feather (S1) along the x axis (extended wing: Fig. 1c, d; folded wing: Fig. 1e, f). Note that this positioning results in a different shoulder angle between each wing configuration and wings with extremely low elbow angles and high wrist angles being positioned at substantially different incidence angles than the body. Each wing was then modelled as a composite of twelve components: bones (humerus, radius, ulna, carpometacarpus/digit, radiale and ulnare), muscles (brachial, antebrachial and manus groups), skin, coverts, and tertiary feathers. In addition, each primary and secondary feather was modelled and positioned individually as a composite structure of five components: calamus, rachis (cortex exterior and medullar interior), and distal and proximal vanes. AvInertia permits a variable number of flight feathers. With our methodology, a bird with 10 primaries and 10 secondaries that flies with an extended neck will be represented by a composite model with 232 individual simple geometric shapes. In our study, we investigated only symmetric wing configurations for a full bird and considered the effects of a single wing independently. We assumed that anisotropic effects such as the air space within the body would have a minimal impact on the overall centre of gravity33.
To calculate the final inertial characteristics of this composite bird, each componentâs shape, mass, and positioning was informed by its corresponding morphological measurements. We began by determining the centre of gravity and I for one of the basic geometric shapes with respect to an origin and frame of reference that simplified the formulation of the centre of gravity and I for that shape. Next, AvInertia computed the mass-weighted summation of the centre of gravity of each object and shifted the origin to the bird reference point, located at the centre of the spinal cord when cut at the clavicle. The centre of gravity was then transformed into the full bird frame of reference, which is defined by Fig. 1câf. We used the parallel axis theorem and the appropriate transformation matrices to transform I to be defined about the final centre of gravity within the full bird frame of reference.
Validating AvInertia
We validated our methodology by comparing the maximum rotational inertia about the roll axis for a single wing (\({{I}_{{xx}}}_{{\rm{wing}}}\), origin at the humeral head) to data from previous experimental studies that measured \({{I}_{{xx}}}_{{\rm{wing}}}\) by cutting an extended wing into strips25,34 (Fig. 1h). Our 95% confidence intervals on the exponent of body mass marginally overlapped with Berg and Raynerâs predictions25 but were significantly lower than Kirkpatrickâs predictions34. However, Kirkpatrick used 10 wing strips while Berg and Rayner later found that at least 15 strips were necessary to minimize systematic error25,34. Next, we directly compared results for the pigeon (Columba livia), the only species in common between the studies, and found \({{I}_{{xx}}}_{{\rm{wing}}}\)(Ã104) was between 1.42 and 1.92âkg m2, which encompasses values from previous studies25,34,35 (1.72 and 1.83âkg m2). The pigeon wingâs maximum centre of gravity position along the y-axis (\({y}_{{\bar{{\rm{C}}}{\rm{G}}}_{{\rm{w}}{\rm{i}}{\rm{n}}{\rm{g}}}}\)) was only 3% of the half span more proximal than Berg and Raynerâs measurement25. We expect minor differences because strip methods enforce that all wing mass is contained within the xây plane while AvInertia accounts for out-of-plane morphology (Fig. 1h).
Agility and stability metrics
We developed a pitch agility metric that estimates the change of the angular acceleration about the y axis (\(\triangle \dot{q}\), known as the time rate of change of the pitch rate) due to a degree change in the angle of attack (Îα) as:
Where m is the body mass, crmax is the maximum root chord for the specimen, Smax is the maximum single wing area for the specimen, xCG is the centre of gravity position on the x-axis measured from the humeral head, and \({\widetilde{x}}_{c/4}\) is the quarter chord of the standard mean chord37 (defined in equation (7)). This equation was derived beginning from the rigid aircraft y axis rotational equation of motion assuming a symmetric configuration undergoing small disturbances5:
From this equation, we estimated the change in pitching moment (ÎÎ) with a Taylor series expansion method assuming that the largest effect is due to angle of attack and then non-dimensionalized as follows5:
Where Ï is air density, V is the freestream scalar velocity, and CM and CL are the coefficients of pitching moment and lift, respectively. Because the pitching moment slope \((\frac{{\rm{\partial }}{C}_{M}}{{\rm{\partial }}{C}_{L}})\) is proportional to static margin16,36, we estimated each configurationâs neutral point (indicated by the subscript NP) using our previous morphing gull wing-body aerodynamic results (see Supplementary Methods). This analysis revealed that the neutral point for a wing-body configuration scaled with:
\({\widetilde{x}}_{c/4}\) is the quarter chord of the standard mean chord defined as37:
Where b is the wingspan and, c and xc/4 are the chord and quarter chord location as a function of the span position (y), respectively. This equation was evaluated numerically for each of the bird wings modelled with 1,000 segments. In addition, we performed a sensitivity analysis on the exponent (see Supplementary Methods). With the estimated neutral point, we calculated the static margin as:
Refer to Supplementary Methods for further details pertaining to the aerodynamic assumptions.
For the pitch agility metric, we incorporated a previously established allometric scaling26 of cruise velocity (\(V\propto \) m0.12). We assumed a constant air density (Ï)Â and constant lift slope \((\frac{{\rm{\partial }}{C}_{L}}{{\rm{\partial }}\alpha })\) across species to obtain the final proportional relationship as:
This result was then returned to equation (4) and rearranged to obtain the pitch agility metric as seen in equation (3).
Phylogenetic and statistical analyses
All phylogenetically informed analyses were carried out using the time-calibrated maximum clade credibility tree from Baliga et al.11, which was pruned to the 22 focal taxa in this study. To determine the linear trends with body mass, we fit first-order phylogenetic generalized linear mixed models (PGLMM) to the data using the R package MCMCglmm38 where the random effects are informed by the phylogeny (Extended Data Table 1). Note that the linear trend of the pitch agility range with body mass remains significant even if the storm petrels are removed from the data. All PGLMM models had priors specified with the inverse Wishart scaling parameters Vâ=â1 and νâ=â0.02 and used 1.3 Ã 107 Markov chain Monte Carlo iterations. As visualized by Fig. 2f, we cannot definitively exclude the possibility that the lower 95% confidence interval on \({x}_{\bar{{\rm{C}}}{\rm{G}}}\) may be positive which would indicate that \({x}_{\bar{{\rm{C}}}{\rm{G}}}\) scales greater than isometric predictions. However, multiple MCMCglmm runs returned an insignificant result. To determine the significance and effect of the elbow and wrist on the centre of gravity and I components, we independently fit first order interactive models to each specimensâ data with a constant scaling on the elbow and wrist angle. We calculated the effect size of the elbow and wrist using the R package effectsize23 and independently fit first order interactive models to each specimensâ data with scaled and mean centred elbow and wrist angles.
Next, to investigate the phenotypic optimum of the pitch agility and stability traits\(,\) we independently fit both Brownian motion and Ornstein Uhlenbeck models to the absolute data using the R package geiger39. We assumed that all species belong to the same regime and thus, fit single-peak evolutionary models. This analysis revealed that there was evidence that the Ornstein Uhlenbeck model was a better model fit for all three of our selected traits (\({x}_{\bar{{\rm{C}}}{\rm{G}}}\) maximum and minimum static margin) due to a lower Akaike information criterion with correction for smaller sample sizes (AICc). Because of the smaller sample size of our study40, we ran a Monte Carlo simulation (nâ=â5,000) with the R package pmc41 to validate that selecting the Ornstein Uhlenbeck model over the Brownian motion model was appropriate (Extended Data Fig. 2). This method returns a distribution of likelihood ratios (twice the difference of the maximum log likelihood for each model) when the traits have been simulated n times under each model. These distributions are then compared to the observed likelihood ratio (black dashed vertical lines in Extended Data Fig. 2). For details, refer to Boettiger et al.41. We found that the likelihood ratio predicted by a Brownian motion model was more extreme than the observed ratio for the minority of simulations (\({x}_{\bar{{\rm{C}}}{\rm{G}}}\):0.2%, maximum static margin: 0.1%, minimum static margin: 1%). Further we had sufficient power to differentiate the two models as the majority of the simulations under the Ornstein Uhlenbeck model fell outside of 95th percentile of the Brownian motion distribution (\({x}_{\bar{{\rm{C}}}{\rm{G}}}\) :73.8%, maximum static margin: 77.2%, minimum static margin: 67.2%). 95% confidence intervals were constructed for each reported metric of each trait (Extended Data Table 2). Together these results provide confidence that the observed likelihood ratio of each trait is more likely to occur under an Ornstein Uhlenbeck model than a Brownian motion model.
Sensitivity analysis
Because both the pitch agility and stability metrics directly depend on \({x}_{{\rm{C}}{\rm{G}}}\), we investigated the sensitivity caused by shifting the combined torso and tail centre of gravity forwards and backwards by up to 15% of the torso. Note that for some species there was a physical limit to the ability to relocate the centre of gravity while maintaining the known morphological properties and if the shifted distance was larger than 4Â cm we removed it from the analysis as that was assumed to be an overestimate. The final estimated shift of the relative maximum and minimum static margin is shown in Extended Data Fig. 4. This sensitivity analysis revealed a minor effect on the parameters.
Finally, we wanted to investigate the potential effect of error in our measured centre of gravity metric on our key evolutionary results. To this end, we used a custom bootstrapping code (nâ=â5,000) and randomly sampled (with replacements) from each specimenâs centre of gravity error range used for the sensitivity analysis to recalculate the mean value of the minimum and maximum static margin for each species. With each of these new trait distributions, we re-fit an Ornstein Uhlenbeck model and extracted the optimal phenotype (Extended Data Fig. 3). We found that even allowing for this substantial centre of gravity error, all minimum static margin cases had an unstable optimum and all maximum static margin cases had a stable optimum (Extended Data Fig. 3). Note that this analysis is equivalent to both accounting for the same magnitude shift in the neutral point with a fixed centre of gravity as well as accounting for possible inter-specific variation within the error bounds shown in Extended Data Fig. 4.
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this paper.
Data availability
All data used in this study have been deposited in public repositories identified at https://doi.org/10.6084/m9.figshare.c.5503989.
Code availability
AvInertia has been publicly released on GitHub (https://github.com/charvey23/AvInertia) and as a CRAN package (https://cran.r-project.org/package=AvInertia). All analysis codes used in this study can be accessed from the public repositories identified at https://doi.org/10.6084/m9.figshare.c.5503989.
References
Dakin, R., Segre, P. S., Straw, A. D. & Altshuler, D. L. Morphology, muscle capacity, skill, and maneuvering ability in hummingbirds. Science 359, 653 (2018).
Warrick, D. R. The turning-and linear-maneuvering performance of birds: the cost of efficiency for coursing insectivores. Can. J. Zool. 76, 1063â1079 (1998).
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).
Gillies, J. A., Thomas, A. L. & Taylor, G. K. Soaring and manoeuvring flight of a steppe eagle Aquila nipalensis. J. Avian Biol. 42, 377â386 (2011).
Pamadi, B. N. Performance, Stability, Dynamics, and Control of Airplanes (American Institute of Aeronautics & Astronautics, 2004).
Smith, J. M. The importance of the nervous system in the evolution of animal flight. Evolution 6, 127â129 (1952).
Dudley, R. Mechanisms and implications of animal flight maneuverability. Integr. Comp. Biol. 42, 135â140 (2002).
Goman, M. G., Khramtsovsky, A. V. & Kolesnikov, E. N. Evaluation of aircraft performance and maneuverability by computation of attainable equilibrium sets. J. Guid. Control Dyn. 31, 329â339 (2008).
Baruh, H. Analytical Dynamics (WCB/McGraw-Hill Boston, 1999).
Obradovic, B. & Subbarao, K. Modeling of flight dynamics of morphing wing aircraft. J. Aircr. 48, 391â402 (2011).
Baliga, V. B., Szabo, I. & Altshuler, D. L. Range of motion in the avian wing is strongly associated with flight behavior and body mass. Sci. Adv. 5, eaaw6670 (2019).
Warrick, D. R., Bundle, M. W. & Dial, K. P. Bird maneuvering flight: blurred bodies, clear heads. Integr. Comp. Biol. 42, 141â148 (2002).
Thomas, A. L. R. & Taylor, G. K. Animal flight dynamics I. Stability in gliding flight. J. Theor. Biol. 212, 399â424 (2001).
Harvey, C., Baliga, V. B., Lavoie, P. & Altshuler, D. L. Wing morphing allows gulls to modulate static pitch stability during gliding. J. R. Soc. Interface 16, 20180641 (2019).
Lentink, D. et al. How swifts control their glide performance with morphing wings. Nature 446, 1082 (2007).
Harvey, C., Baliga, V. B., Goates, C. D., Hunsaker, D. F. & Inman, D. J. Gull-inspired joint-driven wing morphing allows adaptive longitudinal flight control. J. R. Soc. Interface 18, 20210132 (2021).
Hedrick, T. L., Usherwood, J. R. & Biewener, A. A. Wing inertia and whole-body acceleration: an analysis of instantaneous aerodynamic force production in cockatiels (Nymphicus hollandicus) flying across a range of speeds. J. Exp. Biol. 207, 1689â1702 (2004).
Ducci, G., Colognesi, V., Vitucci, G., Chatelain, P. & Ronsse, R. Stability and sensitivity analysis of bird flapping flight. J. Nonlinear Sci. 31, 47 (2021).
Mills, R., Hildenbrandt, H., Taylor, G. K. & Hemelrijk, C. K. Physics-based simulations of aerial attacks by peregrine falcons reveal that stooping at high speed maximizes catch success against agile prey. PLoS Comput. Biol. 14, e1006044 (2018).
Durston, N. E., Wan, X., Liu, J. G. & Windsor, S. P. Avian surface reconstruction in free flight with application to flight stability analysis of a barn owl and peregrine falcon. J. Exp. Biol. 222, jeb185488 (2019).
Prum, R. O. et al. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature 526, 569â573 (2015).
Cheney, J. A. et al. Raptor wing morphing with flight speed. J. R. Soc. Interface 18, 20210349 (2021).
Ben-Shachar, M. S., Lüdecke, D. & Makowski, D. effectsize: estimation of effect size indices and standardized parameters. J. Open Source Softw. 5, 2815 (2020).
Cohen, J. Eta-squared and partial eta-squared in fixed factor Anova designs. Educ. Psychol. Meas. 33, 107â112 (1973).
Berg, C. & Rayner, J. The moment of inertia of bird wings and the inertial power requirement for flapping flight. J. Exp. Biol. 198, 1655â1664 (1995).
Alerstam, T., Rosén, M., Bäckman, J., Ericson, P. G. P. & Hellgren, O. Flight speeds among bird species: allometric and phylogenetic effects. PLoS Biol. 5, e197 (2007).
Neal, D. et al. Design and wind-tunnel analysis of a fully adaptive aircraft configuration. in 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference 1727 (ARC, 2004).
Paranjape, A. A., Chung, S.-J. & Selig, M. S. Flight mechanics of a tailless articulated wing aircraft. Bioinsp. Biomim. 6, 026005 (2011).
Usherwood, J. R. et al. High aerodynamic lift from the tail reduces drag in gliding raptors. J. Exp. Biol. 223, jeb214809 (2020).
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671 (2012).
RStudio Team. RStudio: Integrated Development Environment for R (RStudio, PBC, 2020).
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2020).
Allen, V., Paxton, H. & Hutchinson, J. R. Variation in center of mass estimates for extant sauropsids and its importance for reconstructing inertial properties of extinct archosaurs. Anat. Rec. 292, 1442â1461 (2009).
Kirkpatrick, S. J. Short communication the moment of inertia of bird wings. J. Exp. Biol. 151, 489â494 (1990).
Pennycuick, C. J. & Lock, A. Elastic energy storage in primary feather shafts. J. Exp. Biol. 64, 677 (1976).
Anderson Jr, J. D. Fundamentals of Aerodynamics (Tata McGraw-Hill Education, 2010).
Yates, A. H. Notes on the mean aerodynamic chord and the mean aerodynamic centre of a wing. J. R. Aeronaut. Soc. 56, 461â474 (1952).
Hadfield, J. D. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J. Stat. Softw. 33, 1â22 (2010).
Pennell, M. W. et al. geiger v2.0: an expanded suite of methods for fitting macroevolutionary models to phylogenetic trees. Bioinformatics 30, 2216â2218 (2014).
Cooper, N., Thomas, G. H., Venditti, C., Meade, A. & Freckleton, R. P. A cautionary note on the use of Ornstein Uhlenbeck models in macroevolutionary studies. Biol. J. Linn. Soc. 118, 64â77 (2016).
Boettiger, C., Coop, G. & Ralph, P. Is your phylogeny informative? Measuring the power of comparative methods. Evolution 66, 2240â2251 (2012).
Acknowledgements
We thank I. Szabo for her help in acquiring specimens and E. R. Press for his assistance during data collection. This work is supported in part by the US Air Force Office of Scientific Research (AFOSR) under grant number FA9550-16-1-0087, titled âAvian-Inspired Multifunctional Morphing Vehiclesâ monitored by B. L. Lee, in part by the National Science Foundation (NSF) under grant number 1935216 and in part by the Natural Sciences and Engineering Research Council of Canada (NSERC) under grant number RGPIN-2016-05381. C.H. is further supported by a NSERC PGS-D, Zonta International Amelia Earhart Fellowship, and the François-Xavier Bagnoud Fellowship awarded by FXB International through the University of Michigan Department of Aerospace Engineering.
Author information
Authors and Affiliations
Contributions
C.H., J.C.M.W. and V.B.B. designed the study. J.C.M.W. and V.B.B. collected and processed all morphological measurements including the ROM videos. C.H. wrote AvInertia and analysed output data. D.L.A. and D.J.I. provided conceptual and technical guidance. C.H. wrote the manuscript. All authors edited the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks Aimy Wissa and the other, anonymous reviewer(s) for their contribution to the peer review of this work. Peer review reports are available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 The elbow and wrist angle configurations that yielded the maximum and minimum static margin for each species.
The range of motion investigated for each species with the maximum static margin and minimum static margin identified with a black diamond on each species. Note that the diamonds are coloured by the static margin.
Extended Data Fig. 2 A power analysis confirmed the validity of Ornstein Uhlenbeck models for three key traits.
We used a Monte Carlo-based method to investigate if our phylogeny provides support for the use of the Ornstein Uhlenbeck model. This returns the distribution of likelihood ratios under each model for (a) \({x}_{\bar{{\rm{C}}}{\rm{G}}}\), (b) maximum static margin, and (c) minimum static margin. Comparing the likelihood ratio distributions produced under both a Brownian motion (BM) model (grey) and Ornstein Uhlenbeck (OU) model (yellow) to the observed likelihood ratio (dashed black line) revealed that the Ornstein Uhlenbeck model was a better fit for our three key traits (see Methods).
Extended Data Fig. 3 Bootstrapping our results within the conservative centre of gravity measurement error supported our results.
We used a Monte Carlo method to investigate the impact of centre of gravity measurement error on the phenotypic optima for the maximum static margin (green) and minimum static margin (purple). Our results for a stable maximum and an unstable minimum are confirmed with over a 95% confidence interval.
Extended Data Fig. 4 Conservative measurement sensitivity analysis revealed a minimal effect on pitch stability and agility metrics.
We assumed that experimental error on the centre of gravity measured for the torso and tail was ±15% of the torso length (up to a maximum of 4âcm) and recalculated the (a) maximum pitch agility for stable flight (most negative values in Fig. 4a), (b) minimum and (c) maximum static margin. The estimated error ranges from panels b and c informed the bootstrapping analysis in Extended Data Fig. 3.
Supplementary information
Supplementary Methods
This file contains all Supplementary Methods pertaining to the aerodynamic analysis and implementation of AvInertia and associated Supplementary Figs. 1â7 and Supplementary Table 1.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Harvey, C., Baliga, V.B., Wong, J.C.M. et al. Birds can transition between stable and unstable states via wing morphing. Nature 603, 648â653 (2022). https://doi.org/10.1038/s41586-022-04477-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-022-04477-8
This article is cited by
-
Avian-inspired embodied perception in biohybrid flapping-wing robotics
Nature Communications (2024)
-
Fast ground-to-air transition with avian-inspired multifunctional legs
Nature (2024)
-
Biologically Inspired Pectoral Propulsors with Flapping and Rowing Control for a Specified Stroke Plane Angle
Journal of Nonlinear Science (2024)
-
Decoupling body shape and mass distribution in birds and their dinosaurian ancestors
Nature Communications (2023)
-
Morphological evolution of bird wings follows a mechanical sensitivity gradient determined by the aerodynamics of flapping flight
Nature Communications (2023)