Open access

Species Distribution Patterns, Species-Area and Species- Temperature Relationships in Eastern Asian Plants

Written By

Jianming Deng and Qiang Zhang

Submitted: 26 April 2012 Published: 29 August 2012

DOI: 10.5772/48509

From the Edited Volume

Biodiversity Conservation and Utilization in a Diverse World

Edited by Gbolagade Akeem Lameed

Chapter metrics overview

2,771 Chapter Downloads

View Full Metrics

1. Introduction

It is a fascinating issue for ecologists to develop a general theory or principle to interpret the mechanisms of global gradients and stabilization of biodiversity. This question has perplexed biogeographers and ecologists for about 100 years, and the diverse theories and hypotheses have been put forward to account for latitudinal gradients in biodiversity (Wright 1983; Rohde, 1992; Waide et al., 1999; Colwell & Lees, 2000; Gaston, 2000; Allen et al. 2002; Hawkins et al., 2003; Willig et al. 2003; Ricklefs, 2004; Evans & Gaston, 2005; Evans et al., 2005; Mittelbach et al. 2007; Gillooly& Allen, 2007; Storch et al. 2007; Cardinale, et al., 2009), Recent decade, the metabolic theory of biodiversity (MTB) is developed and attracting a lot of attentions of ecologists as a novel hypothesis based on metabolic theory of ecology (MTE) and the energetic-equivalence rule (West et al. 1997, 1999; Enquist et al. 1998; Allen et al. 2002, 2007; Brown et al. 2004; Deng et al. 2006, 2008). The MTB is recognized as a general principle that can quantify relationships between the dynamic processes of population and biodiversity patterns in ecosystem, and between species richness and environmental factors (see also Allen et al. 2003, 2006; Allen & Gillooly 2006; Gillooly & Allen 2007). The metabolic eco-evolutionary model of biodiversity, the most recent extension of the MTB, has been developed by Stegen et al., (2009).

The MTB considered that species richness, S, in plots of fixed area, A, should be described by a form of equation as followingS=(J/A)(B0/BT¯)eE/kT. In this expression, E is the activation energy of metabolism, -0.6 to -0.7 eV, k is Boltzmann’s constant (8.62×10-5 eV K-1, where K is degrees Kelvin), T is environmental temperature in degrees kelvin,B0=b0M3/4 , where b0 is a normalization constant that varies by taxonomic group, M is individual body size, BT is the total energy flux of a population per unit area, varying with taxonomic group and plot area A, and J/A is the total density of individual per unit area. Apparently, the species richness should vary as a function of abundance, body size and environmental temperature. So, when abundance and size are both presumed constants across geographical space, the relationship between the natural logarithm of species richness (lnS) and the inverse temperature (1/kT) will be re-expressed by a linear equation:lnS=E(kT)1+C. In this equation, the intercept term (C=ln[(B0/BT¯)(J/A)]) incorporates the effect of mean body size of the study taxon, area and total community abundance (Allen et al., 2002; Gillooly & Allen, 2007).

The intense and continuous controversies for the MTB have been focusing on two primary predictions: 1) whether ln-transformed species richness is linearly associated with an inverse rescaling of ambient temperature or not, and 2) if so, whether the slope of the relationship is encompassed in the theoretical value range of -0.6 to -0.7. The proponents argue that this theory accounts for diversity gradients over a range of spatial scales from mountain slopes to continental and global gradients, and for many groups of plants and ectothermic animals (Allen et al. 2002, 2003;Brown et al., 2003, 2004; Gillooly & Allen, 2007). Kaspari et al.’s (2004) and Hunt et al.’s (2005) analyses using respectively ant communities and deep-sea communities datasets were also in favor of the Allen et al.’s (2002) predictions. Concomitantly, the disagreements for MTB emerged in some literature. Hawkins et al. (2007a) tested the predictions of this theory with 46 different data sets compiled from a variety of terrestrial plants, invertebrates, and ectothermic vertebrates, and found that the results were partly deviated from the predictions of the MTB (Allen et al. 2002; Brown et al. 2004), Accordingly, they considered that MTB was a poor predictor for the observed diversity gradients in most terrestrial system. Latimer (2007) subsequently reanalyzed some Hawkins et al.’s (2007a) data sets using a Bayesian approach and supported their conclusions. Algar et al. (2007) recently showed that the relationship between richness and temperature was actually curvilinear for several data sets in North America, and slopes varied systematically in geographical space, which were strongly consistent with Cassemiro et al (2007) analyzing for New World amphibians. As a consequence, they claimed that Allen et al.’s (2002) model did not give an adequate fit to the data.

