Relationships between biodiversity and ecosystem functioning are vital to understand ecosystem properties, and have frequently been addressed in small-scale studies. However, interactions and changes differ at large scales, and should be similarly evaluated to monitor biodiversity and ecosystem functional alterations. In the present study, Mainland China was divided into 241 quadrats of 2° latitude by 2° longitude. Ecosystem function was comprehensively assessed using three indicator variables in each quadrat, primary productivity, bird species richness, and relative humidity. Relationships between each ecosystem function variable were regressed against seed plant species richness. All three indicators exhibited the same change model, a linear model when plant species richness was under 5.000 species, and a hump-back/quadratic model when seed plant richness was over 5.000 species, with an increase in seed plant species richness at a larger scale.
Ecosystem function is the capacity for natural processes and components integral to those processes to provide goods and services, which ultimately meet human needs, either directly or indirectly (de Groot et al. 2002). Indirectly, these services provide refuge and reproductive habitat to wild plants and animals. Many studies have demonstrated the importance of biodiversity to ecosystem functioning, and primary productivity has been representative of ecosystem function (Thompson et al. 2005, Grace et al. 2007). Primary productivity is the rate at which plants and other photosynthetic organisms (e.g., phytoplankton, macroalgae) produce organic compounds i.e., biomass in an ecosystem.
A persistent debate continues regarding the relevance of studies that underlie the conclusions explaining biodiversity and functioning in mature natural ecosystems. Most experimental studies of constructed species assemblages determined that increasing diversity contributes to greater biomass production (e.g., Hooper et al. 2005, Spehn et al. 2005). Alternatively, species reduction experiments in natural systems revealed biomass was resilient to decreased biodiversity (Smith et al. 2004, Suding et al. 2006), while Foster et al. (2004) demonstrated during species addition experiments that species pools or communities were integral to increased biomass. These experiments were comparatively small. For large-scale studies, the productivity gradient was very often associated with plant diversity. In most cases, a unimodal (also called “hump-back” or “humped-backed”) relationship emerged, which exhibited the highest productivity at intermediate diversity levels (Grace 1999, Keddy 2005). This approach rapidly gained enthusiastic popularity among ecologists, and was soon widely recognized as a general empirical relationship (e.g., Begon et al. 1996). Although later studies reported positive or negative response patterns, the unimodal relationship remains the most commonly reported in plants (Hector et al. 1999, Waide et al. 1999, Mittelbach et al. 2001, but see Gillman & Wright 2006).
In the present study, we explored the relationship between plant ecosystem biodiversity, and three potential ecosystem function indicators (primary productivity, bird species richness, and annual mean relative humidity) at a large scale in Mainland China.
Materials and methods
Mainland China was divided into 241 quadrats (including Hainan Island as one quadrat), each 2° latitude by 2° longitude (approximately 36 450 km2). If a quadrat had an area less than half a standard quadrat on national boundaries, it would be merged into its neighboring quadrat (Zhang 2002). A small geographic region such as a 2° latitude by 2° longitude quadrat would typically be characterized as a group of ecosystems (Zhang 2002, Zhang et al. 2006). However, ecosystems in the quadrats were similar in composition, structure, and inter- and intra-environments, and correspondingly similar in function (Zhang & Dong 2010). Therefore, for the purposes of this study, we assumed the relationship between biodiversity and ecosystem functions could be evaluated at the quadrat level. Seed plant and bird species occurring in each quadrat were counted based on available records (e.g., national flora and fauna, provincial flora, and local research literature, among other appropriate reference materials - Wu 1980, Wu 1987, Wu 2004, Chen 1993, Editor Committee of Fauna Sinica 1998, Northwest Research Institute of Botany 2000, Zheng 2005). Primary productivity of each quadrat was estimated by applying the following relationship (eqn. 1):
where L = 300 + 25T + 0.05T2, N is the annual precipitation, L is the annual maximum evaporation, and T the annual mean temperature in each quadrat (see Lieth & Box 1972, Li et al. 1998). This equation was the best fit with least error in simulating primary productivity based on observed and remotely sensed data in China (Li et al. 1998). For each quadrat, annual mean relative humidity, and climatic variables were obtained from 30-year records representing 830 county meteorological stations. Average values were used if more than one meteorological station was located in a quadrat; if no meteorological stations were present in a quadrat, data from the nearest station in its neighbor quadrat was selected.
Regression analysis was chosen to assess the relationships between plant species richness and ecosystem function indicators, and analyzed using the SPSS statistical package (SPSS Inc., Chicago, Illinois, USA). AIC (Akaike’s Information Criteria) was applied to compare different regression models as follows: AIC = 2k - 2 ln(L), where k is the parameter number, and L is the maximized value of the likelihood function in the regression model (Akaike 1974).
Results
Linear and quadratic models were compared in an analysis of the following relationships: ecosystem primary productivity, bird species richness, and relative humidity, each against plant species richness. Results indicated quadratic models were more robust (with smaller AIC) when all quadrat data were used, but linear models were better (with smaller AIC) when the three highest quadrats were excluded from the analysis (Tab. 1). A linear model could be chosen when plant species richness is less than 5000 species, and the quadratic model could be applied when plant species richness exceeds 5000 species (Tab. 1). Ecosystem primary productivity was significantly associated with plant species richness; and a hump-back model relationship between primary productivity and seed plant species richness was evident (R2 = 0.6671, P < 0.001, n = 241 - Fig. 1). Primary productivity exhibited steep increases within the first 2000 seed plants, and continued to increase until reaching a maximum biomass of 2100 g m-2 at 3500 species (Fig. 1). Productivity decreased when the number of seed plants was over 4000 species. Bird species richness was also significantly correlated with plant species richness in mainland China (R2 = 0.2121, P < 0.001, n = 241 - Fig. 2). The number of bird species increased rapidly as seed plant species richness reached 1000 species; and subsequent slow increases occurred until reaching the maximum number of plant species (i.e., 5000). The relationship between bird and plant species richness also followed a loosely defined hump-back model.
Results depicted a close relationship between relative humidity and the number of seed plant species in mainland China (R2 = 0.4727, P < 0.001, n = 241 - Fig. 3). The number of seed plant species increased with a rise in relative humidity, which reached a plateau at 80% relative humidity within 2000 species, and reached a maximum at 3500 seed plant species. Gradually, seed plant species richness over 4000 species was associated with lower relative humidity (Fig. 3). Therefore, the relationship between ecosystem relative humidity, and seed plant species richness was quadratic (hump-back).
Discussion
Primary productivity is the fundamental indicator of ecosystem function, and is at the base of ecosystem function (Díaz & Cabido 2001, Hughes & Petchey 2001, Partel et al. 2007, Zhang & Zhang 2007). Animal species richness is typically associated with food plant availability, and the diversity of vegetation habitats in an ecosystem (Loreau 2001, Grime 1997). Therefore, animal species richness can be used as an indicator of ecosystem functioning, and bird species richness is particularly suitable to this purpose, as bird species diversity is dependent on species composition in a given ecosystem, horizontal and vertical habitat structure, patch mosaic structure, and ecosystem area (Zheng 2005, Hu et al. 2001). Relative humidity is a climatic index, which is related to some but not all climatic factors, including precipitation, wind, and temperature. However, humidity is also related to the ecosystem function, as it is related to ecosystem type, structure, and distribution area (Huston 1993, Tilman et al. 2001, Zhang et al. 2006). In low and medium precipitation regions, the humidity on sunny days depends largely on evaporation and transpiration (Li et al. 1998, Hooper et al. 2005), therefore it is suitable to consider it as a proxy for ecosystem functioning. All three ecosystem function indicators on a large scale generated the same model; a linear model when plant species richness was under 5000 taxa, and quadratic model when plant species richness was over 5000. The regression analysis comparisons of the relationships between plant species richness and relative humidity, bird species richness, and primary productivity using linear and quadratic models with and without the three highest quadrats demonstrated that both models produced reasonable results (Tab. 1). This may suggest a rule for the relationship between biodiversity and ecosystem functioning at a large scale (Gillman & Wright 2006).
The linear and hump-back relationships have been supported by several empirical studies, but showed variability at the small scale. Indeed, positive and negative curve patterns of biodiversity and ecosystem functions have previously been reported (Hector et al. 1999, Waide et al. 1999, Gillman & Wright 2006). This may be due to the problems of scale (Mittelbach et al. 2001, Chase & Leibold 2002), and the species richness gradient length. Experiments addressing a few species, monocultures, and mixtures are not comparable to the complexity of a natural ecosystem (Hillebrand & Matthiessen 2009). Tilman et al. (2001) suggested the biological features of selected experimental species might exhibit one or more important influences on results. Another factor affecting results is species removal from study under natural conditions, mainly where the species is limited in number. In addition, grasses or herbaceous species are frequently used in these experiments, but woody species are excluded due to difficulties in cultivating woody taxa under greenhouse conditions. Most experiments at the small scale are carried out within one- or two-year time frame, and therefore historical and evolutionary effects are not evaluated. Large-scale biodiversity and ecosystem functioning studies are vital, and cannot be substituted by small scale experimental research (Chapin et al. 2000, Chase & Leibold 2002).
Moreover, productivity is one indicator frequently assessed for ecosystem function in empirical studies. It is a viable parameter because productivity is the primary result of ecosystem function (Grace et al. 2007, Partel et al. 2007). However, it emphasizes the goods generated by ecosystems, but neglects to ascertain habitat functions yielded by wild plants and animals, and is not a comprehensive ecosystem measure. In our study, we used three ecosystem indicators, primary productivity, bird species richness, and relative humidity to comprehensively reflect ecosystem function (Loreau 2001). Primary productivity evaluates an ecosystem’s capacity to produce biomass, bird species richness emphasizes the effects of ecosystem structure and the overall species pool on bird biodiversity, and relative humidity is an indicator of ecosystem capacity for balancing habitat conditions, stability, and service functions (Zhang & Zhang 2007, Cadotte et al. 2009). Former studies suggest these three indicators characterize ecosystem function better than any single index. A similar response model for the three indicators to seed plant richness in mainland China ensured that the relationship of biodiversity and ecosystem functions was linear when the biodiversity gradient was comparatively short, and a hump-back model when the biodiversity gradient was at a larger scale.
Acknowledgements
The study was financially supported by the National Natural Science Foundation of China (Grants No. 31170494 and No. 30870399).
ReferencesAkaike HA new look at the statistical model identification. IEEE Transactions on Automatic Control 19 (6): 716-723.1974Begon M, Harper JL, Townsend CREcology: individual, population and communities. Blackwell Scientific, Boston, Massachusetts, USA.1996Cadotte M, Cavender-Bares J, Tilman DUsing phylogenetic, functional and trait diversity to understand patterns of plant community productivity. PLoS One 4 (5): 843-845.2009Chapin FS, Zavaleta ES, Eviners VT, Naylor RL, Vitousek PM, Reynolds HLConsequences of changing biodiversity. Nature 405: 234-242.2000Chase JM, Leibold MASpatial scale dictates the productivity-biodiversity relationship. Nature 416: 427-430.2002Chen LZBiodiversity in China. Science Press, Beijing, China.1993de Groot RS, Wilson MA, Boumans RMJA typology for the classification, description and valuation of ecosystem functions, goods and services. Ecological Economics 41 (3): 393-408.2002Díaz S, Cabido MVive la différence: plant functional diversity matters to ecosystem processes. Trends in Ecology and Evolution 16 (11): 646-655.2001Editor Committee of Fauna SinicaFauna Sinica - Aves, vol. 1-14. Science Press, Beijing, China.1998Foster B, Dickson TL, Murphy CA, Karel IS, Smith VHPropagule pools mediate community assembly and diversity-ecosystem regulation along a grassland productivity gradient. Journal of Ecology 92: 435-449.2004Gillman LN, Wright SDThe influence of productivity on the species richness of plants: a critical assessment. Ecology 87: 1234-1243.2006Grace JBThe factors controlling species density in herbaceous plant communities: an assessment. Perspectives in Plant Ecology, Evolution and Systematics 2: 1-28.1999Grace JB, Anderson TM, Smith MD, Seabloom E, Andelman SJ, Meche G, Weiher E, Allain LK, Jutila H, Sankaran M, Knops J, Ritchie M, Willig MRDoes species diversity limit productivity in natural grassland communities? Ecology Letters 10 (8): 680-689.2007Grime JPBiodiversity and ecosystem function: the debate deepens. Science 277: 1260-1261.1997Hector A, Schmid B, Beierkuhnlein C, Caldeira MC, Diemer M, Dimitrakopoulos PGPlant diversity and productivity experiments in European grasslands. Science 286: 1123-1127.1999Hillebrand H, Matthiessen BBiodiversity in a complex world: consolidation and progress in functional biodiversity research. Ecology Letters 12: 1405-1419.2009Hooper DU, Chapin FS, Ewel JJ, Hector A, Inchausti P, Lavorel SEffects of biodiversity on ecosystem functioning: a consensus of current knowledge. Ecological Monograph 75: 3-35.2005Hu HJ, Jiang ZG, Wang ZWCorrelation between avian and mammalian species richness in different geographic scales in China. Chinese Biodiversity 9 (2): 95-101.2001Hughes JB, Petchey OLMerging perspectives on biodiversity and ecosystem functioning. Trends in Ecology and Evolution 16 (5): 222-223.2001Huston MBiological diversity, soils, and economics. Science 262: 1676-1680.1993Keddy PPutting the plants back into plant ecology: six pragmatic models for understanding and conserving plant diversity. Annals of Botany 96: 177-189.2005Li DQ, Sun CY, Zhang XSModelling the net primary productivity of the natural potential vegetation in China. Acta Botanica Sinica 40 (6): 560-566.1998Lieth H, Box EEvapotranspiration and primary productivity. Publication in Climatology 25 (2): 37-46.1972Loreau MBiodiversity and ecosystem functioning: current knowledge and future challenges. Science 294: 804-808.2001Mittelbach GG, Steiner CF, Scheiner SM, Gross KL, Reynolds HL, Waide RB, Willig MR, Dodson SI, Gough LWhat is the observed relationship between species richness and productivity? Ecology 82: 2381-2396.2001Northwest Research Institute of BotanyFlora of Loess Plateau. Science Press, Beijing, China.2000Partel M, Laanisto L, Zobel MContrasting plant productivity-diversity relationships across latitude: the role of evolutionary history. Ecology 88 (5): 1091-1097.2007Smith MD, Wilcox JC, Kelly T, Knapp AKDominance not richness determines invasibility of tallgrass prairie. Oikos 106: 253-262.2004Spehn EM, Hector A, Joshi J, Scherer-Lorenzen M, Schmid B, Bazeley-White EEcosystem effects of biodiversity manipulations in European grasslands. Ecological Monograph 75: 37-63.2005Suding KN, Miller AE, Bechtold H, Bowman WDThe consequence of species loss on ecosystem nitrogen cycling depends on community compensation. Oecologia 149: 141-149.2006Thompson K, Askew AP, Grime JP, Dunnett NP, Willis AJBiodiversity, ecosystem function and plant traits in mature and immature plant communities. Functional Ecology 19: 355-358.2005Tilman D, Reich PB, Knops J, Wedin D, Mielke T, Lehman CDiversity and productivity in a long-term grassland experiment. Science 294: 843-845.2001Waide RB, Willig MR, Steiner CFThe relationship between productivity and species richness. Annual Review of Ecology and Systematics 30: 257-300.1999Wu ZYVegetation of China. Science Press, Beijing, China.1980Wu ZYFlora of Tibet. Science Press, Beijing, China.1987Wu ZYFlora of China, electronic version. Vol.5-82. Science Press, Beijing, China.2004Zhang JTA study on relations of vegetation, climate and soils in Shanxi province, China. Plant Ecology 162: 23-31.2002Zhang JT, Dong YRFactors affecting species diversity of plant communities and the restoration process in the loess area of China. Ecological Engineering 36: 345-350.2010Zhang JT, Ru WM, Li BRelationships between vegetation and climate on the Loess Plateau in China. Folia Geobotanica 41: 151-163.2006Zhang JT, Zhang FDiversity and composition of plant functional groups in mountain forests of the Lishan Nature Reserve, North China. Botanical Studies 48: 339-348.2007Zheng GMClassification and distribution of birds in China. Science Press, Beijing, China.2005
Regression analysis between primary productivity of ecosystem and species richness of seed plant in the mainland of China. (A): quadratic model; (B): linear model.
Regression analysis between species richness of bird and species richness of seed plant in the mainland of China. (A): quadratic model; (B): linear model.
Regression analysis between annual mean relative humidity of ecosystem and species richness of seed plant in the mainland of China. (A): quadratic model; (B): linear model.
Comparison of regression results of linear and quadratic models in the analysis of relations between primary productivity of ecosystem, species richness of bird and relative humidity with species richness of seed plant with and without the three highest quadrats in the mainland of China. (x): refers to species richness of seed plant; (Y): refers to primary productivity of ecosystem, species richness of bird and relative humidity respectively; (AIC): 2k - 2 ln(L), where k is the number of parameters in the regression model, and L is the maximized value of the likelihood function for the regression model.