Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The Geography of Fear: A Latitudinal Gradient in Anti-Predator Escape Distances of Birds across Europe

  • Mario Díaz ,

    Mario.Diaz@ccma.csic.es

    Affiliation Department of Biogeography and Global Change (BGC), Museo Nacional de Ciencias Naturales, CSIC, Madrid, Spain

  • Anders Pape Møller,

    Affiliation Lab. d'Ecologie, Systématique et Evolution, CNRS UMR 8079, University Paris-Sud, Orsay, France

  • Einar Flensted-Jensen,

    Affiliation Cypresvej 1, DK-9700 Brønderslev, Denmark

  • Tomáš Grim,

    Affiliation Palacký University, Department of Zoology and Lab. Ornithology, Olomouc, Czech Republic

  • Juan Diego Ibáñez-Álamo,

    Affiliation Department of Zoología, Fac. Ciencias, Campus Universitario de Fuentenueva s/n, University Granada, Granada, Spain

  • Jukka Jokimäki,

    Affiliations University Lapland, Arctic Centre, Rovaniemi, Finland, Rovaniemi University of Applied Sciences, Rantavitikka Campus, Rovaniemi, Finland

  • Gábor Markó,

    Affiliations Eötvös Loránd University, Behavioral Ecology Group, Department of Systematics, Zoology and Ecology, Budapest, Hungary, Department of Plant Pathology, Corvinus University of Budapest, Budapest, Hungary

  • Piotr Tryjanowski

    Affiliation Inst. Zoology, Poznań University of Life Sciences, Poznan, Poland

Abstract

All animals flee from potential predators, and the distance at which this happens is optimized so the benefits from staying are balanced against the costs of flight. Because predator diversity and abundance decreases with increasing latitude, and differs between rural and urban areas, we should expect escape distance when a predator approached the individual to decrease with latitude and depend on urbanization. We measured the distance at which individual birds fled (flight initiation distance, FID, which represents a reliable and previously validated surrogate measure of response to predation risk) following a standardized protocol in nine pairs of rural and urban sites along a ca. 3000 km gradient from Southern Spain to Northern Finland during the breeding seasons 2009–2010. Raptor abundance was estimated by means of standard point counts at the same sites where FID information was recorded. Data on body mass and phylogenetic relationships among bird species sampled were extracted from the literature. An analysis of 12,495 flight distances of 714 populations of 159 species showed that mean FID decreased with increasing latitude after accounting for body size and phylogenetic effects. This decrease was paralleled by a similar cline in an index of the abundance of raptors. Urban populations had consistently shorter FIDs, supporting previous findings. The difference between rural and urban habitats decreased with increasing latitude, also paralleling raptor abundance trends. Overall, the latitudinal gradient in bird fear was explained by raptor abundance gradients, with additional small effects of latitude and intermediate effects of habitat. This study provides the first empirical documentation of a latitudinal trend in anti-predator behavior, which correlated positively with a similar trend in the abundance of predators.

Introduction

The global distribution of biodiversity reflects higher degrees of specialization and interspecific interactions at low latitudes. There is a long tradition for linking these trends by viewing high tropical biodiversity as a consequence of interspecific interactions [1], [2]. For example, MacArthur [3] suggested explicitly that southern distribution limits of many species in the northern hemisphere likely depended on biotic interactions, while northern distributions more often were determined by abiotic factors. Studies of latitudinal trends in species diversity have shown consistent decreases with increasing latitude that is also repeated in predators, parasites and numerous other functional groups [e.g. 4–6]. Interestingly, not all groups show the same slope of species richness with latitude, implying that the abundance and species composition of different functional groups, and hence the intensity of interspecific interactions such as predation or competition, may change with latitude (reviewed by [7]).