Figure 1.

Location of the eleven floristic regions and 270 nature reserves of China used in this study. Floristic regions was marked with Arabic numerals as following (1, Daxinganling; 2, NE China; 3, NE China plain; 4, North China region; 5, East China; 6, Lingnan region China; 7, Tsinling Mountains; 8, The region of Hengduan Mountain; 9, Central China; 10, Dian-qian-gui region; 11, The region of Yunnan Plateau)(also see Zhang et al. 2011).

Sample typeVariableNMinimmMaximumMeanSDSkewnessKurtosis
Floristic regionSpecies richness
Seed plant11101979544392.452524.83-0.03-1.47
Climatic variables
MAT (°C)11-2.8020.8812.347.48-0.980.002
Location and area
Mean latitude (°)11235035.339.540.35-1.49
Mean longitude (°)11100132115.19.940.25-0.44
Area (km2)11520009600004016602730700.640.73
Nature reserveSpecies richness
Vascular plant19313845431388.18809.920.860.99
Climatic variables
MAT (°C)270-2.82913.935.71-0.610.44
Location and area
Mean latitude (°)27018.451.630.746.550.710.37
Mean longitude (°)27095130.6111.887.040.12-0.37
Area (km2)2700.646698490.89922.283.9418.69

Table 1.

Summary statistics of species richness in different plant groups, climate variables, and areas of 11 floristic regions and 270 nature reserve used in this paper.

Empirical evaluations of how well observed richness patterns fit the central predictions of the MTB are now appearing in several literature (Allen et al., 2002; Kaspari et al., 2004; Hunt et al., 2005; Algar et al., 2007; Cassemiro et al., 2007; Latimer, 2007; Hawkins et al., 2007a, b; Sanders et al., 2007). Although to date the observed patterns in biodiversity have been taxonomically and geographically limited (Ellison, 2007), the data sets for the detailed plant groups are relatively absent. Wang et al. (2009) showed that magnitude of temperature dependence (i.e. E) of tree species richness in both eastern Asia and North America increases with spatial scale at the large scale of 50×50km to 400×400km. Therefore, we conjectured that the species richness inherent dependence of spatial scale may influence on the successful tests for the predictions of the MTB (also see Zhang et al. 2011). However, it is unclear (i) how the species richness responds to temperature at the variant spatial scale, especially at the small scale level. (ii) what spatial scale range is appropriate to the MTB.

Here we aimed to evaluate how the relationship between species richness and temperature predicted by MTB varied with respect to sampling scales, as well as with respect to different plant taxonomic group using an extensive plant data sets including three divisions in vascular plant at two different sample scales including nature reserve grain and floristic grain.


2. Methods

We compiled species richness and other basic characteristics of 11 floristic regions and 270 natural reserves. All of the plant species richness data sets used in our analysis were collected from the previous reports involving eleven floristic regions and 270 nature reserves across the eastern China (Zhao & Fang, 2006, many others; for details see Zhang et al.2011). All species were compiled and classified into three groups (pteridophyte, gymnosperm, and angiosperm) at both floristic and reserve scales (the details see Zhang et al.2011). Here the alien species were excluded from our data analyses and only the native species retained. The areas of nature reserves and floristic regions were respectively range from 0.64 to 6689 and from 52000 to 960000 square kilometers (km2) between 18.4° N and 51.6° N latitude and between 95° E and 130.6° E longitude covering a total terrestrial area of 132,540 km2 (See Fig1 and Table1). The temperature and the size distribution of the 270 nature reserves also were showed in Table 1.

The mean annual temperature (MAT), assigned to each nature reserve based on its location and used to analyze the relationship between temperature and species richness, was compiled from a 1971-2000 temperature database of China generated from 722 climate stations across China. Flora’s MAT was an average value of all the covered climate stations within each floristic region. Other environmental variables such as geographical range and area were also documented.

Descriptive statistics of plant species richness and environmental variables were produced to interpret the information on the data distributions (Table1). The observed slopes of In-transformed richness versus 1/kT relationships and at the floristic and reserve grains across the taxonomic group (three divisions and five families of the angiosperms) were estimated by reduced major axis (RMA) regression (Brown et al., 2004; Hawkins et al., 2007). The species richness-area relationships were also analyzed to evaluate the effect of the region size on the species richness (See Fig 4 and 5).


