|
|
|
|||
| Home Help Feedback Subscriptions Archive Search Table of Contents | ||||
First published online December 14, 2006
Journal of Experimental Biology 210, 82-90 (2007)
Published by The Company of Biologists 2007
doi: 10.1242/jeb.02612
Analyzing the effect of wind on flight: pitfalls and solutions
1 Computational Biogeography and Physical Geography, Institute of
Biodiversity and Ecosystem Dynamics, University of Amsterdam, Nieuwe
Achtergracht 166, 1018 WV Amsterdam, The Netherlands
2 Swiss Ornithological Institute, 6204 Sempach, Switzerland
* Author for correspondence (e-mail: shamoun{at}science.uva.nl)
Accepted 19 October 2006
| Summary |
|---|
|
|
|---|
Key words: compensation, flight, model, spurious correlation, wind
| Introduction |
|---|
|
|
|---|
Numerous field and theoretical studies have tried to measure or predict how
birds and insects react to wind (e.g.
Liechti, 2006
;
Riley et al., 2003
;
Srygley and Oliveira, 2001
).
One hypothesis regarding migratory flight is that birds should maximize the
distance traveled for a given amount of fuel. In order to fulfill this
hypothesis birds are predicted to increase their air speeds in headwinds and
decrease their air speeds in tailwinds
(Liechti, 1995
;
Pennycuick, 1978
). Pennycuick
(Pennycuick, 1978
) further
proposed that this prediction could be tested by comparing the relationship,
initially assumed to be non-linear, between air speed (Va)
and the difference between ground speed (Vg) and
Va, both scalar quantities. It is noteworthy that this
difference is generally represented in the literature as
Vg-Va and termed the `speed increment
due to wind' or the `wind effect', where positive values of
Vg-Va represent tailwinds and negative
values represent headwinds. Pennycuick states "A `tail wind' is
conventionally defined as the scalar difference between ground speed and true
air speed. The `wind effect' means that a bird whose ground speed is less than
its air speed will normally respond by increasing its air speed, resulting in
a negative correlation between the air speed and `tail wind'"
[(Pennycuick, 2001
) p. 3288].
Perhaps as a result of the simplicity of this particular approach, the linear
relationship between Va and
Vg-Va has been tested in the
literature numerous times for birds
(Alerstam et al., 1993
;
Alerstam and Gudmundsson, 1999
;
Green and Alerstam, 2000
;
Gudmundsson et al., 2002
;
Hedenström and Alerstam,
1996
; Hedenström et al.,
1999
; Hedenström et al.,
2002
; Liechti et al.,
1994
; Pennycuick,
1982
; Pennycuick,
2001
; Rosen and
Hedenström, 2001
;
Wakeling and Hodgson, 1992
)
and for migratory insects (Srygley,
2003
). The results have been used to determine how birds alter
their air speed in relation to tailwinds and headwinds and to predict air
speeds in varying wind conditions. The overwhelming evidence from these
studies has often been used in support of the prediction that birds increase
their air speed in a headwind and decrease their air speed in a tailwind.
However, conflicting results were found when head- and tailwind situations
were separated (Hedenström et al.,
2002
). In several cases, no relationship was found and treated as
potential type II errors (Rosen and
Hedenström, 2001
), or no compensation for wind was made
(Alerstam et al., 1993
;
Rosen and Hedenström,
2001
). Although the initial prediction was for the specific case
of pure headwind or tailwinds, the prediction was expanded to include the
influence of side wind on optimal air speeds
(Liechti et al., 1994
);
however, no solution was provided for analyzing the effects of side and tail
winds simultaneously.
Both wind and flight are composed of two components that can be considered
either in the form of speed and direction or in the form of their x
and y vector components for a given coordinate system. If the
influence of wind on flight is studied in only one of these two dimensions
while either speed or direction vary, information is lost and erroneous
interpretations may result. To our knowledge, all studies investigating the
influence of wind on heading (e.g.
Srygley, 2003
;
Srygley et al., 1996
;
Wege and Raveling, 1984
) or
air speed (e.g. Able, 1977
;
Alerstam et al., 1993
;
Hedenström et al., 2002
;
Pennycuick, 2001
) have adopted
a one-dimensional model. With a one-dimensional model we mean a model with
only one explanatory variable. In this paper we will focus on the analysis of
air speed. We argue that the conventional analysis of air speed in relation to
wind, by testing the linear relationship between the scalar
Va and Vg-Va,
cannot be used to assess the relationship between air speed and wind nor, more
specifically, the prediction that birds should maximize the distance traveled
per fuel cost by increasing their air speeds in headwinds and decreasing it in
tailwinds. This paper provides an alternative approach to test how birds alter
their air speed in relation to wind speed and direction. Three different
datasets are used to illustrate the weakness of the conventional model as well
as the strengths of the newly proposed model for (1) simulated random data,
(2) simulated artificial data including an established influence of wind and
(3) measured autumn passerine migration and corresponding wind conditions.
| Materials and methods |
|---|
|
|
|---|
The relation between vector components, speed and direction
To study bird flight in relation to wind, we need three orthogonal vectors.
The first expresses a displacement per unit time of the bird with respect to
the ground (we will call this the ground vector, g) the second vector
expresses the displacement of the bird with respect to air (the air vector,
a) and the third expresses the displacement of the wind (the wind
vector, w). In this study we consider movement in the horizontal plane
and ignore vertical movement, hence our vectors have two elements only:
displacement in the x- and y-directions.
Fig. 1 gives a graphic
representation of this system, and the three vectors are defined as follows:
![]() | (1) |
|
![]() | (2a) |
![]() | (2b) |
![]() | (2c) |
Angles between the vectors and some reference direction can also be
conveniently calculated on the basis of the x and y
components, using the following equations:
![]() | (3a) |
![]() | (3b) |
![]() | (3c) |
Here
is known as the track direction,
the bird's heading
and
the wind direction. We have chosen to define north as the zero
angle. Following from the definition of positive x and y in
the eastern and northern directions respectively, the angles are positive in
the clockwise direction.
Obviously, the three vectors are not independent: the ground vector is the
sum of the air and wind vectors:
![]() | (4) |
In most studies, air speed (Va) and heading (
)
are not measured directly, nor are the x- and y-components.
Rather, ground speed (Vg) and direction (
) and wind
speed (Vw) and direction (
) are measured. Note that
when studying the influence of wind on flight, we are interested in the
direction wind is blowing to. Wind measurements received from climatological
surface stations often note meteorological wind direction, which is the
direction wind is blowing from.
On the basis of the available information, Va and
can be calculated. It is convenient to first calculate
xa and ya:
![]() | (5) |
The heading
can then be calculated by applying Eqn 3b.
Air speed is finally obtained by applying Eqn 2b using
xa and ya calculated in Eqn 5.
Although apparent from Eqn 4, it is important to note at this point that any
four out of the six variables determine the values of the other two. For
clarity we also show the full expression for Va as a
function of Vg,
, Vw,
in Eqn 6:
![]()
With Eqn 6 in mind we can review the relation between
Va and (Vg-Va),
whose linear relationship has been tested to detect a negative slope in avian
and entomological literature (see Introduction):
![]() | (7) |
By substituting Eqn 6 in Eqn 7, Eqn 7 becomes a rather complex implicit
relation. In this context, implicit means that Va occurs
at both the right and left hand side of the equation. The linearity
hypothesis, expressed by Eqn 7, is not valid in general, since
nVa
n(Vg-Va).
The functional relationship between Va and
(Vg-Va) depends on
Vg and Vw as well as on
and
and is generally not linear as often treated in the literature
(Fig. 2 and Fig. S1 in
supplementary material).
|
-
)=0 or (|
-
|)=180°, a pure
tail- or headwind, respectively, with respect to the ground vector (see
Fig. 2A for pure tailwind). In
this case,
Va=Vg-Vw; note
that although this is quite obvious, it also follows from Eqn 6 by setting
cos(
-
) to 1, and subsequently factoring the equation. This
equality implies that Vg can be eliminated in Eqn 7 so
that one obtains a function with only Vw as explanatory
variable (Va=a+bVw). Clearly,
many different functional relationships between Va and
Vg-Va may be expected, including
positive relationships, for example, with pure side winds in respect to the
ground vector and fairly constant wind speeds
(Fig. 2C). To provide intuitive
insight in the relation between Va and
Vg-Va, the two-dimensional relation
between g, a and w is depicted for a few points from
Fig. 2B
(Fig. 3).
|
-
=45° and
-
=135°
are shown. In this example, Vg and Vw
are independently generated random variables with a mean of 10 m
s-1. The variance of Vw equals 4 m2
s-2 and the variance of Vg equals 0, 1 and 4
m2 s-2 (Fig.
4, top to bottom). Va is calculated on the
basis of Vg,
, Vw and
.
As previously described, when Vg is constant, the relation
between Va and
Vg-Va is a straight line with slope
-1. However, if Vg varies, the amount of correlation
decreases, depending on
-
as well as the variability of
Vg relative to that of Vw.
|
As a more appropriate method of analysis we propose a two-dimensional model. Our method departs from the fact that Va is, by definition, a function of two variables, so any model for Va will have to be explicit (Va only on the left hand side of the equation) and two-dimensional. Secondly, we search for orthogonal variables to make maximum use of the information content in a dataset when fitting our model. In the third place, we try to specify variables that are meaningful in a physical or biological sense. Finally we try to use a modeling framework that is statistically well developed and offers large flexibility to model linear as well as non-linear systems.
The modeling framework that provides the flexibility that we are looking
for is a generalized additive model [GAM; see Guisan et al.
(Guisan et al., 2002
), and
references therein]. For this model we use the
2 statistic to
test whether a predictor explains a significant proportion of the variance and
other graphical diagnostics are used to evaluate the performance of a model
such as residual plots, quantile-quantile plots and Cook's distance to
identify potential outliers. In a GAM, deviance reduction is used as a measure
of model fit and the adjusted D2 value (comparable to the
adjusted R2), which takes into account the number of
predictors and observations, is used to determine the real fit of the model
and to compare models (e.g. Guisan and
Zimmermann, 2000
).
Two orthogonal variables that are related to Va and
have no implicit relationship with Va, or
Vg, are the components of wind along the x and
y axis. The two variables are implemented in a GAM by first
transforming them via a rather constrained LOESS smoother (a locally
weighted regression) with a maximum span of 80% and 1-2 degrees of freedom
(d.f.) (Cleveland, 1979
;
Cleveland and Devlin, 1988
).
The resulting GAM then has the following form:
![]() | (8) |
where the function f() refers to the transformation via a LOESS smoother. If the transformation is not justified, after initial model derivation, then variables are maintained in their linear form.
Testing both analyses
In this study we test the two analyses by performing an experiment using
data sets with known properties. A first data set is generated with no
relation between the air vector and wind vector at all (SRD, simulated random
data), therefore no wind compensation. A second data set is generated with a
very strong relation between the air and wind vectors (SWI, simulated wind
influence), to represent full wind compensation. The two analyses are applied
to both data sets. A correct analysis should obviously identify both
situations correctly.
Both methods are also applied to real data (a situation where the degree of compensation is unknown) to consider the different conclusions from both methods with respect to the real data.
Description of the data sets
Data set SRD, simulated random data, was generated to comprise 880
artificial data records of bird flight and wind. The wind components
(xw and yw) were generated
independently with a pseudo random number generator, using a Gaussian process
with a mean of 0 m s-1 and a standard deviation (s.d.) of 4 m
s-1. Similarly, xa and ya
were generated independently using a Gaussian process with a mean of 0 m
s-1 and s.d.=3 m s-1. Therefore, all four components are
entirely independent of each other and we should not find any relationship
between wind components and flight components. The ground vector is generated
as the sum of the air and wind vectors (viz. Eqn 4).
A second data set SWI, simulated wind influence, was created to represent a
strong influence of wind on Va and
. The dataset
contains 880 artificial data records of birds that adjust
Va and
to maintain constant Vg
and
while winds are variable. The ground components were chosen to be
xg=6 m s-1 and yg=3 m
s-1, with a small Gaussian noise added (using a mean of 0 m
s-1 and a s.d.=1 m s-1). The xw and
yw components were taken from the SRD data. The
xa and ya components were subsequently
calculated by xg-xw and
yg-yw and finally
Vg and Va were calculated (Eqn 2a,
2b).
The third dataset was taken from a radar field study in southern Germany
[for measurement details, see Liechti
(Liechti, 1993
)]. The data set
contains 880 radar tracks of autumn passerine migration (APM) with direct
observations of Vg,
, Vw and
. Va was computed from these using Eqn 6.
Rose plots of the headings and direction of flight for all three data sets are shown in Fig. 5.
|
| Results and discussion |
|---|
|
|
|---|
|
Analyzing air speed in two dimensions, as a function of the x and y wind components
As expected, there is no statistically significant relation between
Va and xw, yw for
the simulated random data (SRD, Table
1) when applying either a GAM
(Fig. 7A) or a linear model. On
the other hand, with the SWI dataset Va is significantly
related to xw and yw
(Fig. 7B,
Table 1). Both components of
wind are equally influential on Va. Va
increases as wind along the x-axis decreases or increases from 6 m
s-1 and decreases or increases from 3 m s-1 along the
y-axis. In other words as winds increasingly deviate from the ground
speed and direction (xg=6 m s-1,
yg=3 m s-1), Va increases.
For complete wind compensation, where Vg and
are
constant, the local minima of each variable indicate the mean ground vector
(xg, yg). Thus, by applying a GAM, we
find the true constant xg, yg
values.
|
|
For the observed autumn passerine migration (APM) the final model shows only a marginal linear influence of wind on Va (Table 1). This model was derived as follows. First a GAM was applied on the LOESS transformed wind components (xw, yw) (Fig. 7C). However, the LOESS transformations of both variables were not significant. A linear model was then fit for both variables. In this model, yw appeared to be insignificant, hence it was excluded.
In the final model, Va is only slightly influenced by wind along the x-axis. As winds blow more strongly towards the east, birds increase Va, whereas birds decrease Va when winds blow to the west (Fig. 7C). As tracks and headings are mainly towards SW (Fig. 5C), the x-component includes a tailwind as well as a side wind component. Although birds appear to increase air speed in headwinds and decrease in tailwinds when applying the conventional analysis, this relationship does not emerge when considering both wind components.
Conclusions
One question addressed by many biologists is: how do flying organisms adapt
their flight behavior to dynamic wind conditions? By way of an artificial
analytical example we have illustrated that the conventional approach to test
the hypothesis that birds maximize their distance per energy consumption by
increasing their airspeed in headwinds and decrease their airspeed in
tailwinds is incorrect. The negative relationship between
Va and Vg-Va can
result from spurious correlations in the data. The two reasons for this
spurious correlation are that: (1) Va is not only a
function of Vg-Va but also a function
of Vw and
-
; (2) Va
depends non-linearly on these variables; a linear model is only applicable
with constant ground speed or pure head/tailwind conditions. Alternative
general forms of analyses, that encompass the multidimensional and non-linear
nature of wind and flight, are necessary.
Our approach, a two-dimensional generalized additive model, provides a
simple, general and straightforward method to analyze the complex relation
between air speed and wind speed and direction without the risks entailed in
information reduction. In this study the new approach was tested on both
synthetic and observed data. Similarly, this approach can be applied to
studying the relationship between heading, rather then air speed, and both
wind components. Although other studies have described the importance of
multidimensional analysis (e.g. Green and
Alerstam, 2002
; Liechti et
al., 1994
) a general approach to studying the influence of wind
speed and direction simultaneously on air speed has not been suggested. See
however Shamoun-Baranes et al.
(Shamoun-Baranes et al.,
2003
), who applied a GAM to study the influence of the tailwind
and side wind component simultaneously on ground speed.
The analytical problems associated with the simplification of observational
data are a recurrent issue in biology (e.g.
Jackson, 1997
) and can have
serious consequences for the foundations of our theories. The broadly accepted
analysis of air speed and wind discussed in this paper has not only been used
to support the current theories on optimal avian flight but also to form them.
The validity of these and other one-dimensional studies in their current form
may therefore also influence our theoretical foundations. Hence, a reanalysis
of previous one-dimensional studies dealing with the effect of wind on air
speed is desirable. The adaptive behavior of birds towards their environment,
particularly wind, has important implications when incorporated into models
of, for example, stopover strategies and take-off decisions during migration
(Liechti and Bruderer, 1998
;
Weber et al., 1998
;
Weber and Hedenström,
2000
) estimations of potential flight range of long distance
migrants (Battley and Piersma,
2005
), of optimal flight speeds of birds
(Hedenström and Alerstam,
1995
), energetic requirements during for migration
(Butler et al., 1997
)
consequences for individual fitness (Clark
and Butler, 1999
) and the evolution of migratory strategies
(Erni et al., 2005
). In order
to properly interpret model results or compare models to measurements we must
ensure that the analysis underlying model assumptions or the predictions
themselves is appropriate.
We hope this paper stimulates new and revisited studies of the influence of wind on flight.



| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
Able, K. P. (1977). The flight behaviour of individual passerine nocturnal migrants: a tracking radar study. Anim. Behav. 25,924 -935.[CrossRef]
Alerstam, T. and Gudmundsson, G. A. (1999). Migration patterns of tundra birds: tracking radar observations along the northeast passage. Arctic 52,346 -371.
Alerstam, T., Gudmundsson, G. A. and Larsson, B. (1993). Flight tracks and speeds of Antarctic and Atlantic seabirds: radar and optical measurements. Philos. Trans. R. Soc. Lond. B Biol. Sci. 340,55 -67.[CrossRef]
Battley, P. F. and Piersma, T. (2005). Body composition and flight ranges of bar-tailed godwits (Limosa lapponica baueri) from New Zealand. Auk 122,922 -937.[CrossRef]
Butler, R. W., Williams, T. D., Warnock, N. and Bishop, M. (1997). Wind assistance a requirement for migration of shorebirds? Auk 114,456 -466.
Clark, C. W. and Butler, R. W. (1999). Fitness components of avian migration: a dynamic model of western sandpiper migration. Evol. Ecol. Res. 1,443 -457.
Cleveland, W. S. (1979). Robust locally weighted regression and smoothing scatterplots. J. Am. Stat. Assoc. 74,829 -836.[CrossRef]
Cleveland, W. S. and Devlin, S. J. (1988). Locally weighted regression: an approach to regression analysis by local fitting. J. Am. Stat. Assoc. 83,596 -610.[CrossRef]
Erni, B., Liechti, F. and Bruderer, B. (2005).
The role of wind in passerine autumn migration between Europe and Africa.
Behav. Ecol. 16,732
-740.
Green, M. and Alerstam, T. (2000). Flight speeds and climb rates of Brent geese: mass-dependent differences between spring and autumn migration. J. Avian Biol. 31,215 -225.[CrossRef]
Green, M. and Alerstam, T. (2002). The problem of estimating wind drift in migrating birds. J. Theor. Biol. 218,485 -496.[Medline]
Gudmundsson, G. A., Alerstam, T., Green, M. and Hedenström, A. (2002). Radar observations of Arctic bird migration at the Northwest Passage, Canada. Arctic 55, 21-43.
Guisan, A. and Zimmermann, N. E. (2000). Predictive habitat distribution models in ecology. Ecol. Modell. 135,147 -186.[CrossRef]
Guisan, A., Edwards, J., Thomas, C. and Hastie, T. (2002). Generalized linear and generalized additive models in studies of species distributions: setting the scene. Ecol. Modell. 157,89 -100.[CrossRef]
Hedenström, A. and Alerstam, T. (1995). Optimal flight speed of birds. Philos. Trans. R. Soc. Lond. Biol. Sci. 348,471 -487.
Hedenström, A. and Alerstam, T. (1996).
Skylark optimal flight speeds for flying nowhere and somewhere.
Behav. Ecol. 7,121
-126.
Hedenström, A., Rosen, M., Akesson, S. and Spina, F. (1999). Flight performance during hunting excursions in Eleonora's falcon Falco eleonorae. J. Exp. Biol. 202,2029 -2039.[Abstract]
Hedenström, A., Alerstam, T., Green, M. and Gudmundsson, G. A. (2002). Adaptive variation of airspeed in relation to wind, altitude and climb rate by migrating birds in the Arctic. Behav. Ecol. Sociobiol. 52,308 -317.[CrossRef]
Jackson, D. A. (1997). Compositional data in community ecology: the paradigm or peril of proportions? Ecology 78,929 -940.[CrossRef]
Liechti, F. (1993). Nachtlicher Vogelzug im Herbst uber Suddeutschland: Winddrift und Kompensation. J. Ornithol. 134,373 -404.[CrossRef]
Liechti, F. (1995). Modelling optimal heading and airspeed of migrating birds in relation to energy expenditure and wind influence. J. Avian Biol. 26,330 -336.
Liechti, F. (2006). Birds: blowin' by the wind? J. Ornithol. 147,202 -211.[CrossRef]
Liechti, F. and Bruderer, B. (1998). The relevance of wind for optimal migration theory. J. Avian Biol. 29,561 -568.
Liechti, F., Hedenström, A. and Alerstam, T. (1994). Effects of sidewinds on optimal flight speed of birds. J. Theor. Biol. 170,219 -225.[CrossRef]
Pennycuick, C. (1978). Fifteen testable predictions about bird flight. Oikos 30,165 -176.
Pennycuick, C. J. (1982). The flight of petrels and albatrosses (Procellariiformes), observed in South Georgia and its vicinity. Philos. Trans. R. Soc. Lond. B Biol. Sci. 300,75 -106.
Pennycuick, C. J. (1989). Bird Flight Performance: A Practical Calculation Manual. Oxford: Oxford University Press.
Pennycuick, C. J. (2001). Speeds and wingbeat
frequencies of migrating birds compared with calculated benchmarks.
J. Exp. Biol. 204,3283
-3294.
Rayner, J. M. V. (2001). Mathematical modelling of the avian flight power curve. Math. Meth. Appl. Sci. 24,1485 -1514.[CrossRef]
Riley, J. R., Greggers, U., Smith, A. D., Stach, S., Reynolds, D. R., Stollhoff, N., Brandt, R., Schaupp, F. and Menzel, R. (2003). The automatic pilot of honeybees. Proc. R. Soc. Lond. B Biol. Sci. 270,2421 -2424.[Medline]
Rosen, M. and Hedenström, A. (2001). Testing predictions from flight mechanical theory: a case study of Cory's shearwater and Audouin's gull. Acta Ethol. 3, 135-140.[CrossRef]
Shamoun-Baranes, J., Baharad, A., Alpert, P., Berthold, P., Yom-Tov, Y., Dvir, Y. and Leshem, Y. (2003). The effect of wind, season and latitude on the migration speed of white storks Ciconia ciconia, along the eastern migration route. J. Avian Biol. 34,97 -104.[CrossRef]
Srygley, R. B. (2003). Wind drift compensation in migrating dragonflies Pantala (Odonata: Libellulidae). J. Insect Behav. 16,217 -232.
Srygley, R. B. and Oliveira, E. G. (2001). Sun compass and wind drift compensation in migrating butterflies. J. Navig. 54,405 -417.[CrossRef]
Srygley, R. B., Oliveira, E. G. and Dudley, R. (1996). Wind drift compensation, flyways, and conservation of diurnal, migrant neotropical Lepidoptera. Proc. R. Soc. Lond. B Biol. Sci. 263,1351 -1357.[CrossRef]
Tucker, C. J. (1975). Flight energetics. Symp. Zool. Soc. Lond. 35, 49-63.
Wakeling, J. M. and Hodgson, J. (1992).
Optimization of the flight speed of the little, common and sandwich tern.
J. Exp. Biol. 169,261
-266.
Weber, T. P. and Hedenstrom, A. (2000). Optimal stopover decisions under wind influence: the effects of correlated winds. J. Theor. Biol. 205,95 -104.[CrossRef][Medline]
Weber, T. P., Alerstam, T. and Hedenstrom, A. (1998). Stopover decisions under wind influence. J. Avian Biol. 29,552 -560.
Wege, M. L. and Raveling, D. G. (1984). Flight speed and directional responses to wind by migrating Canada geese. Auk 101,342 -348.
This article has been cited by other articles:
![]() |
R. B. Srygley and R. Dudley Optimal strategies for insects migrating in the flight boundary layer: mechanisms and consequences Integr. Comp. Biol., July 1, 2008; 48(1): 119 - 133. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||