Levels of defense are predicted to coevolve with levels of offense along productivity gradients such as the gradient from the arctic to the tropics [8]. Although predation constitutes an important interspecific interaction, relatively few studies have investigated latitudinal patterns in predation and anti-predator behavior [7]. Several studies have shown latitudinal changes in palatability [9] or chemical or mechanical defenses of prey [10], [11]. However, many studies relied on a sample of individuals from a single temperate and a single tropical site, thereby offering no replication (i.e., a case of pseudo-replication) and thus increasing the risk that between population differences were caused by factors unrelated to geographical position (e.g. [12]; reviewed in [7]). Furthermore, the generality of latitudinal gradients across taxa implies that processes underlying them should have evolved repeatedly through convergence. Recently, McKinnon et al. [13] showed a northwards decrease in nest predation rates of artificial nests along a latitudinal gradient of more than 3000 km in North America. Likewise, frogs Rana temporaria showed a strong latitudinal cline in anti-predator behavior in Sweden [14], and Møller and Liang [15] interpreted higher levels of fear (longer flight initiation distances -hereafter: FID- when approached by a researcher; see below) in Chinese as compared to European birds partly due to differences in predation risk. Although these studies provide important information, there is a need to extend them to include multiple model species along wide latitudinal gradients to test for the generality of patterns (i.e., meta-replication [16], [17]). Furthermore, it is also necessary to test the underlying assumption that differences in anti-predator behavior are related to differences in predation risk [18].

Animals in general show a range of fear responses that have important implications for understanding habitat use, but also life history, demography, interspecific interactions and conservation [18][21]. Flight initiation distance (hereafter FID) when approached by potential predators and other kinds of anti-predator behavior pose a problem of optimization of benefits associated with escape weighed against costs of disturbance (like stop feeding and/or energetic costs of flight), which may vary with ambient conditions. In general, FID constitutes a reliable measure of risk taking because it correlates with susceptibility to predation, even when the approaching potential predator was the researcher who measured such distances (see validation of flight distances as a meaningful proxy of risk taking in [22][24]). In addition, FID is highly repeatable for individuals among estimates [25], [26].

Another important aspect of risk taking is the effect of urbanization, especially given that several studies have found differences in predator densities between urban and non-urban habitats [27][31]. In fact, birds have considerably shorter FIDs in urban than in rural habitats [32], [33], mainly linked to differences in risk of predation and duration since urbanization [26]. These behavioral responses were largely independent of differences in habitat structure between urban and nearby rural sites [26], [34]. European species of birds with long FIDs have negative population trends as expected if there are costs associated with frequent disturbance by humans, dogs and other potential predators [34]. Thus, it is crucial to take this urbanization effect into account when investigating latitudinal gradients in FID and its relationships with predation risk.

The objective of this study was to test for the existence and predictors of a latitudinal trend in FID of birds when being approached by a human, taken here as a methodologically proper surrogate for the risk of predation experienced by bird populations. We studied FIDs at nine different pairs of nearby rural and urban sites along a latitudinal gradient from Southern Spain to Northern Finland (Figure 1), allowing for repeated tests among independent populations within species. Bird census data from these sites [35] were used to test whether raptor abundance showed a similar latitudinal gradient. Because FID differs consistently between rural and urban populations [32], [33], we also tested if the latitudinal trend was consistent independently of habitat. Finally, we tested whether latitudinal gradients in FID were related to parallel gradients in raptor abundance. Overall, our main aim was to test whether different bird species showed consistent latitudinal patterns of anti-predator behavior in relation to spatial patterns of predation risk. If that was the case, this would be the first empirical documentation of such a latitudinal trend.

thumbnail
Figure 1. Location of the nine study sites for FID in Europe.

Major climatic zones (subtropical, temperate and continental-subarctic, from darker to lighter shading) are also shown. For latitudes and city names see Table S1.

https://doi.org/10.1371/journal.pone.0064634.g001

Methods

Ethic Statement

Recording of flight initiation distance does not require ethics approval or legal permits because it does not involve any important effect on animal welfare. The disturbance produced to birds by our methodology did not differ from standard “background” disturbance caused by humans: birds in our study sites, especially in the urban ones, are regularly flushed by citizens, dogs, cats, cars, etc., and by tourists, pets, agricultural workers and mushroom pickers in rural areas. Moreover, our presence in study sites involved no alteration of the habitat and it was very brief in time. We used public trails and roads, thus, there was no need to ask land managers for approval.

Study Areas

