Skip to main content

To what extent does Tobler's 1st law of geography apply to macroecology? A case study using American palms (Arecaceae)

Abstract

Background

Tobler's first law of geography, 'Everything is related to everything else, but near things are more related than distant things' also applies to biological systems as illustrated by a general and strong occurrence of geographic distance decay in ecological community similarity. Using American palms (Arecaceae) as an example, we assess the extent to which Tobler's first law applies to species richness and species composition, two fundamental aspects of ecological community structure. To shed light on the mechanisms driving distance decays in community structure, we also quantify the relative contribution of geographic distance per se and environmental changes as drivers of spatial turnover in species richness and composition.

Results

Across the Americas, similarity in species composition followed a negative exponential decay curve, while similarity in species richness exhibited a parabolic relationship with geographic distance. Within the four subregions geographic distance decays were observed in both species composition and richness, though the decays were less regular for species richness than for species composition. Similarity in species composition showed a faster, more consistent decay with distance than similarity in species richness, both across the Americas and within the subregions. At both spatial extents, geographic distance decay in species richness depended more on environmental distance than on geographic distance, while the opposite was true for species composition. The environmentally complex or geographically fragmented subregions exhibited stronger distance decays than the more homogenous subregions.

Conclusion

Similarity in species composition exhibited a strong geographic distance decay, in agreement with Tobler's first law of geography. In contrast, similarity in species richness did not exhibit a consistent distance decay, especially not at distances >4000 kilometers. Therefore, the degree to which Tobler's first law of geography applies to community structure depends on which aspect hereof is considered – species composition or species richness. Environmentally complex or geographically fragmented regions exhibited the strongest distance decays. We conclude that Tobler's law may be most applicable when dispersal is a strong determinant of spatial turnover and less so when environmental control predominates.

Background