3. Results

The natural logarithm of species richness was significantly linear with (kT)-1 at the floristic grain. The 95% CI values of all slopes estimated by RMA did not exclude the second primary prediction of MTB, implying that no significant heterogeneity occurred in slopes among taxonomic groups (Fig2; Table2). However, at the nature reserve grain, only two of three large plant groups (angiosperm, Fig3-f and pteridophyte, Fig3-h)showed significantly linear relationship and met the first prediction of MTB, and gymnosperm were not linear and rejected the entire MTB (Fig3-g; Table3). The slope values estimated by RMA regression for all taxonomic groups were significantly exclusive from the second prediction of MTB (Table2).

The species-area relationships for all taxonomic divisions at both the floristic and nature reserve special scales indicated that the area size of community have more impact on the species richness for subdivision (e.g. family) than for division (Fig. 4 and 5). Moreover, the observed slope values were close to or encompass (95% CI) the theoretical values predicted by MBT at the spatial scale range of 50- 6698 km2, excluding the size of area class less than 50 km2 (Fig. 6; Table 4).

Figure 2.

The relationship between natural logarithm of species richness (lnS) and inverse temperature (1/kT) for seven groups in 11 floristic regions: two divisions (gymnosperm and angiosperm) and five families of angiosperm (Compositae, Poaceae, Rosaceae, Liliaceae and Labiatae).

GroupFigureNR2PRMA slope(95%CI)
Compositae2-a90.640.009-0.48(-0.74– -0.23)
Poaceae2-b90.82<0.001-0.55(-0.76– -0.34)
Rosaceae2-c90.730.003-0.66(-0.97– -0.35)
Liliaceae2-d90.670.007-0.46(-0.69– -0.22)
Labiatae2-e90.89<0.001-0.70(-0.92– -0.50)
Angiosperm2-f110.86<0.001-0.64(-0.83– -0.46)
Gymnosperm2-g100.710.002-0.80(-1.15– -0.45)
Seed plant2-h110.89<0.001-0.64(-0.80– -0.48)

Table 2.

Summary of regressions testing Model II (RMA) slopes of richness-temperature relationships for cases with linear relationship between inverse scaled temperature and ln-transformed richness in 11 floristic regions.

Figure 3.

The relationship between natural logarithm of species richness (lnS) and inverse temperature (1/kT) for eight plant groups in 270 nature reserve: three divisions (Pteidophyte, Gymnosperm and Angiosperm) and five families of angiosperm (Compositae, Poaceae, Rosaceae, Liliaceae and Labiatae).

GroupFigureNR2PRMA slope(95%CI)
Compositae3-a710.050.080.73(0.56– 0.91)
Poaceae3-b700.000.58-0.71(-0.89– -0.54)
Rosaceae3-c710.050.050.96(0.73– 1.18)
Liliaceae3-d490.020.360.81(0.58– 1.05)
Labiatae3-e560.010.60-0.81(-1.03– -0.59)
Angiosperm3-f2550.050.09-0.90(-1.01– -0.79)
Gymnosperm3-g2340.000.326-0.91(-1.03– -0.79)
Pteridophyte3-h1890.21<0.001-1.55(-1.75– -1.35)
Vascular plant3-i1930.09<0.001-0.82(-0.93– -0.71)

Table 3.

Summary of regressions testing Model II (RMA) slopes of richness-temperature relationships for cases with linear relationship between inverse rescaled temperature and ln-transformed richness in nature reserve.

Figure 4.

Species richness-area relationships for the seven plant groups in floristic regions: two divisions (gymnosperm and angiosperm) and five families of angiosperm (Compositae, Poaceae, Rosaceae, Liliaceae and Labiatae).

Figure 5.

Species richness-area relationships for the eight plant groups in 270 nature reserves: three divisions (Pteidophyte, Gymnosperm and Angiosperm) and five families of angiosperm (Compositae, Poaceae, Rosaceae, Liliaceae and Labiatae).

Figure 6.

Temperature–richness relationships for the vascular plant group along an area classes. Note that the classification of area is following: a)<50, b)50-100, c)100-200, d)200-400, e)400-800, f)800-1600, g)>1600.