Fieldwork was performed in nine cities, each paired with a nearby rural area, located along a large latitudinal gradient across Europe (Fig. 1). The distance between urban and rural study sites was 1–20 km. The benefits of this approach are two-fold: (1) neighboring study sites will share most potentially confounding environmental characteristics including weather, altitude, soil and many others; and (2) birds will not be prevented from moving between neighboring urban and rural habitats for distance reasons, although other factors such as philopatry and assortative mating may prevent such movement. All urban study sites included areas with multi-storey buildings, single family houses, roads and parks, while nearby rural areas had open farmland and woodland and did not contain continuous urban elements like multi-storey buildings, one family houses, roads and parks. This simple operational definition was also adopted in other studies [36][39], and our definitions of urban (percent of built-up area >50, building density >10/ha and residential human density >10/ha) and rural habitats (percent of built-up area 5–20, building density 2.5–10/ha, and residential human density 1–10/ha) follow the suggestion made by Marzluff et al. [39].

Flight Initiation Distance

We recorded FIDs of birds when approached by a human during the breeding seasons 2009–2010 using a modified technique of that used by Blumstein [21]. A full description of the procedures and three different cross-validations of the data are reported elsewhere [24], [33], [34]. We tried to obtain data from as many species and individuals as possible by systematic searches of the study areas; at the same time, we tried to avoid sampling the same individual twice by moving to another site right after samples were taken. If the same general area was visited, only individuals of different species, sex or age than those sampled before were tested. In brief, when an individual bird had been located with a pair of binoculars, an observer moved at a normal walking speed towards the individual bird, while recording the number of 1 m steps [22]. The distance from the observer to the bird when it first took flight was recorded as the FID, while the ‘starting distance’ was the distance from where the observer started walking up towards the bird and the location of the bird. If the individual bird was positioned in the vegetation, the height above ground was recorded to the nearest meter (birds on the ground were assigned a height of zero). While recording these flight initiation distances, we also recorded date and time of day, and the sex and the age of the individual if external characteristics allowed sexing and aging with binoculars. Starting distances and FIDs were estimated as the Euclidian distances that equals the square root of the sum of the squared horizontal distances and the squared height above ground level [21].

Previous studies have shown that starting distance is strongly positively correlated with FID [21], thereby causing a problem of collinearity. We eliminated this problem of collinearity by searching habitats for birds with a pair of binoculars when choosing an individual for estimating flight initiation distance. In this way we assured that most individuals were approached from a distance of at least 30 m, thereby keeping starting distances constant across species.

We extracted information on mean body mass of adult birds of each species from Cramp & Perrins [40], as FIDs are consistently affected by bird size [27]. FID was weakly negatively related to starting distance in a Generalized Linear Model that included species, age, habitat, country and body mass as factors (partial F = 34.99, df = 1, 12493, P<0.0001; FID and body mass log-transformed), explaining less than 0.3% of the variance (see [26] for a similar procedure). None of the results presented in this paper changed statistically when including starting distance as an additional variable, and we thus excluded this variable from all subsequent analyses for simplicity.

FIDs varied significantly among species, using one-way analysis of variance with log10-transformed flight distance as the response variable and species as a factor (F = 29.58, df = 158, 12337, R2 = 0.33, P<0.0001, repeatability (R) following Becker [41] R = 0.27, SE = 0.03). Thus there was consistency in flight distance within species.

Raptor Abundance

The breeding bird community at study sites was censused by means of standard point counts with unlimited recording distance [42]. Censuses were made during the springs 2009–2010 in both urban and rural habitats in all study locations at the same sites where FID information was recorded [35]. From these data we derived an index of abundance of raptors (mostly hawks Accipiter spp., kestrels Falco spp. and buzzards Buteo buteo) by summing the numbers of individuals detected at each site divided by the number of sampling points (data in [35]). Briefly, we placed 25–50 points in each urban and rural site at distances of at least 100 m between two consecutive points. The exact location of each point was determined with a GPS, allowing us to make the second census in exactly the same sites as the first census. The first census was made in early April in Southern Spain, delaying the census at higher latitudes so it was completed in Northern Finland in late May, and the second census was carried out three-four weeks later. The census started at local sunrise, while remaining five minutes at each point recording all birds seen or heard. Censuses started on separate days in urban and rural study plots ensuring that there was no difference in timing of censuses between habitats. The same observer made all surveys in each particular city and its paired rural site. Vegetation cover (trees, shrubs, herbs and grass) and cover with buildings and other man-made structures were evaluated in the field within 50 m of each survey point. We obtained very similar abundance indices when controlling or not controlling for differences in coverage for the three vegetation layers [35].