Tobler's first law of geography, 'Everything is related to everything else, but near things are more related than distant things' [1] (see review in [3]; hereafter referred to as Tobler's law), was first applied to urban growth systems, but it also applies to biological systems as illustrated by a general occurrence of distance decays in ecological community similarity [2]. Its applicability to ecology is closely related to key theoretical issues such as what determines species diversity [4] and the distribution and abundance of species [51], as well as central to the way analyses in ecology are performed [5, 65]. A negative relationship between community similarity and geographic distance is often attributed to environmental gradients [2, 20]. However, the 300-years old observation that environmentally similar, but non-contiguous regions harbour distinct assemblages of vertebrates and plants (Buffon's law or 'the first principle of biogeography' [6]) suggests that other factors play a role, too. Traditional explanations have emphasized dispersal limitation due to geographic barriers [20], but spatially limited dispersal can generate distance decays in community similarity even in the absence of barriers [7, 8]. A negative relationship is therefore expected between community similarity and geographic distance not only as a consequence of environmental gradients, but also due to dispersal limitation [7–9]. The latter notion is strongly contrasted by the view that 'everything is everywhere, but the environment selects' (Baas-Becking's or Beijerick's law), which suggests that dispersal limitation is unimportant [10, 11]. At a global scale, this view clearly does not apply to larger organisms, as epitomized in Buffon's law. Nevertheless, it has often been argued that species distributions are largely in equilibrium with environmental conditions within continents or smaller regions [12, 13]. The issue is controversial, however [14], and other authors have emphasized the role of non-environmental range constraints [16], notably dispersal limitation [8, 15].

When applied to ecological communities, Tobler's law has been used to refer to community similarity in terms of species composition, but communities are characterized by many other features, e.g., species richness. Large-scale variability in species richness is often argued to largely depend on climate [21, 22], but many competing explanations exist [15, 21, 23–30]. Therefore, it becomes relevant to ask whether Tobler's law can be extended to also cover other macroecological features such as community similarity in terms of species richness and to understand the underlying drivers as well.

Here, we use American palms to test the applicability of Tobler's law to macroecology. Palms are common in warm parts of the New World [31–33], and are particularly species-rich close to the equator [34]. Climatic water-related factors appear to be a major control of palm species richness patterns in the Americas, but nonetheless there are also historical and unexplained broad-scale spatial patterns [34, 35]. Previous studies of distance decays in palm species composition have focused on local to regional scales [36, 37]. In this study, we use distribution data on palm species richness and composition across the Americas to investigate the general applicability of Tobler's law to palm macroecology. Specifically, to obtain a deeper understanding of the mechanisms controlling distance decays in similarity of species composition and richness, we assess the following three key hypotheses: (1) If species composition is more strongly influenced by dispersal limitation than species richness, a stronger, more regular distance decay is expected for similarity in species composition. (2) As a further corollary, geographic distance will have a stronger impact than environmental distance on the distance decay in similarity in species composition, whereas the opposite will be true for species richness. (3) Comparing different regions within the Americas (Figure 1), the strength of the distance decay in community similarity will be positively correlated with the heterogeneity and complexity of the region, i.e., strongest in environmentally complex (e.g., mountainous regions) or geographically fragmented regions (e.g., island archipelagos). The former may reflect either the direct effect of the environmental gradients or the many barriers to dispersal in environmentally complex regions, while the latter more unambiguously reflect limited dispersal.

Figure 1
figure 1

Similarity of palm species richness and composition across the Americas. Data compiled in 1° × 1° grid cells across the Americas. The four smaller subregions Amazon, Andes, Caribbean and Central America are marked.

Results

Distance decay in palm species richness and composition

The distance decay for palm species richness is weaker and less consistent than the decay for palm species composition across the Americas. The similarity of species richness declines over the first 4000 kilometers, but then increases again (Fig. 2), reflecting that species richness is high in the central, equatorial part of the Americas and low towards the northern and southern limits of our study area (Fig. 1). In contrast, similarity in species composition decreases approximately exponentially with geographic distance over the entire study area (Fig. 2). The decrease is very steep over the first 4000 km, where after the similarity slowly approaches zero.

Figure 2
figure 2

Distribution of palm species richness. Similarity as a function of geographic distance between 1° × 1° grid cells. Fits are quadratic Gaussian loess fits with automatic span selection (S-PLUS 7.0). Only every 2000th data point is shown.

Within the four subregions (Table 1), both aspects of community similarity exhibited distance decay (Fig. 3 &4), but it was less regular for species richness than for species composition in the Andes, Caribbean and Central American subregions (Fig. 3 &4). At small distances, the distance decay was always strongest for similarity in species composition, as shown by the lower initial similarity values (Table 2). The same was true at larger distances, as indicated by lower quartile distances, with the exception of the Amazon subregion (Table 2; see also Fig. 3 &4). The geographically and environmentally least complex Amazon subregion (Table 1) had the highest initial similarity and greatest quartile distance for species composition indicating a low beta diversity and a low species turnover even at large distances (Table 2). The Amazon subregion also had the lowest initial similarity for species richness, but, in contrast, also the lowest quartile distance for this measure (Table 2), possibly reflecting greater regularity of the distance decay for similarity in species richness (Fig. 3 &4).

Figure 3
figure 3

Similarity of palm species richness and composition in the four subregions. Similarity as a function of geographic distance between 1° × 1° grid cells. Fits are quadratic Gaussian loess fits with automatic span selection (S-PLUS 7.0). Data points are only shown for the Amazon subregion.

Figure 4
figure 4

Similarity in species richness and composition per 1° grid cell in the four subregions. Percentage of similarity in species richness (4 maps to the left) and composition (4 maps to the right) between one single grid cell in the center of each subregion and all other grid cells within the study area. The subregions are indicated on the individual maps.

Table 1 Descriptions of the four subregions
Table 2 Initial similarity† and quartile distance†† in species richness and composition

Environmental and geographic distance as controls of community similarity

Which model that best described the variation of palm community similarity varied among community measures and areas (Table 3). Across the Americas and in the subregions, similarity in species richness depended more on environmental distance than on geographic distance, whereas similarity in species composition depends more on geographical distance than on environmental distance. This is clear from both the partial regression coefficients of the best regression models (Table 3) and from the variation partitioning (Table 4). There were two exceptions to this pattern: Geographic distance was more strongly related to richness similarity and explains more of its variation in the Amazon subregion (Tables 3, 4). Conversely, environmental distance had the strongest relationship to similarity in species composition and explained more of its variation in the Andes subregion (Tables 3, 4).

Table 3 Multiple regression analyses of species richness (r) and species composition (c)
Table 4 Partial regression analyses on species richness (R) and species composition (C)

Discussion

Applicability of Tobler's first law of geography to macroecology

Species richness and species composition constitute two fundamental aspects of community structure [38, 39]. With respect to species composition, we found a strong geographic distance decay at the bi-continental scale (Fig. 2) and though more variable, within the four smaller regions (Fig. 3 &4). Several previous studies of similarity in species composition have shown variation with geographic distance, e.g., for palms and other tropical plants at local to landscape-scales [37, 40] and large regional scales [8, 41], boreal and temperate plants at regional to continental scales [20], terrestrial and stream invertebrates at landscape-scales [42, 43], parasites on vertebrate hosts at continental scales [17, 18], and terrestrial microbial eukaryotes from local to continental scale [44] (for a recent meta analysis see [2]. Since species composition so consistently exhibits distance decay, this aspect of community structure clearly conforms to Tobler's law.

Large-scale geographic variation in species richness is one of the most studied topics in biogeography (e.g., [21, 45–48], but, in contrast to species composition, little attention has been given to the possible existence and nature of geographic distance decays in species richness. To some extent, we expect patterns of species richness and species composition to co-vary. However, since it is clearly possible for species richness to remain constant despite a complete change in species composition a tight relationship is not expected. Here, we found that similarity in species richness did not decline monotonically with geographic distance at the bi-continental scale (Fig. 2). Hence, it can be argued that geographic distance decay does not really exist for species richness at the bi-continental scale, and that, consequently, this aspect of community structure does not conform to Tobler's first law of geography. A phenomenological explanation for this result is found in the well-known latitudinal diversity gradient [49], which is also conspicuous in the American palm flora [35].

The greater applicability of Tobler's law to species composition than to species richness was further confirmed by the weaker and less regular distance decays for similarity in species richness than for species composition in three of the four subregions. A potential explanation may be that dispersal is the dominant control of similarity in species composition, while environmental conditions (in ecological and/or evolutionary time [50]) provide the main control of species richness. Distant regions can contain similar environmental conditions, e.g., on the northern and southern hemispheres. As a consequence, there need not be any distance decay for similarity in species richness. In contrast, given a single place of origin for each species and limited subsequent dispersal, a consistent distance decay for similarity in species composition is expected. Had species composition also been primarily determined by the environment, following Baas-Becking's law, patterns similar to those for richness would have been expected, i.e., generally less consistent and weaker or even absent distance decays. We note that consistent distance decays for similarity in species composition are also expected from the phenomenological perspective that species-range size frequency distributions are generally right-skewed, i.e., most species ranges are small [51].

Stronger distance decays in environmentally complex or geographically fragmented regions

Differences in distance decays of similarity may be caused by several environmental factors, taxa related characteristics such as dispersal properties of the species, spatial configuration, extent, and grain size [14, 17, 20]. These are not mutually exclusive, but likely to interact [20]. In spatially heterogeneous environments, the frequent occurrence of highly unsuitable environmental conditions (e.g., high mountain ridges) may act as barriers to dispersal and generate particularly strong distance decays in community composition. In geographically fragmented regions such as archipelagos, sea areas constitute strong barriers to dispersal for many terrestrial organisms, again resulting in strong distance decays in community composition. The hypothesis that the distance decay in community similarity would be strongest in environmentally complex or geographically fragmented regions was confirmed by our results (Table 2) supporting the view that dispersal can be limited by geographic barriers, and hence that community similarity is not alone 'selected by the environment' [10, 11].

The importance of environmental and geographic distance

The relative importance of dispersal limitation and environmental determination is a key issue in studies of species distributions and beta diversity [8]. A similar discussion is also a key focal point in studies of large-scale gradients in species richness although, in this case, the alternative to environmental control is considered to be historical factors in general [30, 52]. Time effects (time-for-speciation, time-for-immigration) are prominent among historical explanations of species richness patterns, and clearly involve dispersal limitation at the species or above-species levels [53–55]. Nevertheless, as stated in our third study hypothesis and discussed earlier, dispersal is expected to pose a stronger constraint on species composition than on species richness, while the opposite is true with respect to environmental conditions. Our results for New World palms generally provide support for this hypothesis. Hereby, additional evidence is provided for the greater importance for dispersal as a control of species composition and a greater importance of the environment as a control of species richness.

Environmental distance was always the dominant control for similarity in species richness (Tables 3, 4), except in the Amazon region. In contrast, the relative importance of geographical and environmental distance for similarity in species composition seems to depend on scale. We found geographical distance to be a stronger control of similarity in species composition at the bi-continental scale than in the smaller regions (especially in terms of variation explained, Table 4), except in the geographically fragmented Caribbean regions, where dispersal limitation would expected to especially strong. The weak role played by geographic distance in the Andes can be expected by the close juxtaposition of highly divergent environments and strong longitudinal barriers in this region. In a previous study of palm communities in a small subregion of Amazonia, the relative importance of geographic and environmental distance was also scale-dependent, with geographic distance dominating at the regional scale, while environmental distance dominated within single localities [37]. Including somewhat larger distances, a study on palm communities in the western Amazon basin reported that geographic distance was more important than environmental distance as a control of similarity in species composition [40], while environmental distance predominated in a local-scale (50 ha) study of Amazon palm species composition [56]. Similarly, Harrison et al. [57] found that in 15 taxa (including plants, vertebrates and invertebrates) beta diversity was determined by the spatial structure of the environment, and argued that the influence of distance would only be important at larger distances. Our results corroborate this idea, suggesting the distance, and by inference dispersal, becomes more important as the spatial extent increases.

Conclusion

We conclude that the applicability of Tobler's first law of geography differs among different aspects of community structure, i.e., it is strongly applicable to species composition and only partially applicable to species richness. It appears that Tobler's law is most applicable when dispersal limitation is a strong determinant of community structure and less applicable when environmental control predominates. Corroborating this interpretation, the applicability of Tobler's law to species composition appears to increase with increasing spatial extent, i.e., with increasing likelihood of dispersal limitation. As a general hypothesis, we propose that Tobler's law is highly applicable to aspects of macroecology that depend on the single place of origin of each species and the limited dispersal abilities of most macroscopic organisms. In contrast, we expect Tobler's law to be much less applicable to aspects of macroecology that are largely driven by the abiotic environment, as abiotic conditions are often similar in highly distant locations.

Methods

Study species

Distributional data was obtained by scanning all 550 palm species distribution maps from Henderson et al. [33] Field Guide to the palms of the Americas. These maps, the only data on palm distributions currently available for all of the Americas, were digitized and georeferenced in ArcView 9.0, ESRI Inc., Redlands, California, USA at a 1° × 1° grid square resolution.

Study area

Our analyses were done for the entire tropical to warm-temperate parts of the Americas (34°N – 34°S; 33°W – 120°W; 1567 grid cells) and for four subregions (700 km × 1800 km, covering 110 grid cells each) in contrasting geographic and environmental settings and placed as parallel pairs at two latitudes. Grid cells with less than 25% land cover or without palm records were excluded (Table 1). The four subregions and their geographic and environmental setting were:

1. The Amazon subregion, which has a weak north-south gradient in temperature, precipitation, and topography and has not been exposed to major tectonically events for millions of years [58]. Geographically and climatically it is the least complex among the four studied subregions.

2. The Andean subregion, which includes portions of the Ecuadorian and Peruvian cordillera and its foreland stretching into the Amazonian basin (Fig. 1, Table 1). This complex region spans a broad range of temperatures and precipitation and is geologically young, resulting from a major uplift in Late Miocene about 5 million years ago [59].

3. The Caribbean subregion, which covers the Greater Antillean archipelago formed during the Eocene 55–35 millions years ago [60]. This geographically fragmented and topographically diverse island region (Fig. 1, Table 1) located just south of the Tropic of Capricorn has a more seasonal and less humid climate than the equatorial regions.

4. The Central American subregion, which covers large parts of Mexico including most of the Yucatan peninsula, Guatemala, Belize and part of Honduras (Fig. 1, Table 1). It is climatically and topographically complex.

Environmental variables

For each grid cell nine explanatory environmentally related variables were computed: (1) mean annual temperature (°C); (2) annual precipitation (mm yr-1); (3) number of wet days per year (variables 1–3 were obtained from [67]); (4) topographical range (maximum – minimum elevation, extracted from the Digital Elevation Model from United States Geological Survey [68]; (5) number of vegetation types, computed from a vegetation map with a resolution of 1:20,000,000 [61] using the majority type option in the Zonal Statistics function in Spatial Analyst [62]; (6) soil pH; (7) percentage of sand; (8) soil cation exchange capacity; (9) percentages of CaCO3 in the soil (variables 6–9 describes 0–30 cm topsoil properties and were obtained from FAO's Digital Soil Map of the World, Version 3.5, November 1995). The variable land cover describes the percentage of land in each grid cell. The residuals from a regression between land cover and number of species per grid were used in parallel analyses. However, the influence of land cover turned out to be negligible (results not shown).

Distance matrices

All distance matrices were computed in R-package version 4.0 d6 [63]. All environmental variables were standardized and converted into Euclidean distance matrices. For species richness analyses, we used two different environment matrices, one based on all nine environmental variables (environmental distance) and one based on three climate variables (climatic distance). For species composition analyses, topographic range and number of vegetation types were excluded from the computation of environmental distance, as species composition is not expected to be related to measures of environmental heterogeneity.

Geographic distance between grid cells was calculated as the distance in kilometers between the grid cell centroids. Two geographic distance matrices were used, one based on the linear distance and one based on the ln-transformed distance. Dispersal limitation is expected to cause logarithmic distance decay according to Hubbell's neutral model [8].

Similarity in species composition was computed using the Sørensen index, while similarity in species richness was based on the Euclidean distance (D), converted to a similarity (S) using the formula S = 1 - D/Dmax, where Dmax is the maximum distance observed. Community similarity was analyzed directly or after ln-transformation [17, 19, 20].

Data analyses

To obtain an estimate of the strength of the distance decay in community similarity, we calculated initial similarity following Soininen et al. [2]. In our case, initial similarity was defined as the similarity at a distance of 150 km, to ensure that we did not calculate the similarity within just one 1° × 1° grid cell (approximately 110 km * 110 km close to the Equator) (Table 4). Furthermore, we calculated the distance at which the initial similarity was 75% of its original value (the quartile distance). This measure was inspired by Soininen et al.'s [2] halving distance, but we were not able to measure the halving distance in all subregions as the similarity sometimes did not drop below 50%. We used two different calculations depending on the form of the original regressions, linear-linear (y = α + β × x) or log-linear (y = α + β × lnx) (y = similarity at the distance x, α and β being the regression parameters) (Table 2). Initial similarity reflects turn-over of species richness or composition at relatively small spatial distances, while the quartile distance describes turn-over at broad spatial distances [2].

The importance of geographic and environmental distance as controls of community similarity was analyzed using multiple regression analyses on distance matrices [64]. Multiple regressions were run for the entire study region (the Americas) and the four subregions, separately. Four combinations of explanatory distance matrices were used: (A) environmental and linear geographic distance, (B) environmental and ln-transformed geographic distance, (C) climatic and linear geographic distance, and (D) climatic and ln-transformed geographic distance. The best model was selected as the model with the highest R2. The multiple regression analyses on distance matrices were done using Permute 3.4! with levels of significance assessed by a permutation procedure (999 permutations) that take into account the non-independence of the similarity values [64].

We partitioned the community similarity variation into its pure environmental distance (RPE), pure geographic distance (RPG), mixed geographic-environmental distance (RMX), and unexplained (RUN) fractions using partial regressions [14, 19, 36, 65, 66]. Variation partitioning was done for both measures of community similarity and for the entire study region as well as each subregion. For each data set, the best of the four models described above was used as the basis for the partitioning. Multiple regressions on both the environmental and the geographic distance matrices, the environmental distance matrices alone, and the geographic distance matrices alone were computed to obtain the total explained variation (R2 = RT), the variation explained by geographic distance (RS), and the variation explained by environmental distance (RE). Based on these values, the pure geographic distance, pure environmental distance, mixed geographic-environmental distance, and the unexplained fractions of the variation in community similarity were calculated as RPG = RT - RE; RPE = RT - RG; RMX = RT - (RE + RG) and RUN = 1 - RT [61].

References

  1. Tobler WR: A computer movie simulating urban growth in the Detroit Region. Economic Geography. 1970, 46: 234-240. 10.2307/143141.

    Article  Google Scholar 

  2. Soininen J, McDonald R, Hillebrand H: The distance decay of similarity in ecological communities. Ecography. 2007, 30: 3-12.

    Article  Google Scholar 

  3. Sui DZ: Tobler's first law of geography: A big idea for a small world?. Annals of the Association of American Geographers. 2004, 94: 269-277. 10.1111/j.1467-8306.2004.09402003.x.

    Article  Google Scholar 

  4. Pennisi E: What determines species diversity?. Science. 2005, 309: 90-10.1126/science.309.5731.90.

    Article  CAS  PubMed  Google Scholar 

  5. Legendre P: Spatial autocorrelation – trouble or new paradigm?. Ecology. 1993, 74: 1659-1673. 10.2307/1939924.

    Article  Google Scholar 

  6. Lomolino MV: Elevation gradients of species-density: historical and prospective views. Global Ecology and Biogeography. 2001, 10: 3-16. 10.1046/j.1466-822x.2001.00229.x.

    Article  Google Scholar 

  7. Bell G: Neutral Macroecology. Science. 2001, 293: 2413-2418. 10.1126/science.293.5539.2413.

    Article  CAS  PubMed  Google Scholar 

  8. Condit R, Pitman N, Leigh EG, Chave J, Terborgh J, Foster RB, Núñez VP, Aguilar S, Valencia R, Villa G: Beta-diversity in tropical forest trees. Science. 2002, 295: 666-669. 10.1126/science.1066854.

    Article  CAS  PubMed  Google Scholar 

  9. Qian H, Ricklefs RE, White PS: Beta diversity of angiosperms in temperate floras of eastern Asia and eastern North America. Ecology Letters. 2005, 8: 15-22. 10.1111/j.1461-0248.2004.00682.x.

    Article  Google Scholar 

  10. Telford RJ, Vandvik V, Birks HJB: Dispersal limitation matter for microbial morphospecies. Science. 2006, 312: 1015-10.1126/science.1125669.

    Article  CAS  PubMed  Google Scholar 

  11. Whitaker RJ, Grogan DW, Taylor JW: Geographic barriers isolate endemic populations of hyperthermophilic archaea. Science. 2003, 301: 976-978. 10.1126/science.1086909.

    Article  CAS  PubMed  Google Scholar 

  12. Guisan A, Zimmerman NE: Predictive habitat distribution models in ecology. Ecological Modelling. 2000, 135: 147-168. 10.1016/S0304-3800(00)00354-9.

    Article  Google Scholar 

  13. Pearson R, Dawson TP: Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful?. Global Ecology and Biogeography. 2003, 12: 361-371. 10.1046/j.1466-822X.2003.00042.x.

    Article  Google Scholar 

  14. Tuomisto H, Ruokolainen K, Yli-Halla M: Dispersal, environment, and floristic variation of western Amazonian forests. Science. 2003, 299: 241-244. 10.1126/science.1078037.

    Article  CAS  PubMed  Google Scholar 

  15. Holt RD, Keitt TH, Lewis MA, Maurer BA, Taper ML: Theoretical models of species' borders: single species approaches. Oikos. 2005, 108: 18-27. 10.1111/j.0030-1299.2005.13147.x.

    Article  Google Scholar 

  16. Svenning J-C, Skov F: Could the tree diversity in Europe be generated by postglacial dispersal limitation?. Ecology Letters. 2007, 10: 453-460. 10.1111/j.1461-0248.2007.01038.x.

    Article  PubMed  Google Scholar 

  17. Poulin R: The decay of similarity with geographical distance in parasite communities of vertebrate hosts. Journal of Biogeography. 2003, 30: 1609-1615. 10.1046/j.1365-2699.2003.00949.x.

    Article  Google Scholar 

  18. Krasnov BR, Shenbrot GI, Mouillot D, Khokhlova IS, Poulin R: Spatial variation in species diversity and composition of flea assemblages in small mammalian hosts: geographical distance or faunal similarity?. Journal of Biogeography. 2005, 32: 633-644. 10.1111/j.1365-2699.2004.01206.x.

    Article  Google Scholar 

  19. Qian H, Ricklefs RE: A latitudinal gradient in large-scale beta diversity for vascular plants in North America. Ecology Letters. 2007, 10: 737-744. 10.1111/j.1461-0248.2007.01066.x.

    Article  PubMed  Google Scholar 

  20. Nekola JC, White PS: The distance decay of similarity in biogeography and ecology. Journal of Biogeography. 1999, 26: 867-878. 10.1046/j.1365-2699.1999.00305.x.

    Article  Google Scholar 

  21. Hawkins BA, Field R, Cornell HV, Currie DJ, Guégan J-F, Kaufman DM, Kerr JT, Mittelbach GG, Oberdorff T, O'Brian EM: Energy, water, and broad-scale geographic patterns of species richness. Ecology. 2003, 84: 3105-3117. 10.1890/03-8006.

    Article  Google Scholar 

  22. Francis AP, Currie DJ: A globally consistent richness-climate relationship for angiosperms. American Naturalist. 2003, 161: 523-536. 10.1086/368223.

    Article  PubMed  Google Scholar 

  23. Latham RE, Ricklefs RE: Global patterns of tree species richness in moist forests: energy-diversity theory does not account for variation in species richness. Oikos. 1993, 67: 325-333. 10.2307/3545479.

    Article  Google Scholar 

  24. Gentry AH: Tree species richness of upper Amazonian forests. Proceedings of the National Academy of Science USA. 1988, 85: 156-159. 10.1073/pnas.85.1.156.

    Article  CAS  Google Scholar 

  25. Phillips OL, Hall P, Gentry AH, Sawyer SA, Vásquez R: Dynamics and species richness of tropical rain forests. Proceedings of the National Academy of Sciences USA. 1994, 91: 2805-2809. 10.1073/pnas.91.7.2805.

    Article  CAS  Google Scholar 

  26. Clinebell RR, Philips OL, Gentry AH, Stark N, Zuuring H: Prediction of neotropical tree and liana species richness from soil and climatic data. Biodiversity and Conservation. 1995, 4: 56-90. 10.1007/BF00115314.

    Article  Google Scholar 

  27. Kahn F: The distribution of palms as a function of local topography in Amazonian terra-firme forests. Experientia. 1987, 43: 251-259. 10.1007/BF01945548.

    Article  Google Scholar 

  28. Ricklefs RE, Latham RE: Global patterns of tree species richness in moist forests: distinguishing ecological influences and historical contingency. Oikos. 1999, 86: 369-373. 10.2307/3546454.

    Article  Google Scholar 

  29. Ricklefs RE: A comprehensive framework for global patterns in biodiversity. Ecology Letters. 2004, 7: 1-15. 10.1046/j.1461-0248.2003.00554.x.

    Article  Google Scholar 

  30. Turner JRG: Explaining the global biodiversity gradient: energy, area, history and natural selection. Basic and Applied Ecology. 2004, 5: 435-448. 10.1016/j.baae.2004.08.004.

    Article  Google Scholar 

  31. Uhl NW, Dransfield J: Genera Palmarum: a classification of palms. 1987, Kansas: Allen Press

    Google Scholar 

  32. Kahn F, de Granville JJ: Palms in the forest ecosystems of Amazonia. 1992, Berlin: Springer Verlag

    Book  Google Scholar 

  33. Henderson A, Galeano G, Bernal R: Field guide to the palms of the Americas. 1995, Princeton, New Jersey: Princeton University Press

    Google Scholar 

  34. Bjorholm S, Svenning J-C, Skov F, Balslev H: Environmental and spatial controls of palm (Arecaceae) species richness across the Americas. Global Ecology and Biogeography. 2005, 14: 423-429. 10.1111/j.1466-822x.2005.00167.x.

    Article  Google Scholar 

  35. Bjorholm S, Svenning J-C, Baker WJ, Skov F, Balslev H: Historical legacies in the geographical diversity patterns of New World palm (Arecaceae) subfamilies. Botanical Journal of the Linnean Society. 2006, 151: 113-125. 10.1111/j.1095-8339.2006.00527.x.

    Article  Google Scholar 

  36. Vormisto J, Svenning J-C, Hall P, Balslev H: Diversity and dominance in palm (Arecaceae) communities in terra firme forests in the western Amazon basin. Journal of Ecology. 2004, 92: 577-588. 10.1111/j.0022-0477.2004.00904.x.

    Article  Google Scholar 

  37. Normand S, Vormisto J, Svenning J-C, Grández C, Balslev H: Geographical and environmental controls of palm beta diversity in paleo-riverine terrace forrests in Amazonian Peru. Plant Ecology. 2006, 186: 161-176. 10.1007/s11258-006-9120-9.

    Article  Google Scholar 

  38. Whittaker RH: Vegetation of the Siskiyou Mountains, Oregon and California. Ecological Monographs. 1960, 30: 279-338. 10.2307/1943563.

    Article  Google Scholar 

  39. Whittaker RH: Evolution and measurement of species diversity. Taxon. 1972, 21: 213-251. 10.2307/1218190.

    Article  Google Scholar 

  40. Vormisto J, Tuomisto H, Oksanen J: Palm distribution patterns in Amazonian rainforests of NE Peru: what is the role of topographic variation?. Journal of Vegetation Science. 2004, 15: 485-494. 10.1658/1100-9233(2004)015[0485:PDPIAR]2.0.CO;2.

    Article  Google Scholar 

  41. Tuomisto H, Ruokolainen K, Aguilar M, Sarmiento A: Floristic patterns along a 43-km long transect in an Amazonian rain forest. Journal of Ecology. 2003, 91: 743-756. 10.1046/j.1365-2745.2003.00802.x.

    Article  Google Scholar 

  42. Cleary DFR, Moores AØ, Eichlorn KAO, van Tol J, de Jong R, menken SBJ: Diversity and community composition of butterflies and odonates in an ENSO-induces fire affected habitat mosaic: a case study from East Kalimantan, Indonesia. Oikos. 2004, 105: 426-446. 10.1111/j.0030-1299.2004.12219.x.

    Article  Google Scholar 

  43. Thompson RM, Townsend CR: A truce with neutral theory: local deterministic factors, species traits and neutrality determine patterns of diversity in streams. Journal of Animal Ecology. 2006, 75: 476-484. 10.1111/j.1365-2656.2006.01068.x.

    Article  PubMed  Google Scholar 

  44. Green JL, Holmes AJ, Westoby M, Oliver I, Briscoe D, Dangerfield M, Gillings M, Beattie AJ: Spatial scaling of microbial eukaryote diversity. Nature. 2004, 432: 747-750. 10.1038/nature03034.

    Article  CAS  PubMed  Google Scholar 

  45. Currie DJ: Energy and large-scale patterns of animal and plant species richness. American Naturalist. 1991, 137: 27-49. 10.1086/285144.

    Article  Google Scholar 

  46. Whittaker RJ, Willis KJ, Field R: Scale and species richness: towards a general, hierarchical theory of species diversity. Journal of Biogeography. 2001, 28: 453-470. 10.1046/j.1365-2699.2001.00563.x.

    Article  Google Scholar 

  47. Rahbek C: The role of spatial scale and the perception of large-scale species-richness patterns. Ecology Letters. 2005, 8: 224-239. 10.1111/j.1461-0248.2004.00701.x.

    Article  Google Scholar 

  48. Mittelbach G, Steiner CF, Scheiner SM, Gross KL: What is the observed relationship between species richness and productivity?. Ecology. 2003, 84: 3390-3395. 10.1890/03-3077.

    Article  Google Scholar 

  49. Hillebrand H: On the generality of the latitudinal diversity gradient. The American Naturalist. 2004, 163: 192-211. 10.1086/381004.

    Article  PubMed  Google Scholar 

  50. Svenning JC, Borchsenius F, Bjorholm S, Balslev H: High tropical net diversification drives the New World latitudinal gradient in palm species richness. Journal of Biogeography. 2008, 35: 394-406. 10.1111/j.1365-2699.2007.01841.x.

    Article  Google Scholar 

  51. Gaston KJ: The structure and dynamics of geographic ranges. 2003, Oxford: Oxford University Press

    Google Scholar 

  52. Hawkins BA, Diniz-Filho JAF, Soeller S: Water links the historical and contemporary components of the Australian bird diversity gradient. Journal of Biogeography. 2005, 32: 1035-1042. 10.1111/j.1365-2699.2004.01238.x.

    Article  Google Scholar 

  53. Stephens PR, Wiens JJ: Explaining species richness from continents to communities: the time-for-speciation effect in emydid turtles. The American Naturalist. 2003, 161: 112-128. 10.1086/345091.

    Article  PubMed  Google Scholar 

  54. Stephens PR: Historical processes enhance patterns of diversity along latitudinal gradients. Proceedings of the Royal Society of London, series B. 2006, 273: 2283-2289. 10.1098/rspb.2006.3596.

    Article  Google Scholar 

  55. Wiens JJ, Donoghue MJ: Historical biogeography, ecology and species richness. Trends in Ecology and Evolution. 2004, 19: 639-644. 10.1016/j.tree.2004.09.011.

    Article  PubMed  Google Scholar 

  56. Svenning J-C: Microhabitat specialization in a species-rich palm community in Amazonian Ecuador. Journal of Ecology. 1999, 87: 55-65. 10.1046/j.1365-2745.1999.00329.x.

    Article  Google Scholar 

  57. Harrison S, Ross SJ, Lawton JH: Beta diversity on geographic gradients in Britain. Journal of Animal Ecology. 1992, 61: 151-158. 10.2307/5518.

    Article  Google Scholar 

  58. Morley RJ: Origin and Evolution of Tropical Rain Forests. 2000, Chichester, England: John Wiley & Sons, Ltd

    Google Scholar 

  59. Morley RJ: Interplate dispersal paths for megathermal angiosperms. Perspectives in Plant Ecology, Evolution and Systematics. 2003, 6: 5-20. 10.1078/1433-8319-00039.

    Article  Google Scholar 

  60. Brown JH, Lomolino MV: Biogeography. 1998, Massachusetts: Sinuaer Associats

    Google Scholar 

  61. Bletter N, Janovec J, Brosi B, Daly DC: A digital base map for studying the Neotropical flora. Taxon. 2004, 53: 469-477. 10.2307/4135624.

    Article  Google Scholar 

  62. McCoy J, Johnston K: Using ArcGIS Spatial Analyst. 2002, California: Environmental Systems Research Institute

    Google Scholar 

  63. Casgrain P, Legendre P: The R Package for multivariate and spatial analysis, version 4.0. 2001, Department de Sciences biologiques, Universite de Montreal

    Google Scholar 

  64. Legendre P, Lapointe F-J, Casgrain P: Modelling brain evolution from behavior: A permutational regression approach. Evolution. 1994, 48: 1487-1499. 10.2307/2410243.

    Article  Google Scholar 

  65. Legendre P, Legendre P: Numerical Ecology. 1998, Amsterdam, the Netherlands: Elsevier

    Google Scholar 

  66. Duivenvoorden JF, Svenning J-C, Wright SJ: Beta diversity in tropical forests. Science. 2002, 295: 636-637. 10.1126/science.295.5555.636.

    Article  CAS  PubMed  Google Scholar 

  67. Intergovernmental Panel on Climatic Change. [http://www.ipcc.ch]

  68. U.S. Geological Survey. [http://edc.usgs.gov]

Download references

Acknowledgements

The Faculty of Natural Sciences at the University of Aarhus supports Stine Bjorholm's PhD study. The Danish Natural Science Research Council for economic support (grant 21-04-0346 to Jens-Christian Svenning and grant 272-06-0476 to Henrik Balslev).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Henrik Balslev.

Additional information

Authors' contributions

SB compiled the GIS-data with help from FS. Data analysis was performed by SB and J-CS. The manuscript was written by SB, J-CS and HB. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Bjorholm, S., Svenning, JC., Skov, F. et al. To what extent does Tobler's 1st law of geography apply to macroecology? A case study using American palms (Arecaceae). BMC Ecol 8, 11 (2008). https://doi.org/10.1186/1472-6785-8-11

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1472-6785-8-11

Keywords