Area classes(km2)FigureNR2PRMA slope(95%CI)
50-1007-b250.040.32-0.72(-1.02– -0.41)
100-2007-c440.25<0.001-0.63(-0.80– -0.46)
200-4007-d330.060.19-1.24(-1.69– -0.80)
400-8007-e220.360.003-0.60(-0.82– -0.37)
800-16007-f120.250.08-0.96(-1.55– -0.38)
"/>16007-g160.180.10-0.84(-1.27– -0.40)

Table 4.

Summary of regressions testing Model II (RMA) slopes of richness-temperature relationships for vascular plant group with linear relationship between inverse rescaled temperature and ln-transformed richness along an area classes.


4. Discussion

Hawkins’ et al. (2007a) suggested that the relationship of logarithm transformed species richness and inverse temperature was nonlinear through analyzing the datasets of Chinese angiosperm taken from nature reserves with a range of area from 100 km2 to 247 km2. Here we similarly failed to observe significantly linear relationships between them at the nature reserve grain with the regions ranging from 0.64 km2 to 6689 km2, excepting for two large groups (angiosperm and pteridophyte). Moreover, almost all slope values were exclusive from the predictive range of MTB (Table 2) as the pattern of tree species distribution in eastern Asia (Wang et al., 2009). However, when we analyzed these data sets at the floristic regions ranging from 52000 km2 to 960000 km2, not only this linear relationship was observed, but also the slopes is highly in agreement with the theoretical values of MTB (Allen et al. 2002; Brown et al. 2004). Therefore, the plant species richness patterns predicted by MTB apparently depended on the grain size (Ellison, 2007). This scenario may be due to the fact that the number of species at the large scale overwhelmed the number of species at the relative small sample scale (e.g. nature reserve). However our analysis of species richness-area relationships showed no significant relations at floristic grain (Fig 4). The adjacent nature reserves frequently have the similar annual temperature, but the other environmental factors (i.e. water, elevation and nutrition) may exhibit a lot of variations between them that can also strongly influence the local plant species richness (Storch et al., 2007). The large-scale (floristic region) patterns are not simply explicable in terms of knowledge of small-scale (nature reserve) processes (Storch and Gaston, 2004). On the contrary, despite the habitat heterogeneity including annual temperature is large between plant flora, it is usually overwhelmed within plant flora because of the enormous sample scale (Field et al., 2009).

For the purpose of evaluating the MBT’s robustness, Hawkins et al. (2007a) show the relationship between the inverse of temperature and the natural log of richness in terrestrial ectotherms (including amphibians, reptiles), invertebrates, mammals and plant around the world. However, in their plant data sets, detailed taxonomic unit (e g, pteridophyte, gymnosperm and family unit) were not contained. In their 46 data sets, 14 had no significant relationship; 9 of the remaining 32 were linear, meeting the first prediction of the MBT, but the slope values against its second prediction. So, they contended that it was important to use appropriate taxonomic ranges for accepting or refusing the prediction of MBT (see also, Ellison, 2007).

Our results clearly showed that the significant taxonomic dependence in the nature reserve data sets. Pteridophyte unit which potentially supported the first prediction of MTB dominantly differs from the other groups in particular. Pteridophytes have a reproductive strategy based on the high dispersibility of spores, and have a strong moisture dependence of the sexual reproduction (Pausas & Sáez, 2000; Lehmann, et al., 2002; Castán & Vetaas, 2005). Thus, the life history and growth cycle for pteridophytes are probably more directly and tightly linked to abiotic factors than many other groups of plants because of the lack of co-evolved relationships with animal vectors (Barrington, 1993; Lwanga et al., 1998; Pausas & Sáez, 2000; Castán & Vetaas, 2005). So our results possibly support the perspective that the ability of MBT to predict richness patterns will also depend on dispersal ability (Latimer, 2007).

The plant species were not subdivided into division group to test the slopes converge around the predicted value -0.65 by MTB (Allen et al. 2002; Brown et al. 2004). Whereas the significant heterogeneity of slopes were observed at both floristic region and reserve scale among the different taxonomic groups as the most recently reported by Hawkins et al. (2007a,b) and Wang et al.(2009), indicating that the plant groups may hold variable activation energies rather than an invariant value. Our more recently research showed that validity of the MTB lies on if the area size of the community has no significant effect on species richness (Zhang et al. 2011). Therefore we believe that the slope value for each taxonomic group should be co-influenced by the restriction of distribution range, the area size of sampling community and other abiotic factors, as well as the inherent activation energy differences.


5. Conclusion