Statistical Analyses

We log10-transformed FID, body mass and raptor abundance to achieve distributions that did not differ from normality. Differences in raptor abundance along the latitudinal gradient and between urban and rural sites were tested by means of General Linear Models (GLM). Relationships between FID and predictor (latitude, urban-rural habitat, and raptor abundance) and confounding (body size) variables were also tested by a GLM approach using the mean FID obtained for all individuals of the same species sampled in each study site. However, taxonomic units such as species cannot be considered statistically independent observations due to effects of common ancestry [43]. As species occupy a variable number of study sites [35], latitudinal or habitat effects could be partly due to phylogenetic effects mediated by species composition. To control for the phylogenetic relationship among the sampled species we used phylogenetic generalized least square regression (PGLS) models as implemented in R statistical environment (see [44], [45] for a similar approach). We used the R libraries ape, MASS and mvtnorm and the function pglm3.1.r. First, we estimated the phylogenetic scaling parameter lambda (λ), that varies between 0 (phylogenetic independence) and 1 (species’ traits covary in direct proportion to their shared evolutionary history; [43]). Then, we calculated the phylogenetically corrected correlation between the variables of interest after adjusting for phylogenetic effects through the estimated λ. We have used the most recent and comprehensive bird phylogeny available [46], after editing it to include different populations of the same species as polytomies with a constant small genetic distance of 1·10−10 between conspecific populations. This phylogeny considers Corvus corone and C. cornix as conspecifics, so that we have also assigned a genetic distance between them of 1·10−10. The R script and the edited phylogeny file are supplied as Files S1 and S2, respectively. Sequential (type I) models were fitted, including first the effect of body size to test for effects of latitude and habitat, and the sequential effects of body size and raptor abundance to ascertain whether latitudinal and habitat effects changed when raptor abundance effects were included first.

We evaluated the magnitude of associations between FID and predictor variables based on effect sizes and its 95% CI, as computed from P values of t-tests according to Cohen [47]. Cohen’s criteria were small (Pearson r <0.10, explaining <1% of the variance), intermediate (r = 0.11–0.49, 9–24% of the variance) and large (r >0.50, 25% or more of the variance) effect sizes.

Results

Summary Statistics for FID and Latitude

We obtained information on FID from 714 populations belonging to 159 species, with the total number of individual estimates of flight distance being 12,495 (Table S1). The number of estimates per species ranged from 1 to 1474, mean (SE) = 79 (9). There were 430 estimates from rural and 284 from urban populations. Latitude ranged from 37°N to 66°N, or 29°×111 km per degree = 3,220 km, with climate ranging from subtropical to subarctic. Mean FID in the 714 populations ranged from 1 to 180 m, mean (SE) = 11.46 m (1.03; back-transformed from log10-transformed data), with a highly significant effect of species (F = 3.06, df = 158, 555, R2 = 0.47, P<0.0001).

Raptor Density and Latitude

The abundance of raptors decreased significantly with increasing latitude when analyzing abundance at point counts (GLM; F = 12.08, df = 1, 2425, P = 0.0005) and was smaller in urban than in rural habitats (F = 5.51, df = 1, 2425, P = 0.019; least square mean for rural habitats: 0.0659 birds/census point (SE = 0.0080); least square mean for urban habitats: 0.0397 birds/census point (0.0077)). Mean density of raptors at the nine sites showed a significant trend with decreasing abundance at higher latitudes (Fig. 2; GLM: F = 5.67, df = 1, 7, P = 0.032) with a large effect size (r = 0.506; 95% CI: 0.052–0.787). This decreasing trend did not differ among paired urban and rural sites (F = 1.19, df = 1, 7, P = 0.293; Fig. 2).

thumbnail
Figure 2. Mean (± SE) index of abundance of raptors in relation to latitude.

Data are for paired rural (closed circles, continuous line) and urban (open circles, pointed line) sites. The index of abundance was estimated as the mean number of individuals detected at 25–50 point counts (depending on urban area size) each lasting five minutes. Lines are linear regression lines based on mean values. Symbols have been slightly displaced above (urban) and below (rural) real latitudes to improve clarity.

https://doi.org/10.1371/journal.pone.0064634.g002

FID and Latitude

