Locomotor preferences in terrestrial vertebrates: An online crowdsourcing approach to data collection

Understanding how animals move within their environment is a burgeoning field of research. Despite this, relatively basic data, such as the locomotor speeds that animals choose to walk at in the wild, are sparse. If animals choose to walk with dynamic similarity, they will move at equal dimensionless speeds, represented by Froude number (Fr). Fr may be interpreted from simple limb kinematics obtained from video data. Here, using Internet videos, limb kinematics were measured in 112 bird and mammal species weighing between 0.61 and 5400 kg. This novel method of data collection enabled the determination of kinematics for animals walking at their self-selected speeds without the need for exhaustive fieldwork. At larger sizes, both birds and mammals prefer to walk at slower relative speeds and relative stride frequencies, as preferred Fr decreased in larger species, indicating that Fr may not be a good predictor of preferred locomotor speeds. This may result from the observation that the minimum cost of transport is approached at lower Fr in larger species. Birds walk with higher duty factors, lower stride frequencies and longer stance times compared to mammals at self-selected speeds. The trend towards lower preferred Fr is also apparent in extinct vertebrate species.

terrestrial locomotion, it is unclear if the preferred walking speeds and associated kinematics of wild animals are coincident with similar Fr.
Kinematic parameters are intimately linked to the work requirements of walking and therefore to the CoT (assuming a relationship between mechanical and metabolic cost) but are rarely quantified outside of the laboratory. Such field data, in which animals are self-selecting their walking speeds, are critical in understanding the pattern of energy use of animals in their natural environments. Contra to the predicted kinematics when animals are walking at dynamically similar speeds (Fr), self-selected DF increases with size in mammals with a slope proportional to M b 0.02 (Fig. 1a). Equivalent data for birds are unavailable. This increase in DF appears to result from larger species taking relatively longer stances as opposed to shorter swing durations (Fig. 1c,d). f stride scales as approximately M b −0.16 in freely walking mammals 1,3,8 and as M b −0.4 in birds 2 (Fig. 1b), a slope lower than that predicted if animals were walking at the same Fr 7 . The latter exponent, however, may be influenced by low sample sizes of species with diverse morphologies 9 and a low number of steps analysed (only 14 of the 25 species investigated had more than 5 steps analysed) 2 . Together, the available kinematic data suggest that animals may not self-select speeds that are dynamically similar, however these data are sparse.
To date, a comprehensive study of the self-selected walking speeds of terrestrial species and the associated kinematics is lacking. Here, using user-generated content from an online video database, the walking kinematics of a range of birds and mammals was quantified to test the hypothesis that at self-selected walking speeds, animals move in dynamically similar ways. Additionally, because the idea of dynamic similarity is underpinned by geometric similarity, the hypothesis that dynamic similarity would only hold within class groups was also investigated.

Results
DF increased at a faster rate with M b in birds (proportional to M b 0.03 ) than in mammals (proportional to M b 0.02 ) and was higher in birds at any given M b (Fig. 2a, Table 1 and Supplementary Table S2). In contrast, the incremental change in self-selected f stride with M b was similar for birds and mammals (Fig. 2b)

