Figures
Abstract
Aim
Understanding and predicting ecosystem functioning such as biomass accumulation requires an accurate assessment of large-scale patterns of biomass distribution and partitioning in relation to climatic and soil environments.
Methods
We sampled above- and belowground biomass from 26 sites spanning 1500 km in Inner Mongolian grasslands, compared the difference in aboveground, belowground biomass and below-aboveground biomass ratio (AGB, BGB, and B/A, respectively) among meadow steppe, typical steppe, and desert steppe types. The relationships between AGB, BGB, B/A and climatic and soil environments were then examined.
Results
We found that AGB and BGB differed significantly among three types of grasslands while B/A did not differ. Structural equation model analyses indicated that mean annual precipitation was the strongest positive driver for AGB and BGB. AGB was also positively associated with soil organic carbon, whereas B/A was positively associated with total soil nitrogen.
Conclusions
These results indicated that precipitation positively influence plant production in Inner Mongolian grasslands. Contrary to the prediction from the optimal partitioning hypothesis, biomass allocation to belowground increased with soil total nitrogen, suggesting that more productive sites may increase belowground allocation as an adaptive strategy to potentially high fire frequencies.
Citation: Kang M, Dai C, Ji W, Jiang Y, Yuan Z, Chen HYH (2013) Biomass and Its Allocation in Relation to Temperature, Precipitation, and Soil Nutrients in Inner Mongolia Grasslands, China. PLoS ONE 8(7): e69561. https://doi.org/10.1371/journal.pone.0069561
Editor: James F. Cahill, University of Alberta, Canada
Received: February 23, 2013; Accepted: June 10, 2013; Published: July 25, 2013
Copyright: © 2013 Kang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This research was supported by the Project (2008-KF-03) funded by the State Key Laboratory of Earth Surface Processes and Resource Ecology, Beijing Normal University, the National Science Foundation of China (number 41271059), and Natural Science and Engineering Council of Canada (DG 283336-09). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: Han Y. H. Chen is a PLOS ONE Board member. This does not alter the authors’ adherence to all the PLOS ONE policies on sharing data and materials.
Introduction
Grassland ecosystems, accounting for about a quarter of global land area and 10 percent of global carbon storage, play an important role in regulating the global carbon cycle [1], [2]. The area of Chinese grasslands is approximately 1.18×108 hectares [3], accounting for approximately 8 percent of total area of global grasslands [4]. Of the Eurasian grasslands, Inner Mongolian grasslands are the vastest.
Previous studies have shown that grassland biomass is influenced by both climate and soil characteristics [5]. However, studies about the influence of climate on grassland biomass have so far focused mostly on aboveground biomass (AGB). By contrast, belowground biomass (BGB), which accounts for a substantially higher portion of total ecosystem biomass in grasslands [6], [7] is insufficiently studied [8], largely because of lack of a simple and efficient method to accurately determine BGB [7], [9]. In particular, our understanding of soil influence on BGB and belowground biomass allocation is limited. For Chinese temperate grasslands, several studies indicate that biomass is strongly influenced by precipitation [10], [11]. However, the response of biomass and its allocation to multiple environmental drivers, especially soil conditions, remains unclear.
Biomass allocation, typically assessed by belowground to aboveground biomass ratio (B/A) at the ecosystem level or shoot to root ratio at the individual plant level [12], reflects the strategy of plants or ecosystems to partition photosynthate in belowground and aboveground tissues [8]. Plants or ecosystems typically increase B/A ratio to take up limited belowground resources under nutrient-poor and/or water deficit conditions [13]–[16]. The idea has received support from global meta-analyses that consider a wide range of terrestrial ecosystems, where vegetation types co-vary with resource availability [12], [17]. However, the B/A ratio has been reported to be insensitive to climatic variations in Chinese grasslands [10], [11].
In this study, we examined the patterns of AGB, BGB and B/A along temperature and precipitation gradients in Inner Mongolian grasslands. We further examined whether patterns of AGB, BGB and B/A are associated with soil resource availability along temperature and precipitation gradients using structural equation models. To understand casual relationships between AGB, BGB or B/A and climate and soil variables, we used structural equation models to account for direct, indirect, and total effects of one variable on another [18]–[22].
Materials and Methods
Study Area and Sites
The study area is located in the middle and east of Inner Mongolia, the core part of the Inner Mongolian grasslands. The longitude extends from 112°40′ E to 121°10′ E and the latitude from 42°48′ N to 49°20′ N (Fig. 1). No specific permissions are required for our conducting field survey in this area, since land in China belongs to the public and our field studies did not involve any endangered or protected species within. The transect line extended from the Ulan Qab Plateau, through the Xilin Gol Grassland and the Hulun Buir Prairie to Horqin Grassland. This area is dominated by arid and semi-arid temperate continental to continental monsoon climate, which is characterized by cold, long winters and hot, rainy, and short summers. From southwest to northeast, mean annual temperature (MAT) ranged from −2.4 to 5.1°C, and mean annual precipitation (MAP) ranged from 176 to 376 mm. The soil types in these sites included sierozem, brown calcic soil, kastanozem (chestnut soil), chernozem, dark brown soil, and black soil [23]. According to China’s vegetation classification system, the grassland was divided into three types: desert steppe, typical steppe, and meadow steppe [24]. The desert steppe was dominated by Stipa glareosa and S. klemenzii. The dominant species in the typical steppe were S. grandis, S. krylovii, and Artemisia frigida, whereas in the meadow steppe S. Baicalensis, Leymus chinensis and Filifolium sibiricum were dominant species.
Data Collection
Samples were collected in August 2010, when grassland biomass was at its peak in Inner Mongolian grasslands. In order to sample a wide range of climates and soil conditions, we selected 26 sites from Sonid Right Banner in the southwest to Hulun Buir in the northeast using a transect with approximately 50 km apart between adjacent sites. The distance between the Hulun Buir Prairie and the Xilin Gol Steppe sites, separated by the Greater Khingan Range was about 200 km. The selected sites were characterized by flat relief, protected for at least 10 years by fences constructed by local governments to monitor long-term ecosystem dynamics without human disturbances such as grazing and mowing. In each site, three plots with an area 15×15 m2 were randomly assigned with a distance of at least 20 m between the two adjacent plots.
Aboveground Biomass (AGB).
Four 0.5×0.5 m2 quadrats were randomly placed within each plot. The constituent species were recorded, and their height and cover were respectively measured. All plants were clipped at the soil surface and taken to the laboratory. After removing sand and gravel, the biomass samples were oven-dried at 65°C to a constant mass and then weighed.
Belowground Biomass (BGB).
Roots were extracted with a soil corer (8 cm in diameter) at 10 cm intervals to a depth of 30 cm using a power auger, a similar method used by Brassard et al. [25]. In each plot, four cores were taken, with a total number of 12 cores in each site. The core samples were soaked fully in the water to remove soil. After washing through a 0.2 mm mesh sieve, roots were put in paper bags, taken to the laboratory, oven-dried at 65°C to a constant mass and then weighed.
Soil characteristics.
In each site, three 30 cm-deep soil profiles were excavated, each separated into three layers with a depth of 10 cm to collect soil samples. After being air-dried in a cool, well-ventilated place, the samples were passed through a 1 mm sieve and roots were removed. Soil organic carbon (SOC) and total nitrogen (TN) were measured using Vario EI elemental analyzer (Elementar Company Inc., Hanau, Germany). Potentiometric method was used to measure pH. Soil characteristics were averaged for each site in final analysis.
Climate data.
Climatic variables including mean annual temperature (MAT) and mean annual precipitation (MAP) from 1960 to 2009 were obtained from 58 meteorological stations located in or around sample sites from China Meteorological Administration. The climate data were interpolated by Kriging interpolation method using ArcGIS10.0 (ESRI Company Inc., Redlands, California, USA) for each of our 26 sample plots according to their geographic coordinates (latitude, longitude, and elevation).
Data Analysis
Differences in AGB, BGB, total biomass (i.e., sum of AGB and BGB), and B/A among three grassland types were tested by one-way analysis of variance (ANOVA), followed by Tukey’s post hoc test when grassland type was significant. The assumption of normal distribution was met according to Shapiro-Wilk test. A logarithmic transformation was needed for AGB, BGB, and TB to meet the assumptions of normality and homogenous variances. Because inherently complex causal connections exist among environmental variables in natural environments [18]–[22], we developed structural equation models (SEM) to determine the direct, indirect, and total effects of environmental variables on AGB, BGB, or B/A. For each SEM model, we hypothesized paths between endogenous variable (i.e., AGB, BGB or B/A) and exogenous variables (i.e., MAT, MAP, TN, SOC, and pH) and casual paths of climatic variables on soil variables. We specified correlations between MAT and MAP and among soil variables (TN, SOC, and pH). To better meet SEM’s assumption of linearity between dependent and independent variables, we applied logarithmic transformation to both dependent and independent variables. The linearity assumption was verified by plotting the residuals after fitting linear regression between the log-transformed dependent variable and each of log-transformed independent variable [26]. Since the values in MAT were smaller than zero in a few sites, all MAT values were transformed by adding 3 to eliminate negative values. Statistical analyses were carried out by using the SPSS 19.0 package (IBM Company Inc., Armonk, New York, USA) and structural equation models were analyzed using AMOS package (expansion pack of SPSS).
Results
AGB, BGB and B/A among Grassland Types
AGB, BGB and TB differed significantly among the three grassland types (Table 1). AGB, BGB and TB decreased from meadow steppe to typical steppe and to desert steppe (P<0.001). Mean AGB, BGB and TB of meadow steppe were 228.7, 2511.7, and 2740.5 g·m−2, respectively (Table 1). In typical steppe, they were 162.8, 1556.1, and 1718.9 g·m−2, respectively (Table 1). In desert steppe, they were 28.5, 240.7, and 269.2 g⋅m−2, respectively. On average, AGB, BGB and TB of Inner Mongolia grasslands were 154.8, 1537.5, 1692.3 g·m−2, respectively (Table 1). However, B/A ratio, with an average value of 10.6, tended to increase from meadow, typical, to desert steppe (Table 1), but did not differ significantly among the three steppe types (P = 0.47). With all data pooled, AGB and BGB were strongly correlated, as the relationship between BGB and AGB (r2 = 0.78, P<0.001, Fig. 2).
The fitted regression is y = 1.064+0964x (R2 = 0.77, P<0.001).
Casual Effects of Climate and Soil Variables on AGB, BGB and B/A
Climatic variables, MAT and MAP, were negatively correlated (Fig. 3, r = −0.66, P<0.001). Soil variables were also strongly correlated (Fig. 3). MAP affected all soil variables: direct effects of MAP on SOC, TN, and pH scored at 0.52 (P<0.01), 0.57, (P<0.01); and −0.42 (P<0.001), respectively (Table 2, Fig. 3). MAT also had direct effects on pH (path coefficient = 0.56, P<0.001) (Table 2). Additionally, there was a strong positive correlation between SOC and TN (r = 0.90, P<0.001), and significant negative correlations were found between pH and SOC (r = −0.48, P<0.05) and between pH and TN (r = −0.55, P<0.01). Climatic variables strongly influenced soil variables (Fig. 3).
Each arrow represents a direct linear causal relationship. The arcs show the correlation between two variables. Values on arrows are path coefficients and on arcs are standardized correlation coefficients. Italic values indicate the coefficients are significant at P<0.05 (*), P<0.01 (**), and P<0.001 (***). The coefficients that are not statistically significant are shown by dashed arrows. All values are log-log transformed.
MAP had the strongest positive effects on both AGB and BGB among all predictors (Table 2, Figs. 3a, b), whereas no significant effect of MAT on AGB or BGB was found. AGB and BGB responded similarly to MAP. Standardized total effect of MAP on AGB was 0.89, consisting of direct effect (path coefficient = 0.65, P<0.001) and indirect effects through soil variables (path coefficient = 0.89–0.65 = 0.24). The total effect of MAP on BGB was 0.91, consisting of direct effect (path coefficient = 0.62, P<0.001) and indirect effects through soil variables (path coefficient = 0.91–0.62 = 0.29). Among soil variables, only SOC had a significant positive effect on AGB (path coefficient = 0.49, P<0.05) (Fig. 3a). The structural equation models explained 89% and 91% of the variation in AGB and BGB, respectively.
The structural equation model explained 37% of the variation in B/A (Fig. 3c). Among soil variables, only TN had a significant positive effect on B/A (Standardized direct effect of TN on B/A, path coefficient = 1.27, P<0.01). However, neither direct nor indirect effects of MAT and MAP on B/A were significant.
Discussion
The present study is, to our knowledge, the first to demonstrate that both climatic and soil characteristics simultaneously affect both above- and belowground biomass in Inner Mongolian grasslands, China. Furthermore, we show that biomass allocation to below- and aboveground did not change along the climate gradient in Inner Mongolian grasslands, but positively associated with total nitrogen availability along the studied environmental gradients.
Variation in Biomass and its Allocation
The mean AGB and BGB of our study area (154.84 g m−2 and 1537.49 g m−2, respectively) were higher than those reported by Yang et al. [11] in the same study area (AGB = 116.6 g m−2, BGB = 553.3 g m−2). The higher biomass observed in this study than that by Yang et al. [11] may be attributed to the fact that our sample sites were protected from human disturbances such as grazing and mowing, and also that the proportion of desert sites in our study was smaller. Our mean BGB was similar to, but our mean AGB was smaller than the respective global averages reported by Jackson et al. [6]. Consequently, our B/A (10.6) was higher than the global mean [3.7 reported by Jackson et al. [6] and 4.5 reported by Mokany et al. [12]] and the overall mean of Chinese grasslands [5.7 reported by Yang et al. [11]]. However, the B/A in our study was within the top range of those reported by Jackson et al. [6] and Mokany et al. [12], and was lower than the values reported by a previous study in Inner Mongolia grasslands [10].
Relationships of Biomass and Environmental Factors
Our structural equation models explained about 90% of the variation in AGB and BGB, supporting the idea that water, temperature, and soil are the main environmental factors influencing grassland biomass [27]. Our results show that precipitation had the strongest effect on AGB and BGB. Although temperature, a surrogate of energy available to plants, is positively associated with production at a global scale [28], our results show no significant influence of temperature on AGB and BGB in our study area, indicating that there is a tradeoff effect between temperature and precipitation on plant growth. Our findings reinforce that precipitation is a major limiting factor influencing biomass in arid and semi-arid grassland ecosystems [29], [30].
Soil nutrients, in particular nitrogen, have been found to be a limiting factor for both above- and belowground production in most terrestrial ecosystems [31]–[33]. Nevertheless, others have reported limited or no influence of soil nitrogen on ecosystem production [34]–[36]. It appears that the different responses of plants to increasing nitrogen content are a result of different nitrogen deficiency levels of the local systems. However, when TN and SOC were treated as casual predictors in our structural equation models, only SOC had a positive effect on aboveground biomass.
Our results show support for the notion that plant growth in natural grasslands is primarily limited by water, resulting in high amounts of biomass allocation to belowground in order to capture these resources [37], [38]. The effect of nitrogen availability on grassland biomass is conditional to water availability [39]. Dry regions with low water availability such as Inner Mongolian grasslands may also prevent nitrogen from becoming available, as nutrients are not in solution [40].
Factors Influencing Biomass Allocation
We showed a complex web of causations for B/A in Inner Mongolian grasslands. However, our results indicated that B/A did not respond to climatic gradients, but responded positively to soil nitrogen. These results are in disagreement with the prediction from the optimal partitioning hypothesis [41]. This disagreement can be resulted from several reasons. First, despite the large spatial scale in our study and other similar studies in dry regions [10], [11], [42], the climatic variations in these studies represent only a small fraction of the global climatic variations that are considered in global meta-analyses [12], [17]. Therefore, the differences in ecosystem-specific results and global syntheses reflect different ecological scales. Second, Mccarthy and Enquist [43] suggest that the patterns of optimal partitioning can vary from species to species, obscuring large-scale patterns such as our study. The positive association between B/A and soil nitrogen may reflect a shift in plant composition along the soil nitrogen gradient [44]. Third, mean annual precipitation and temperature are negatively correlated along the climate gradients in our study. The biomass in dry and warm sites may not be sufficient to support frequent fires, whereas sites with more precipitation may have high fire frequencies due to more frequent lightning ignitions and available biomass fuels. As hypothesized by Bhattachan et al. [42], high allocation to belowground due to the need for additional storage as the risk of fire increases may be an important life-history strategy for sites with high percipitation and high fire frequencies.
Conclusions
Our results show that precipitation strongly affects aboveground and belowground biomass in semi-arid ecosystems. Moreover, our results show that belowground to aboveground biomass ratio is positively associated with soil total nitrogen, but this ratio is not related to climatic variables. Future work shall attempt to partition the influences of environmental variations and vegetation types on biomass allocation.
Supporting Information
File S1.
Table S1, Total, direct, and indirect effects of mean annual temperature (MAT), mean annual precipitation (MAP), soil organic carbon (SOC), soil total nitrogen, and soil pH on aboveground biomass (AGB), belowground biomass (BGB), and below- to aboveground biomass ratio (B/A); Table S2, Effects of MAT and MAP on soil variables. Abbreviations are same as in Table S1; Table S3, Correlations among exogenous variables. Abbreviations are same as in Table S1.
https://doi.org/10.1371/journal.pone.0069561.s001
(DOC)
Author Contributions
Conceived and designed the experiments: MYK HYHC ZYY YJ. Performed the experiments: CD WYJ. Analyzed the data: CD WYJ. Contributed reagents/materials/analysis tools: MYK YJ. Wrote the paper: MYK CD WYJ YJ ZYY HYHC.
References
- 1. Shantz H (1954) The place of grasslands in the Earth’s cover. Ecology 35: 143–145.
- 2. Scurlock JMO, Johnson K, Olson RJ (2002) Estimating net primary productivity from grassland biomass dynamics measurements. Global Change Biology 8: 736–753.
- 3. Ma W, Yang Y, He J, Zeng H, Fang J (2008) Above- and below-ground biomass in relation to environmental factors in temperate grasslands, Inner Mongolia. Science in China Series C-Life Sciences 51: 263–270.
- 4. Ni J (2002) Carbon storage in grasslands of China. Journal of Arid Environments 50: 205–218.
- 5. Schenka HJ, Jackson RB (2005) Mapping the global distribution of deep roots in relation to climate and soil characteristics. Geoderma 126: 129–140.
- 6. Jackson RB, Canadell J, Ehleringer JR, Mooney HA, Sala OE, et al. (1996) A global analysis of root distributions for terrestrial biomes. Oecologia 108: 389–411.
- 7. McNaughton SJ, Banyikwa FF, McNaughton MM (1998) Root Biomass and Productivity in a Grazing Ecosystem: The Serengeti. Ecology 79: 587–592.
- 8. Titlyanova AA, Romanova IP, Kosykh NP, Mironycheva-Tokareva NP (1999) Pattern and process in above-ground and below-ground components of grassland ecosystems. Journal of Vegetation Science 10: 307–320.
- 9. Vogt KA, Vogt DJ, Bloomfield J (1998) Analysis of some direct and indirect methods for estimating root biomass and production of forests at an ecosystem level. Plant and Soil 200: 71–89.
- 10. Han B, Fan JW, Ping ZH (2006) Grassland biomass of communities along gradients of the Inner Mongolia Grassland transect. Journal of Plant Ecology 35: 553–562.
- 11. Yang YH, Fang JY, Ma WH, Guo DL, Mohammat A (2010) Large-scale pattern of biomass partitioning across China’s grasslands. Global Ecology and Biogeography 19: 268–277.
- 12. Mokany K, Raison RJ, Prokushkin AS (2006) Critical analysis of root: shoot ratios in terrestrial biomes. Global Change Biology 12: 84–96.
- 13. Friedlingstein P, Joel G, Field CB, Fung IY (1999) Toward an allocation scheme for global terrestrial carbon models. Global Change Biology 5: 755–770.
- 14. Wilson JB (1988) Shoot Competition and Root Competition. Journal of Applied Ecology 25: 279–296.
- 15. Gedroc JJ, McConnaughay KDM, Coleman JS (1996) Plasticity in root shoot partitioning: Optimal, ontogenetic, or both? Functional Ecology 10: 44–50.
- 16. Vinton MA, Burke IC (1997) Contingent effects of plant species on soils along a regional moisture gradient in the Great Plains. Oecologia 110: 393–402.
- 17. Schenk HJ, Jackson RB (2002) Rooting depths, lateral root spreads and below-ground/above-ground allometries of plants in water-limited ecosystems. Journal of Ecology 90: 480–494.
- 18. Grace JB, Anderson TM, Olff H, Scheiner SM (2010) On the specification of structural equation models for ecological systems. Ecological Monographs 80: 67–87.
- 19. Stomp M, Huisman J, Mittelbach GG, Litchman E, Klausmeier CA (2011) Large-scale biodiversity patterns in freshwater phytoplankton. Ecology 92: 2096–2107.
- 20. Clough Y (2012) A generalized approach to modeling and estimating indirect effects in ecology. Ecology 93: 1809–1815.
- 21. Matthews B, Pomati F (2012) Reversal in the relationship between species richness and turnover in a phytoplankton community. Ecology 93: 2435–2447.
- 22. Laliberte E, Tylianakis JM (2012) Cascading effects of long-term land-use changes on plant traits and ecosystem functioning. Ecology 93: 145–155.
- 23.
Gong Z (2001) Chinese Soil Taxonomy: Science Press, Beijing, China.
- 24.
Chinese Academy of Sciences (2001) Vegetation Atlas of China: Science Press, Beijing, China.
- 25. Brassard BW, Chen HYH, Cavard X, Laganiere J, Reich PB, et al. (2013) Tree species diversity increases fine root productivity through increased soil volume filling. Journal of Ecology 101: 210–219.
- 26. Harrell F, Lee KL, Mark DB (1996) Tutorial in biostatistics multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Statistics in medicine 15: 361–387.
- 27.
Coupland RT (1979). Grassland ecosystems of the world: analysis of grasslands and their uses. New York: Cambridge University Press. 41–45.
- 28. Melillo JM, McGuire AD, Kicklighter DW, Moore B, Vorosmarty CJ, et al. (1993) Global climate change and terrestrial net primary production. Nature 363: 234–240.
- 29. Epstein HE, Burke IC, Lauenroth WK (2002) Regional Patterns of Decomposition and Primary Production Rates in the U.S. Great Plains. Ecology 83: 320–327.
- 30. Munkhtsetseg E, Kimura R, Wang J, Shinoda M (2007) Pasture yield response to precipitation and high temperature in Mongolia. Journal of Arid Environments 70: 94–110.
- 31. LeBauer DS, Treseder KK (2008) Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed. Ecology 89: 371–379.
- 32. Yuan ZY, Chen HYH (2010) Fine root biomass, production, turnover rates, and nutrient contents in boreal forest ecosystems in relation to species, climate, fertility, and stand age: literature review and meta-analyses. Critical Reviews in Plant Sciences 29: 204–221.
- 33. Yuan ZY, Chen HY (2012) A global analysis of fine root production as affected by soil nitrogen and phosphorus. Proc Biol Sci 279: 3796–3802.
- 34. Lee EH, Tingey DT, Beedlow PA, Johnson MG, Burdick CA (2007) Relating fine root biomass to soil and climate conditions in the Pacific Northwest. Forest Ecology and Management 242: 195–208.
- 35. Darby FA, Turner RE (2008) Below- and aboveground biomass of Spartina alterniflora: Response to nutrient addition in a Louisiana salt marsh. Estuaries and Coasts 31: 326–334.
- 36. Heggenstaller AH, Moore KJ, Liebman M, Anex RP (2009) Nitrogen Influences Biomass and Nutrient Partitioning by Perennial, Warm-Season Grasses. Agronomy Journal 101: 1363–1371.
- 37. Clair SBS, Sudderth EA, Fischer ML, Torn MS, Stuart SA, et al. (2009) Soil drying and nitrogen availability modulate carbon and water exchange over a range of annual precipitation totals and grassland vegetation types. Global Change Biol 15: 3018–3030.
- 38. Gill RA, Kelly RH, Parton WJ, Day KA, Jackson RB, et al. (2002) Using simple environmental variables to estimate below-ground productivity in grasslands. Global Ecology and Biogeography 11: 79–86.
- 39.
Xia J, Wan S (2008) Global response patterns of terrestrial plant species to nitrogen addition. New phytol 179.
- 40.
Lee M, Manning P, Rist J, Power SA, Marsh C (2010) A global comparison of grassland biomass responses to CO2 and nitrogen enrichment. Phil Trans R Soc Lond B Biol Sci 365.
- 41. Bloom AJ, Chapin FS, Mooney HA (1985) Resource limitation in plants - an economic analogy. Annual Review of Ecology and Systematics 16: 363–392.
- 42. Bhattachan A, Tatlhego M, Dintwe K, O’Donnell F, Caylor KK, et al. (2012) Evaluating ecohydrological theories of woody root distribution in the Kalahari. PLoS One 7: e33996.
- 43. Mccarthy MC, Enquist BJ (2007) Consistency between an allometric approach and optimal partitioning theory in global patterns of plant biomass allocation. Funct Ecol 21: 713–720.
- 44. Chen Q, Hooper DU, Lin S (2011) Shifts in species composition constrain restoration of overgrazed grassland using nitrogen fertilization in Inner Mongolian steppe, China. PLoS One 6: e16909.