Log-transformed mean FID of bird populations decreased with increasing latitude (Fig. 3), with an intermediate effect size (ca. 5%) after accounting for significant positive effect of body mass on FID as well as for phylogenetic effects (Table 1). Thus, anti-predator behavior decreased in intensity at higher latitudes. There was a highly significant difference in mean FID between populations of the same species from rural and urban habitats (mean (SE) rural = 13.34 m (1.05); urban = 7.12 m (1.05); paired t-test, t = −13.99, df = 229, P<0.0001), with a large effect size of 27% after accounting for body mass and phylogeny (Table 1). Finally, the difference in mean FID between rural and urban habitats decreased with increasing latitude (Fig. 3), with a small effect size (0.8%) after adjusting for the effect of body mass and phylogeny (Table 1).

thumbnail
Figure 3. Mean FID (m) in relation to latitude for different populations and species of birds.

Data are for paired rural (closed circles, continuous line) and urban (open circles, pointed line) sites. Lines are linear regression lines. Symbols have been slightly displaced above (urban) and below (rural) real latitudes to improve clarity.

https://doi.org/10.1371/journal.pone.0064634.g003

thumbnail
Table 1. Flight initiation distance in relation to latitude, habitat (coded as 0 = rural or 1 = urban) and the interaction between latitude and habitat in European birds, after correcting for phylogenetic relationships among the sampled species by means of phylogenetic generalized least square regression (PGLS) models and for body size effects by mean of sequential (Type I) Sum of Squares (SS).

https://doi.org/10.1371/journal.pone.0064634.t001

Log-transformed mean FID of bird populations increased with increasing raptor abundance (F1,712 = 76.79, P<0.0001), with an intermediate effect size (ca 10%; r = 0.315; 95% CI: 0.247–0.380), after accounting for phylogenetic and body size effects by means of a sequential PGLS model including body size and raptor abundance. This model had the statistics F2,711 = 49.54, Adjusted R2 = 0.12, P<0.0001, AIC = 301.99, and phylogenetic scaling parameter (λ) = 0.472, P<0.0001. Latitudinal and habitat effects were still significant when accounting for the effects of raptor abundance (Table 2), but effect sizes decreased to small and intermediate, respectively (1.4 and 24%, respectively). The model including raptor abundance, latitude and habitat fitted significantly better to the data than the model not including this factor (ΔAIC = 113.03–106.84 = 6.19; Tables 1 and 2) and the model including only raptor abundance effects (ΔAIC = 301.99–106.84 = 195.15).

thumbnail
Table 2. Flight initiation distance in relation to latitude, habitat (coded as 0 = rural or 1 = urban) and the interaction between latitude and habitat in European birds, after correcting for phylogenetic relationships among the sampled species by means of phylogenetic generalized least square regression (PGLS) models, and for body size and raptor abundance effects by means of sequential (Type I) Sum of Squares (SS).

https://doi.org/10.1371/journal.pone.0064634.t002

Discussion

We have shown that flight initiation distance (FID) of birds when approached by a human, which constitutes an important component of anti-predator behavior, showed a clear-cut latitudinal gradient in Europe, after accounting statistically for factors influencing FID such as body size, habitat (urban or rural) or phylogenetic effects. This behavioral gradient was consistent among bird species and matched partially local changes in avian predator abundance. Overall, our study thus provides the first empirical documentation of a latitudinal trend in predator-prey interactions. Such a gradient may potentially influence the composition and structure of prey communities [7].

The latitudinal trend in anti-predator behavior paralleled latitudinal changes in raptor abundance, and varied depending on whether the study population was urban or rural, with smaller differences between rural and urban habitats at high latitudes, where overall raptor abundance was smaller. These findings show that anti-predator behavior consistently varied with latitude, apparently in response to latitudinal differences in predation risk [12]. Additionally, they highlight that the observed changes in anti-predator behavior detected in other studies due to the process of urbanization [30], [31] could depend on latitude too. However, predation risk, as measured by local estimates of raptor abundance, did not fully explain latitudinal and habitat differences in anti-predator behavior, as the most parsimonious model also included additional small effects of latitude per se and intermediate effects of habitat. In fact, the relationship between predator abundance and FID seemed to hold mainly for rural habitats along the sampled latitudinal gradient, whereas latitudinal changes in both FID and raptor abundance were much smaller in nearby urban sites (Figs. 2 and 3). We suggest that effects of other predators may account for additional variance.