Our results suggested that the relationship predicted by MTB between the plant richness and temperature can be tested at the larger regional scale (e.g. floristic region) well. However, at the small scale (e.g. nature reserve), the predicted relationships were easily influenced by the many other factors such as area size of community, taxonomic divisions, seed dispersal and so on. Allen et al. (2003) claimed that the theory of biodiversity proposed by themselves is not complete and comprehensive. Here we consider that the theory must integrate the fundamental influences of multifactor involving temperature, area size, water, elevation and nutrition on the species richness patterns in small scale regions where the disturbance of environmental factors easily result in change of the species diversity. At the same time, we should also seek the more biological interpretation for the noticeable differences among taxonomic groups in the future.


This study were supported by the Natural Science Foundation of China (31000286), the Program for New Century Excellent Talents in University to J.M.D. and Key Project of Ministry of Education of China (no. 109152).


  1. 1. AlgarA. C.KerrJ. T.CurrieD. J.2007A test of metabolic theory as the mechanism underlying broad-scale species-richness gradients. Global Ecol. Biogeogr 16170178
  2. 2. AllenA. P.BrownJ. H.GilloolyJ. F.2002Global biodiversity, biochemical kinetics, and the energetic-equivalence rule. Science 29715451548
  3. 3. AllenA. P.GilloolyJ. F.2006Assessing latitudinal gradients in speciation rates and biodiversity at the global scale. Ecol Lett 9947954
  4. 4. AllenA. P.GilloolyJ. F.BrownJ. H.2003Response to comment on “Global biodiversity, biochemical kinetics and the energetic-equivalence rule”. Science 299, 346c.
  5. 5. AllenA. P.GilloolyJ. F.SavageV. M.BrownJ. H.2006Kinetic effects of temperature on rates of genetic divergence and speciationProc. Natl Acad. Sci. USA 10391309135
  6. 6. BrownJ. H.GilloolyJ. F.AllenA. P.SavageV. M.WestG. B.2004Toward a metabolic theory of ecology. Ecology 8517711789
  7. 7. CardinaleB. J.HillebrandH.HarpoleW. S.GrossK.PtacnikR.2009Separating the influence of resource ‘availability’ from resource ‘imbalance’ on productivity-diversity relationshipsEcology letters12475487
  8. 8. CassemiroF. A. S.BarretoB. S.RangelT. F. L. V. B.Diniz-FilhoJ. A. F.2007Non-stationarity, diversity gradients and the metabolic theory of ecology. Global Ecol. Biogeogr 16820822
  9. 9. ColwellR. K.LeesD. C.2000The mid-domain effect: geometric constraints on the geography of species richnessTrends in Ecology and Evolution, 157076
  10. 10. DengJ. M.LiT.WangG. X.LiuJ.ZhaoC. M.JiM. F.ZhangQ.LiuJ. Q.2008Trade-offs between the metabolic rate and population density of plants. Plos One, 3 (3), 1799.
  11. 11. DengJ. M.WangG. X.MorrisE. C.WeiX. P.LiD. X.ChenB. M.ZhaoC. M.LiuJ.WangY.2006Plant mass-density relationship along a moisture gradient in north-west China. Journal of Ecology 94953958
  12. 12. EllisonA. M.2007Metabolic theory and patterns of species richness. Ecology1889 EOF
  13. 13. EnquistB. J.BrownJ. H.WestG. B.1998Allometric scaling of plant energetics and population densityNature395163165
  14. 14. EvansK. L.GastonK. J.2005Can the evolutionary-rates hypothesis explain species-energy relationships? Functional Ecology, 19899915
  15. 15. EvansK. L.WarrenP. H.GastonK. J.2005Species-energy relationships at the macroecological scale: a review of the mechanismsBiological Reviews80125
  16. 16. FieldR.HawkinsA. B.CornellH. V.CurrieD. J.Diniz-FilhoJ. A. F.GuéganJ. F.KaufmanD. M.KerrJ. T.MittelbachG. C.OberdorffT.O’BrienE. M.TurnerJ. R. G.2009Spatial species-richness gradients across scales: a meta-analysis. Journal of Biogeography, 36132147
  17. 17. GastonK. J.2000Global patterns in biodiversity. Nature, 405220227
  18. 18. GilloolyJ. F.AllenA. P.2007Linking global patterns in biodiversity to evolutionary dynamics using metabolic theory.Ecology8818901894
  19. 19. HawkinsB. A.AlbuquerqueF. S.AraújoM. B.BeckJ.BiniL. M.Cabrero-SańudoF. J.Castro-PargaI.Diniz-FilhoJ. A. F.Ferrer-Castá al.2007aA global evaluation of metabolic theory as an explanation for terrestrial species richness gradients. Ecology 8818771888
  20. 20. HawkinsB. A.Diniz-FilhoJ. A. F.BiniL. M.AraújoM. B.FieldR.Horta,lJ.KerrJ. T.RahbekC.RodríguezM. Á.SandersN. J.2007bMetabolic theory and diversity gradients: where do we go from here? Ecology, 8818981902
  21. 21. HawkinsB. A.FieldR.CornellH. V.CurrieD. J.GueganJ. F.KaufmanD. M.KerrJ. T.MittelbachG. G.OberdorffT.O’BrienE. M.PorterE. E.TurnerJ. R. G.2003Energy, water, and broad-scale geographic patterns of species richnessEcology8431053117
  22. 22. HuntG.CroninT. M.RoyK.2005Species-energy relationship in the deep see: a test using the Quaternary fossil record. Ecol Lett 8739747
  23. 23. KaspariM.WardP. S.YuanM.2004Engery gradients and the geographic distribution of local ant diversity. Oecologia 140407413
  24. 24. LatimerA. M.2007Geography and resource limitation complicate metabolism-based predictions of species richness.Ecology8818951898
  25. 25. MittelbachG. G.SchemskeD. W.CornellH. V.AllenA. P.BrownJ. M.BushM. B.HarrisonS. P.HurlbertA. H.KnowltonN.LessioH. al.2007Evolution and the latitudinal diversity gradient: speciation, extinction and biogeography. Ecol Lett, 10315331
  26. 26. RicklefsR. E.2004A comprehensive framework for global patterns in biodiversityEcol Lett 7115
  27. 27. RohdeK.1992Latitudinal gradients in species diversity: the search for the primary causeOikos65514527
  28. 28. RosenzweigM. L.1995Species Diversity in Space and TimeCambridge University Press, Cambridge, UK.
  29. 29. SandersN. S.LessardJ. P.FitzpatrickM. C.DunnR.2007Temperature, but not productivity or geometry, predicts elevational diversity gradients in ants across spatial grains. Global Ecol. Biogeogr 16640649
  30. 30. 2009 RStegenJ. C.EnquistB. J.Ferriere20.R.Advancing the matabolic theory of biodiversity. Ecology Letters 1210011015
  31. 31. StorchD.2003Comment on “Global biodiversity, biochemical kinetics, and the energetic-equivalence rule”. Science 299, 346.
  32. 32. StorchD.GastonK. J.2004Untangling ecological complexity on different scales of space and timeBasic and Applied Ecology5389400
  33. 33. StorchD.MarquetA.BrownJ. H.2007Scaling BiodiversityCambridge University Press, Cambridge, UK.
  34. 34. WaideR. B.WilligM. R.SteinerC. F.MittelbachG.GoughL.DodsonS. I.JudayG. P.ParmenterR.1999The relationship between productivity and species richness. Annual Review of Ecology and Systematics, 30257300
  35. 35. WangZ.BrownJ. H.TangZ.FangJ.2009Temperature dependence, apatial scale, and tree species diversity in eastern Asia and North America. Proc. Natl Acad. Sci. USA 1061338813392
  36. 36. WestG. B.BrownJ. H.EnquistB. J.1997A general model for the origin of allometric scaling laws in biology.Science 276122126
  37. 37. WestG. B.BrownJ. H.EnquistB. J.1999A general model for the structure, function, and allometry of plant vascular systems. Nature 400664667
  38. 38. WilligM. R.KaufmanD. M.StevensR. D.2003Latitudinal gradients of biodiversity: pattern, process, scale, and synthesisAnnu Rev Ecol Evol Syst 34273309Wright, D. H. 1983 Species-energy theory: an extension of species-area theory. Oikos 41, 496-506.
  39. 39. ZhangQ.WangZ. Q.JiM. F.FanZ. X.DengJ. M.2011Patterns of species richness in relation to temperature, taxonomy and spatial scale in eastern ChinaActa Oecologica307 EOF313 EOF
  40. 40. ZhaoS.FangJ.2006Patterns of species richness for vascular plants in China’s nature reservesDiversity Distrib 12364372

Written By

Jianming Deng and Qiang Zhang

Submitted: 26 April 2012 Published: 29 August 2012