Discussion
The present study is the first to quantify locomotor kinematics at self-selected speeds across such a diversity of species, particularly birds. We reject the hypothesis that animals prefer to walk at dynamically similar speeds as defined by Fr. Previous data indicate that large animals have a lower DF than small animals walking at the same Fr ( Fig. 1a) 4,6 . Here, preferred DF increases with M b in both mammals and birds (Fig. 2a), providing strong evidence that larger animals walking relatively more slowly than small species. A preference for lower Fr with a coincident increase in DF has been demonstrated in 9 species of cat and the present data suggest that this pattern may be observed broadly across terrestrial vertebrates.
, an exponent lower than the lowest value for birds when moving at an equal Fr of 0.1 (Fig. 1b, ) again suggesting that larger species walk at relatively slower speeds. This exponent is not as low as that found previously in birds 2 or mammals 1,3,8 . These studies, however, contain much smaller and less diverse datasets compared to the data presented here. The relatively higher DF and relatively lower f stride in larger species indicating a lower Fr appears to be due to a combination of both relatively longer t stance and relatively shorter t swing in large animals. This is contra to the previous data in mammals that suggested a role of t stance alone and may result from differences in the passive swing properties of the legs with increasing size as well as differences in the active propulsion of the limbs. Regardless of the underlying mechanisms, the preference for relatively slower speeds at increased size appears to be a central feature of locomotion in terrestrial vertebrates despite large variations in morphology and life histories. This preference is even observed in extinct bipedal and quadrupedal dinosaurs when their preferred Froude is estimated from fossilized trackways (Fig. 3, Supplementary Table S3). In general, the rate of oxygen consumption (  Vo 2 , ml O 2 s −1 kg −1 ) increases linearly with walking speed and the rate of this increase (the CoT min ) is negatively correlated with M b 10 .  Vo 2 may be used to determine the CoT for species of different size (Fig. 4), by conversion to metabolic power (assuming an energetic equivalent of 20.1 J per ml of O 2 10 ). When plotted against speed, the asymptotes of the calculated curves occur at higher absolute speeds with increasing M b , although this pattern is weak (Fig. 4a). There is, however, a decrease in the Fr at which the CoT min for walking is approached in animals at larger sizes (Fig. 4b). This downward shift in optimal Fr can be observed in humans walking with artificially elongated limbs 11 and may result from the relatively longer effective limb lengths of large species for a given weight 12,13 . This downward shifting energetically optimal Fr may, in part, account for the lower relative speed preferences of larger birds and mammals. It remains unclear, however, why animals don't walk at their maximum walking speed at which DF approaches 0.5 and CoT is at its theoretical minimum. Either the CoT is not minimal at the fastest walking speeds, as is the case in larger species of bird 14 and mammal 15 which show U-shaped CoT curves, or there are non-energetic components of speed preference. For example, animals may have a size dependent optimal intrinsic frequency above and below which locomotor reflex responses are compromised 16 . In small species, behavioral adaptations towards reducing exposure to predation and the time travelling between foraging sites (i.e. moving relatively quickly) may be additional non-energetic drivers of higher preferred   Froude numbers estimated from dinosaur trackways over a large size range of quadrupeds and bipeds mirror our findings that smaller vertebrate species walk relatively more quickly (i.e. with a higher Froude number) than larger species. The references for these data can be found in Supplementary Table S3. walking speeds. Alternatively, larger species may walk relatively more slowly in order to minimise the chances of suboptimal foot placement, decreasing the chance of falling, which may be more costly to larger species. Although many locomotor studies focus on stance phase mechanics, inertial forces on the limb during the swing phase of walking may also influence preferred locomotor speeds 17 . Given the cost associated with swinging the limbs 18,19 , minimizing this source of work may contribute to self-selected locomotor speeds and kinematics in walking species. In order to mitigate the cost of swinging the limbs, animals may move with swing durations matching or proportional to the natural pendular periods (NPP) of their limbs 3 . This passive swing is predicted by compass gait mechanics and is seen, for example in waddling birds 20 . The preferred mammalian swing durations reported here are significantly lower than the half NPP of freely swinging mammalian limbs of similar size (Fig. 2d) 21 . Animals are therefore actively swinging their legs, as has been previously found in African ungulates whose limbs swing faster than if modelled as a pendulum 1 . Faster swing phases than predicted for a compound pendulum are also observed in humans and may be associated with a ballistic model in which muscles are active early in swing to propel the leg forward in an otherwise passive manner 22 . Our findings suggest that such a ballistic model of walking may be broadly applicable. Although swinging the limb at a rate faster than its NPP may be energetically costly, this cost may be mitigated by the geometric advantages of a smoother centre of mass path achieved with relatively shorter steps 23 . Although t swing values don't match the expected NPP values for limbs, the f stride data for both mammals and birds are in close agreement with the values of Pennycuick 3 and are proportional to the natural frequency of a simple pendulum with a length of one third of shoulder height (Fig. 2b). Care must be taken, however, in interpreting these data as in reality, limb morphology is complex and diverse. More data regarding the inertial and mass properties of the limbs of birds and mammals would allow refinement of the model, improving its predictive power.
With the exception of DF, the kinematics of both birds and mammals change with the same exponent of M b when walking at self-selected speeds. This suggests commonalities in the underlying mechanics driving preferred gaits regardless of the number of legs. Throughout the range of M b , however, birds have a higher t stance than mammals at their self-selected walking speed. t swing is similar between the two classes resulting in higher DF and lower f stride in birds. Higher DFs in birds are consistent with their more crouched posture, resulting from their compliant limb morphology. The energetic consequences of these kinematic differences between quadrupeds and bipeds at their preferred speeds are unknown. From DF alone the present data suggest that birds would have a lower CoT at their self-selected speed compared to mammals of the same M b as a result of increased foot contact time. Birds, however, may use 1.7 times more energy than quadrupeds for a given rate of force generation, making energetic predictions problematic 24 .
In summary, the hindlimb kinematics of animals walking at self-selected speeds suggest that larger animals (both extant and extinct) prefer to walk relatively more slowly and not at equal values of Fr. Thus, despite being a central tenet of terrestrial locomotor research, Fr and dynamic similarity cannot explain the movement of birds or mammals in the wild. As suggested previously, at their self-selected walking speed, animals walk with t swing faster than expected if the limbs were swinging as a simple pendulum. No animals choose to walk at their theoretical minimum DF (0.5) at which the cost of transport reaches its minimum value, indicating energy expenditure may not be the only driver of preferred walking speed. Finally, collecting video data across a broad range of animals is time consuming. Extracting data from web based video footage has limitations (for example, controlling camera orientation relative to the subject and a lack of a scale), yet, nonetheless, this 'crowdsourcing' approach to data collection offers a useful tool for the study of animal locomotion outside of the lab, which will no doubt improve in parallel with the camera technology available to the public.