Urban populations had consistently shorter flight distances than nearby rural populations, as found in previous studies of birds [32], [33], with the mean estimate being 7.1 m for urban birds and 13.3 m for rural birds, or an almost two-fold difference. There was an independent positive relationship between flight distance and body mass, with a large effect size after accounting for phylogenetic effects (Table 1). While many traits are correlated with body mass, thus making interpretation of this effect difficult, a greater difference in flight distance in large species may be attributed to smaller local population sizes, larger founder effects and elevated probability of local adaptation to urban conditions. Alternatively, large species have longer flight distances than small species [33], [48], and such large flight distances may be particularly energetically costly due to the disproportionately high costs of short flights that may reach a level 23 times higher than basal metabolic rate [49]. Thus, selection for adaptation to urban conditions will be particularly intense in large species due to frequent disturbance and flight, and relatively large reductions in flight distance of urban compared to rural populations in large species are consistent with this interpretation. Incidentally, this body size effect contributes to explain why urban areas are safer for prey, as avian predators are more affected by human disturbance in cities than their smaller-sized prey [48].

We predicted that interspecific interactions as reflected by the abundance of raptors are more important at southern latitudes, as the relative importance of abiotic (i.e. climatic) vs. biotic factors (i.e. competition, predation) on prey would decrease at lower latitudes. The observation that differences in FID between rural and urban habitats decreased towards the north (in the northern hemisphere) supports this prediction. Thus, our findings suggest that differences in abundance of raptors contributed to explain the difference in FID between rural and urban habitats along the latitudinal gradient across Europe. These findings are also consistent with the original suggestions by Janzen, Connell and MacArthur [1][3] that biotic interactions are important at lower latitudes, while abiotic factors are more important at higher latitudes.

In conclusion, we have shown that flight distance in birds, when approached by a human (an important component of anti-predator behavior) decreased with latitude, although varying among species and between rural and urban sites locally. Differences in flight distance between rural and urban habitats decreased with increasing latitude Overall, we have demonstrated for the first time a consistent spatial pattern of anti-predator behavior over a wide latitudinal and land-use gradient, with potential relevant effects on prey communities with different susceptibilities to predation.

Supporting Information

Table S1.

Summary statistics. Information on bird species and populations, study cities and countries, mean FID (m), SD FID (m), sample size (N), habitat: urban (1) or rural (0), latitude of study site (° N) and body mass (g) for the 714 sampled populations. See Material and methods for sources. Nomenclature follows [46] and population codes match those in file S2, where such populations are included as polytomies with a distance of 1·10−10.

https://doi.org/10.1371/journal.pone.0064634.s001

(DOC)

File S1.

R script for the phylogenetic generalized least square regression (PGLS) models. They were fitted to analyze relationships among FID, body mass, raptor abundance, latitude and habitat.

https://doi.org/10.1371/journal.pone.0064634.s002

(DOCX)

File S2.

Phylogeny of all European bird species edited to include the sampled populations as polytomies. Polytomies were indicated with a constant small genetic distance of 1·10−10 between conspecific populations.

https://doi.org/10.1371/journal.pone.0064634.s003

(TXT)

Acknowledgments

J. J. Cuervo and J.J. Soler helped us with the phylogenetic analyses. Comments by D. Blumstein and an anonymous referee were very helpful. This paper is a contribution by MD to the project MONTES (CSD2008–00040) and RiskDisp (CGL2009–08430) of the Spanish Ministerio de Ciencia e Innovación and to the thematic network GlobiMed (www.globimed.net).

Author Contributions

Conceived and designed the experiments: MD APM. Performed the experiments: MD APM EFJ TG JDIA JJ GM PT. Analyzed the data: MD APM. Contributed reagents/materials/analysis tools: MD APM EFJ TG JDIA JJ GM PT. Wrote the paper: MD APM EFJ TG JDIA JJ GM PT.

References

  1. 1. Janzen DH (1970) Herbivores and the number of tree species in tropical forests. American Naturalist 104: 501–528.
  2. 2. Connell JH (1971) On the role of natural enemies in preventing competitive exclusion in some marine animals and in rain forest trees. In: den Boer PJ, Gradwell GR editors. Dynamics of numbers in populations. Wageningen, The Netherlands: Centre for Agricultural Publication and Documentation, 298–312.
  3. 3. MacArthur RH (1972) Geographical ecology. New York: Harper & Row.
  4. 4. Rohde K (1998) Latitudinal gradients in species diversity. Area matters, but how much? Oikos 82: 184–190.
  5. 5. Guégan J-F, Thomas F, Hochberg ME, de Meeus T, Renaud F (2003) Disease diversity and human fertility. Evolution 55: 1308–1314.
  6. 6. Guernier V, Hochberg ME, Guégan J-F (2004) Ecology drives the worldwide distribution of human diseases. PLoS Biology 2: 740–746.
  7. 7. Schemske DW, Mittelbach GG, Cornell HV, Sobel JM, Roy K (2009) Is there a latitudinal gradient in the importance of biotic interactions? Annual Review of Ecology, Evolution and Systematics 40: 245–269.
  8. 8. Hochberg ME, Van Baalen M (1998) Antagonistic coevolution over productivity gradients. American Naturalist 152: 620–634.
  9. 9. Bolser RC, Hay ME (1996) Are tropical plants better defended? Palatability and defenses of temperate vs tropical seaweeds. Ecology 77: 2269–2286.
  10. 10. Vermeij GJ (1978) Biogeography and adaptation: Patterns of marine life. Cambridge: Harvard University Press.
  11. 11. Dyer LA, Coley PD (2002) Tritrophic interactions in tropical and temperate communities. In: T Tscharntke, B Hawkins, editors. Multitrophic level interactions. Cambridge: Cambridge University Press, 67–88.
  12. 12. Jokimäki J, Kaisanlahti-Jokimäki ML, Sorace A, Fernández-Juricic E, Rodríguez-Prieto I, et al. (2005) Evaluation of the "safe nesting zone" hypothesis across an urban gradient: a multi-scale study. Ecography 28: 59–70.
  13. 13. McKinnon L, Smith PA, Nol E, Martin JL, Doyle FI, et al. (2010) Lower predation risk for migratory birds at high latitudes. Science 327: 326–327.
  14. 14. Laurila A, Lindgren B, Laugen AT (2008) Antipredator defenses along a latitudinal gradient in Rana temporaria. Ecology 89: 1399–1413.
  15. 15. Møller AP, Liang W (2013) Tropical birds take small risks. Behavioral Ecology 24: 267–272.
  16. 16. Johnson DH (2002) The importance of replication in wildlife research. Journal of Wildlife Management 66: 919–932.
  17. 17. Grim T, Samaš P, Moskát C, Kleven O, Honza M, et al. (2011) Constraints on host choice: why do parasitic birds rarely exploit some common potential hosts? Journal of Animal Ecology 80: 508–518.
  18. 18. Zanette LY, White AF, Allen MC, Clinchy M (2011) Perceived predation risk reduces the number of offspring songbirds produced per year. Science 334: 1398–1401.
  19. 19. Ydenberg RC, Dill LM (1986) The economics of fleeing from predators. Advances in the Study of Behaviour 16: 229–249.
  20. 20. Caro T (2005) Antipredator defenses in birds and mammals. Chicago: University of Chicago Press.
  21. 21. Blumstein DT (2006) Developing an evolutionary ecology of fear: How life history and natural history traits affect disturbance tolerance in birds. Animal Behavior 71: 389–399.
  22. 22. Møller AP, Nielsen JT, Garamszegi LZ (2008) Risk taking by singing males. Behavioral Ecology 19: 41–53.
  23. 23. Møller AP, Erritzøe J, Nielsen JT (2010a) Causes of interspecific variation in susceptibility to cat predation on birds. Chinese Birds 1: 1–15.
  24. 24. Møller AP (2008c) Flight distance and blood parasites in birds. Behavioral Ecology 19: 1305–1313.
  25. 25. Carrete M, Tella JL (2010) Individual consistency in flight initiation distances in burrowing owls: a new hypothesis on disturbance-induced habitat selection. Biology Letters 23: 167–170.
  26. 26. Møller AP (2010) Interspecific variation in fear responses predicts urbanization in birds. Behavioral Ecology 21: 365–371.
  27. 27. Jerzak L (2001) Synurbanization of the Magpie in the Palearctic. In: Marzluff JM, Bowman R, Donelly R, editors. Avian Ecology and Conservation in an Urbanizing World. New York: Kluwer, 403–425.
  28. 28. Lepczyk CA, Mertig AG, Liu J (2003) Landowners and cat predation across rural-to-urban landscapes. Biological Conservation 115: 191–201.
  29. 29. Sorace A, Gustin M (2009) Distribution of generalist and specialist predators along urban gradients. Landscape and Urban Planning 90: 111–118.
  30. 30. Møller AP (2011) Song post height in relation to predator diversity and urbanization. Ethology 117: 529–538.
  31. 31. Møller AP, Ibáñez Álamo JD (2012) Escape behaviour of birds provides evidence of predation being involved in urbanization. Animal Behavior 84: 341–348.
  32. 32. Cooke AS (1980) Observations on how close certain passerine species will tolerate an approaching human in rural and suburban areas. Biological Conservation 18: 85–88.
  33. 33. Møller AP (2008a) Flight distance of urban birds, predation and selection for urban life. Behavioral Ecology and Sociobiology 63: 63–75.
  34. 34. Møller AP (2008b) Flight distance and population trends in European breeding birds. Behavioral Ecology 19: 1095–1102.
  35. 35. Møller AP, Díaz M, Flensted-Jensen E, Grim T, Ibáñez-Álamo JD, et al. (2012) High urban population density of birds reflects their timing of urbanization. Oecologia 170: 867–875.
  36. 36. Klausnitzer B (1989) Verstädterung von Tieren. Wittenberg Lutherstadt: Neue Brehm-Bücherei.
  37. 37. Gliwicz J, Goszczynski J, Luniak M (1994) Characteristic features of animal populations under synurbanization: The case of the Blackbirds and the striped field mouse. Memorabilia Zoologica 49: 237–244.
  38. 38. Stephan B (1990) Die Amsel. Wittenberg-Lutherstadt: Neue Brehm-Bücherei.
  39. 39. Marzluff JM, Bowman R, Donnelly R (2001) A historical perspective on urban bird research: trend, terms, and approaches. In: Marzluff JM, Bowman R, Donnelly R, editors. Avian Ecology and Conservation in an Urbanizing World. New York: Kluwer, 20–47.
  40. 40. Cramp S, Perrins CM, editors (1977–1994) The birds of the Western Palearctic. Vols. 1–9. Oxford: Oxford University Press.
  41. 41. Becker WA (1984) Manual of quantitative genetics. Pullman: Academic Enterprises.
  42. 42. Voříšek P, Klvaňová A, Wotton S, Gregory RD (2010) A best practice guide for wild bird monitoring schemes. Brussels: European Union.
  43. 43. Freckleton RP, Harvey PH, Pagel M (2002) Phylogenetic analysis and comparative data: A test and review of evidence. The American Naturalist 160, 712–726.
  44. 44. Garamszegi LZ, Møller AP (2007) Prevalence of avian influenza and host ecology. Proceedings of the Royal Society B 274: 2003–2012.
  45. 45. Soler JJ, Avilés JM (2010) Sibling competition and conspicuousness of nestling gapes in altricial birds: A comparative study. PLoS ONE 5(5): e10509
  46. 46. Thuiller W, Lavergne S, Roquet C, Boulangeat I, Lafourcade B, et al. (2011) Consequences of climate change on the tree of life in Europe. Nature 470: 531–534.
  47. 47. Cohen J (1988) Statistical power analysis for the behavioral sciences. 2nd edn. Hillsdale: Lawrence Erlbaum.
  48. 48. Møller AP (2012) Urban areas as refuges from predators and flight distance of prey. Behavioral Ecology 23: 1030–1035.
  49. 49. Tatner P, Bryant DM (1986) Flight cost of a small passerine measured using doubly labeled water: implications for energetic studies. The Auk 103: 169–180.
  50. 50. Lipsey MW, Wilson DB (2001) Practical Meta-analysis. Thousand Oaks: SAGE Publications Inc. George mason University website. Available: http://gemini.gmu.edu/cebcp/EffectSizeCalculator/index.html. Accessed 2013 Apr 26.