Methods
Data collection. Videos of 112 walking bird (n = 32) and mammal species (n = 80) spanning 14 orders (M b = 0.61 to 5400 kg) from wild and captive (zoo) individuals were found on YouTube (the entire list of species and video sources can be found in Supplementary Table S1). The kinematic parameters of DF, f stride , t stance and t swing were calculated from the relative durations of the stance and swing phases of a single limb using Tracker 4.85 © video analysis software (Open Source Physics, Java framework). None of these measures require, either the image scale or the exact angle of the camera lens relative to the animal to be known. The criteria for inclusion in the analyses were that animals were walking whilst not feeding or ground foraging. Strides were analysed for between 1 and 5 individuals per species and for between 2 and 20 strides per individual depending on the availability and length of usable video data. M b for all individuals were estimated using published values. In some cases it was possible to determine an individual's sex. Where sex was indeterminate, the species' median M b was used.

Statistical analyses. Data were normally distributed (by investigation with quantile-quantile plots and
Shapiro-Wilk tests). Analysis of covariance (ANCOVA, two-sided, with a significance level of p < 0.05) was performed on log-transformed data to identify differences in the observed walking kinematics between class groups across all body masses. If the interaction term (class group x M b ) was non-significant indicating similar slopes, they were removed from the ANCOVA model and the test was rerun assuming parallel slopes (testing for differences in the intercepts only. All statistics were conducted in Matlab R2013a (The MathWorks, Inc., Natick, MA, USA). The lines of best fit were calculated from the coefficients table generated by the ANCOVA. The ANCOVA outputs are reported in Tables 2 and 3 of the Supplementary Information.