Biogeographical patterns of forest biomass allocation vary by climate, soil and forest characteristics in China

To explore whether the large-scale patterns of biomass allocation vary by climate, soil, and forest characteristics in terrestrial ecosystems, on the basis of the national forest inventory data (2004–2008) and our previous field measurements (2011–2012), we investigated the variation of four biomass allocation fractions (BAFs), and their relationship with environmental factors (e.g. climate and soil chemistry) and forest characteristics (e.g. stand age and stand density) across 11 of China’s forest types. Our results revealed that BAFs have significant latitudinal, longitudinal and altitudinal trends. Stepwise multiple regression models that involve the climate, soil and forest stand properties account for a part of the biogeographical variation in BAFs, and the stand age, stand density and mean growing season temperature mainly explain these variations. Reduced major axis regression models showed that BAFs differ in their sensitivity (slope of their response to environmental gradients) to climate, soil and forest characteristics among different forest types. The results of the current study do not support the isometric allocation hypothesis, which suggests that component biomass scales equivalently as total biomass across different plant species along environmental gradients.


Introduction
Forests cover about 30% of the land surface of the Earth and 283 Gt carbons (C) are stored in global forest biomass; therefore, they play an important role in global C cycling (Dixon et al 1994, Houghton et al 2000, Piao et al 2003. When new biomass (net primary production) produces via photosynthesis, the proportional allocation of new biomass to leaves, branches, stems and roots results in the proportional distribution of standing biomass C among these biomass components (Reich et al 2014). Such processes can be influenced by plant age, resource supply, and/or climate (Poorter et al 2009, Hui et al 2012. The general lack of knowledge about biomass allocation is critical to the accuracy of global C cycle modeling and accounting (Fan et al 1998, Brown et al 1999. Root/shoot ratios (R/S) have been used to calibrate and estimate C storage from easily measurable aboveground biomass, and have also been incorporated into terrestrial ecosystem C modeling (Sack et al 2002, Scurlock et al 2002, Milchunas 2009). Compared with R/S, biomass allocation fractions (BAFs, defined as the ratios of plant biomass components (stems, branches, foliage, and roots) and total plant biomass), are biologically clear, and are well grounded in plant growth theory (Körner 1994, Poorter et al 2009. The use of BAFs retains the maximum amount of information because all ratios can be calculated from BAFs (Poorter et al 2012), and it might be useful to separate the plant bodies of larger woody species into even more categories (Körner 1994, Poorter et al 2012. Moreover, why the biomass allocation patterns vary along environmental gradients were continuously debated in China offers a unique opportunity for examining the relationship between BAFs and environmental factors across different terrestrial ecosystems in large scale because it contains complex forest characteristics (e.g. age, density, and forest types) and soil nutrients under various climatic conditions (e.g. temperature and precipitation) (Zhao and Zhou 2006). In addition, biotic and abiotic factors such as climate, soil conditions, and forest age are believed to account for a large amount of forest C stocks (Brown 2002). Therefore, a better understanding of their relative contributions to forest carbon storage is fundamentally important to make environmental policies and manage ecosystem C to enhance the forest C sink. However, despite the progress made to date, the large-scale pattern of the biomass allocation fractions (BAFs) in China's forest remains poorly defined.
In this study, we investigated the BAFs of forest leaves (LMF), branches (BMF), stems (SMF), and roots (RMF) by using a large-scale biomass survey across the forest communities in China. Our objectives were as follows: (1) to document the biogeographical patterns of BAFs in plants at a national scale, (2) to examine the sensitivity of BAFs to forest characteristics (e.g. age, density, and forest types) and soil nutrients and climatic factors (e.g. temperature and precipitation), and (3) to explore whether component biomass scales isometrically with total biomass across different forest types at the community level.

Material and methods
2.1. Large-scale forest biomass data The forest biomass, stand age and density at 1022 sites across China (figure 1) were obtained from national forest inventory data at 912 sites during from 2004 to 2008 (Luo et al 2012) and our field survey data for forest carbon at 110 random sites across southern China from 2011 to 2012 . For each site in our fields, three 20 × 50 m replication plots were first established in one sampling site for all field surveys. Then, each plot was divided into ten 10 × 10 m quadrats. Height and tree diameter at breast height (DBH) of each individual tree for stems with DBH ⩾ 2 cm and the total number of individuals in each plot were recorded. Based on the height and DBH measured above, five to seven live trees with no defect of different diameter class sizes for each species within a plot were randomly selected for cutting and weighing of their component parts (stems, branches, leaves, coarse root and fine root). Representative samples of the stems, branches, leaves and roots of standard trees were then taken back to the laboratory, dried, and used for calculation of the relationship between dry and fresh weights. Leaf (M ), National forest inventory data were also complied and tree biomass components (stems, branches, leaves, and roots) for stems with DBH ⩾ 2 cm in each site were estimated with the method described as above (Luo et al 2012). A total of 1022 average values of all tree components (leaves, branches, stems, roots, and total plants) per hectare in each site were included in the dataset. Site-related information, site conditions (longitude, latitude, altitude), forest characteristics (stand age, stand density, and forest type), soil chemistry (pH, total nitrogen (TN), and total phosphorus (TP), and climate factors (e.g. mean annual temperature (MAT), mean annual precipitation (MAP), growing season temperature (GST), growing season precipitation (GSP), potential evapotranspiration (PET, mm), and aridity index (AI, PET/MAP)) were also documented in the dataset (see appendix data S1 and table 1 for details).

Climatic variables, soil data, and forest types
Growing season temperature and precipitation are important factors that affect the large-scale pattern of grass biomass allocation (Yang et al 2010). In the current study, we selected GST, GSP, PET (calculated using the Penman-Monteith method), and AI as climate variables. GST, GSP, PET, AI latitude, longitude, and altitude were recorded at each site, and the data were used in the analyses. For records that lacked detailed altitude data, we used three-dimensional topographical maps to acquire this information based on latitude/longitude coordinates. For sites where climate variables were not recorded, the relevant data were estimated with records of 680 climatic stations in China using a Kriging extrapolation method (Luo et al 2012).
Data regarding soil pH, total nitrogen (TN), and total phosphorus (TP) were obtained from the national second soil survey and our field measurements (Han et al 2011). Forests in the dataset included all the major forest types in China, and were grouped into eleven forest types as follows: boreal/temperate Larix forest (BTLF), boreal/alpine Picea abies forest (BAPF), temperate Pinus tabulaeformis forest (TPTF), temperate/subtropical montane Populus-Betula deciduous forest (TSPF), subtropical montane Pinus yunnanensis and P. khasya forest (SPPF), subtropical Pinus massoniana forest (SPMF), subtropical montane Pinus armandii, P. taiwanensis, and P. densada forest (SMPF), subtropical Cunninghamia lanceolata forest (SCLF), subtropical evergreen broadleaved forest (SEBF), tropical rainforest and monsoon forest (TRMF), and desert riverside woodland (DRW) (Li et al 2005). Stepwise multiple regression (SMR) was used to identify the effects of climate (e.g. GST, GSP, PET, and AI), soil chemistry (e.g. pH, TN and TP) and forest characteristics (e.g. stand age and stand density) on BAFs. To demonstrate the relative effects of climate, soil and forest characteristics, partial general linear model (GLM) analyses were used to separate the variance explained by different factors into the independent effects of each individual factor and interaction effects between factors (Heikkinen et al 2005). In addition, because the reduced major axis (RMA) regression slopes of the relationship between BAFs and these variations can be used to indicate the response of biomass allocation to variation in climate, soil chemistry and forest characteristics (Sokal and Rohlf 1995), we calculated these slopes for all of the four BAFs for the 11 forest types to analyze the sensitivities of the variations in forest BAFs among forest types. Positive RMA slopes at the forest-type level indicated an increase in BAFs with increasing temperature, precipitation, soil elements etc, and vice versa. To eliminate effects of the different units of these variables on the slope values, each variable was converted to relative values (observation value/mean value) and the RMA slopes were Table 1. Site conditions (longitude, latitude, altitude), forest characteristics (stand age, stand density), soil chemistry (pH, total nitrogen (TN), and total phosphorus (TP), and climate factors (e.g., mean annual temperature (MAT), mean annual precipitation (MAP), growing season temperature (GST), growing season precipitation (GSP), potential evapotranspiration (PET, mm), and aridity index (AI, PET/MAP)) across different forest types in China. . LMF and RMF increased from west to east, whereas SMF increased from east to west. When latitudinal patterns were analyzed, BMF and RMF decreased and increased from south to north, respectively. In addition, LMF and BMF decreased significantly with increasing altitude, while SMF increased significantly with increasing altitude.

Relationship between BAFs and climatic, soil and forest characteristics
Stepwise multiple regression (SMR) analysis showed that LMF and RMF were significantly and negatively correlated with GSP and GST (table 4), respectively. SMF were significantly and positively correlated with GST, and BMF were significantly and positively correlated with both of GSP and GST (table 4). The four BAFs also exhibited different sensitivities to climate gradients (table 4). SMF and RMF responded to GST with steeper regression slopes than did GSP.
LMF significantly decreased with soil pH, soil TN, and stand age, and BMF significantly decreased with stand age (table 4). SMF significantly increased with soil TN and stand age, but decreased significantly with stand density due to lot of competition and less lateral growth (table 4). RMF significantly decreased with stand age, but increased significantly with soil pH, soil TP and stand density (table 4). In addition, SMF responded to stand age with steeper regression slope than did other BAFs, and RMF responded to soil pH with steeper regression slope than did LMF and SMF.
3.3. Sensitivities of the variations in forest BAFs among forest types SMR indicated that LMFs in BTLF, SEBF and TRMF were significantly and negatively correlated with GST, soil TP and soil pH, respectively, while that in SPPF was significantly and positively correlated with GSP (table 5). Both LMFs in TPTF and SMPF decreased with increasing stand age, and these in TPTF and SMPF decreased and increased with increasing soil TP and soil TN, respectively. LMF in TSPF was significantly and positively correlated with stand density and GSP (table 5). LMF in BAPF was significantly and positively correlated with stand density, while negatively correlated with GST, stand age, soil pH and soil TN. Likewise to LMF, the other three BAFs in the 11 forest types were also affected by different variations of climate, soil chemistry and forest characteristics (table 5).
The RMA analyses revealed that BAFs exhibited different sensitivities to climate gradients across forest types (figure 3). The sensitivity of LMF (assessed using the RMA regression slope) to GST in BTLF was the highest of all forest types. BMF in SMPF and SCLF displayed steeper slopes against the GST than those in other forest types. SMF in SPPF showed more sensitivity to GST than those in other forest types (figure 3). Similarly, LMF in TSPF, SPPF, and SCLF displayed steeper slopes against the GSP than those in other forest types. RMF in BAPF had more sensitivity to GSP than those in other forest types ( figure 3). Broadly, BAFs in BTLF, BAPF, TPTF, SPPF, SMPF and SCLF showed more sensitivity to climate than those in other forest types.
BAFs had different sensitivities to soil chemistry among the different forest types (figure 3). BAFs in BAPF, TPTF, TSPF, SMPF, SCLF and TRMF showed more sensitivity to soil pH than those in the other forest types ( figure 3). Likewise, BAFs in BTLF, SPMF and TRMF displayed shallower slopes against the soil TN or soil TP than those in the other forest types. In addition, BAFs in BAPF, TPTF, SPPF and SPMF had more sensitivity to stand age than those in other forest types (figure 3). BAFs in BTLF, BAPF, TSPF, SMPF, SPPF and SEBF displayed more sensitivity to stand density than those in other forest types.

Biogeographical patterns of forest biomass allocation
Climatically, the north-to-south and west-to-east gradients in China both reflect shifts from cold and dry to warm and moist conditions, although the thermal gradient is steeper in the former and the moisture gradient more pronounced in the latter (Han et al 2011). In this study, large variations in the biomass allocation were identified among forests in China (figure 2). LMF and SMF exhibited inverse longitudinal and altitudinal trends, but no latitudinal trends. BMF and RMF displayed a difference in latitudinal directions (table 3). The different changes in the BAFs along the geographical gradients suggest that different climate controls shape the biogeographical patterns of plant biomass allocation. Moreover, both the latitudinal and longitudinal gradients of plant biomass allocation are associated with pervasive geographical patterns in the structure and function of terrestrial ecosystems (such as functional type, biodiversity, soil property, vegetation primary production, and the ecological traits of plants) (Hedin 2004, Wright et al 2004, which themselves reflect responses to climatic conditions.

Patterns of forest biomass allocation and environmental control
The influence of climate, soil, and species composition on patterns of plant biomass allocation is complex (Enquist and Niklas 2002, Reich andOleksyn 2004).
Stepwise multiple regression (SMR) that involves the climate, soil, stand age and stand density reveals that these factors together account for a part of the biogeographical variation in BAFs (the full models are shown in  (Zhang et al 2011(Zhang et al , 2012. Our data indicated that BAFs varied markedly across different forest types (table 2), which was consistent with previous studies (Wang et al 2008). Moreover, our study also showed that BAFs had different sensitivities to climate, soil chemistry and stand characteristics among different forest types (figure 3, table 5). The integrated analysis on the responses of BAFs to multiply abiotic and biotic factors may provide comprehensive understanding on the variations of BAFs and their response of BAFs to environmental factors (Heikkinen et al 2005, Hui et al 2014. Partial general linear model (GLM) regressions could separate the variance explained by multiple factors into independent effects of all individual factors and interactive effects among factors (Legendre andLegendre 1998, Heikkinen et al 2005). The partial GLM performed in the current study revealed that the independent effects of climate were greater than those of soil chemistry and forest characteristics for BMF ( figure 4). However, the independent effects of climate on LMF, SMF, and RMF were the largest, which suggests that the geography pattern of plant biomass allocation was largely controlled by forest characteristics.
In addition, previous studies demonstrated how environmental factors affected forest biomass allocation; however, the role of interactions between environmental factors and forest characteristics has rarely been analyzed (Enquist and Niklas 2002, Hui et al 2014. In the current study, the interactions between forest characteristics, climate, and soil accounted for obviously portions of the variation in BAFs ( figure 4).
The results of the current study do not support the isometric allocation hypothesis, which suggests that component biomass scale equivalently as total biomass across a diverse range of individuals from different plant species (Euquist andNiklas 2002, Niklas 2005) and community types (Cheng and Niklas 2007). In contrast, the BAFs varied with climate, soil chemistry, and forest types (tables 4 and 5). These results were consistent with some previous studies. For example, temperature drives global patterns in forest biomass distribution in leaves, stems, and roots (Reich et al 2014). R/S in forests varied with stand age, height, shoot biomass, precipitation, and temperature (Mokany et al 2006, Wang et al 2008, Luo et al 2012. However, the RMF data in this research are also inconsistent with some previous studies. For example, the belowground biomass scaled nearly isometrically with the aboveground biomass, and the scaling exponent did not vary by tree age, density, latitude, or longitude (Cairns et al 1997, Hui et al 2014. The results of the current study can be explained in terms of optimal biomass allocation in response to environment (e.g. light, nutrients, water, and temperature). In forest communities, biomass allocation to stems, branches, and leaves all potentially increase the capture of light by trees (Dudley andSchmitt 1996, Poorter et al 2006). However, the costs (i.e. biomass allocation) and benefits (i.e. light capture) of stem, branch, and leaf growth differ (Poorter et al 2006). Nutrients and drought are the most complicated factors for the biomass allocation because there is no simple but objective way to characterize the severity of treatments (Poorter et al 2012). The response to severe drought is remarkable because plant biomass can be reduced by >50% compared with control plants, and RMF can be increased strongly (Deng et al 2006, Padilla et al 2009. Comparing for resource competition (e.g. light, nutrients, and water), low temperatures decrease the number of stems and leaves, and increase RMF (Zhou et al 2002). A variety of plant functions are also impaired by low temperature, Figure 3. Reduced major axis (RMA) regression slopes for four biomass allocation fractions (BAFs; LMF, BMF, SMF, and RMF) against climate, site conditions, soil chemistry, and stand characteristics across eleven forest types in China. RMA slopes for the eight variables (MAT, MAP, soil pH, soil TN, soil TP, altitude, stand age, and stand density) were transformed to eliminate the effects of different units (see Methods). The segmental lengths of the bars directly represent the slopes of the regression lines between BAFs and the total variables. LMF, leaf mass fraction; BMF, branch mass fraction; SMF, stem mass fraction; RMF, root mass fraction; BTLF, Boreal/temperate Larix forest; BAPF, boreal/alpine Picea abies forest; TPTF, temperate Pinus tabulaeformis forest; TSPF, temperate/ subtropical montane Populus-Betula deciduous forest; SPPF, subtropical montane Pinus yunnanensis and P. khasya forest; SPMF, subtropical Pinus massoniana forest; SMPF, subtropical montane Pinus armandii, P. taiwanensis and P. densada forest; SCLF, subtropical Cunninghamia lanceolata forest; SEBF, subtropical evergreen broadleaved forest; TRMF, tropical rainforest and monsoon forest.
In addition, biomass allocation patterns are controlled by plant types and their interaction with stand dynamics and stand conditions, e.g. age, species mixture, competition etc (Zhang et al 2010(Zhang et al , 2011. Within the group of woody species, there are different patterns of biomass allocation between deciduous and evergreen species (Poorter et al 2012). Higher allocation to roots and lower allocation to leaves in an evergreen forest will occur if higher leaf longevity is associated with infertile and/or dry sites (Givnish 2002). In a previous study, we also found that root biomass accumulation in Cyclobalanopsis glauca was more rapid in younger stands than in older stands , which is consistent with the current study (table 4). Moreover, plant biomass allocations can be affected by the interspecific or intraspecific competition as the result of the depletion of resources (Poorter et al 2012). For example, the research on plants competing for nutrients in splitroot designs has shown that RMF increases when two plants rather than one are grown in twice the volume (Gersani et al 2001, O'Brien et al 2005. Consistent with the above-mentioned studies, patterns of biomass allocation were affected by multiple environmental factors and forest characteristics in the current study (table 5). In addition to the above-mentioned influencing factors, management practices such as thinning, pruning, and fertilization might change biomass allocation in trees (Luo et al 2012. However, we did not analyze the effects of management practices in the current study, because management is a complex process with many . Variation partitioning (r 2 ) of environmental factors in accounting for the (a) LMF leaf mass fraction, (b) BMF branch mass fraction, (c) SMF stem mass fraction, and (d) RMF root mass fraction across forests in China. The symbols a, b, c represent the independent effects of climate, soil, and stand characteristics, respectively; ab is the interactive effect of climate and soil; ac, the interactive effect of climate and stand characteristics; bc, the interactive effect of soil and stand characteristics; and abc, the interactive effect of climate, soil and stand characteristics. unexplained variables. Moreover, although our data revealed that BAFs varied with climate, soil conditions, and forest types on a national level, the mechanism by which allocation is regulated remains poorly understood and source/sink relationships of all organs should be investigated in future studies (Poorter et al 2012).

Implications for estimates of C storage in China's forests
On a larger spatial scale, mean R/S is commonly calculated from published studies to quantify root biomass at the regional scale when only aboveground biomass is known (Brown et al 1993, Schroeder and Winjum 1995, Luo et al 2012. In this study, we compared forest type-specific R/S values that were calculated from the RMF of China's forests with those of global forests (Jackson et al 1996). The R/S of China's forests, which ranged from 0.11 to 0.35 (Luo et al 2012), were smaller than the global values of 0.18-0.70 (Cairns et al 1997). This suggests that China's forest C storage may be heavily overestimated if the global R/S values were used in calculations. Moreover, the large variation coefficients for BAFs suggest that the use of constant BAF values would lead to uncertain estimates of forest biomass and carbon.
The results of the current study demonstrated that high variation in BAFs was related to climate (GST and GSP), soil chemistry, stand age, stand density, and forest types (figure 3; table 4), and these predictive BAF models for specific forest types by integrating climatic, soil and forest factors may help to predict accurately regional C sequestration and cycling under climate change. Finally, we are also aware that these predictive BAF models for specific forest types by integrating climatic, soil and forest factors is the first step toward the goal of improving accurately regional C sequestration and cycling under climate change. It is still needed to incorporate these predictive BAF models into terrestrial ecosystem C modeling in further research.

Conclusions
The current study reports that the geographical variability of biomass allocation fractions (BAFs) varies by climate conditions, soil chemistry and forest characteristics in China's major forest types. Whether the component biomass scales isometrically with total biomass across different forest ecosystems and along environmental gradients was tested at the community level. Data revealed that BAFs exhibit significant latitudinal, longitudinal and altitudinal trends, which are driven by significant influences of climate, soil chemistry, and forest characteristics.
Stepwise multiple regression (SMR) models that involve the climate, soil, stand age and stand density show that these factors together account for a part of the biogeographical variation in BAFs and the explanatory power of these factors for different BAFs varied significantly. Reduced major axis (RMA) regression models indicated that BAFs differ in their environmental sensitivity (slope of their response to environmental gradients) to climate, soil and forest characteristics among the major forest types. Therefore, climate, soil and forest stand characteristics have a marked impact on biogeographical patterns of forest biomass allocation in China' forests. Our findings will help to better understand the relationship between BAFs, climate, soil chemistry, and stand characteristics, and are fundamentally important to policymaking and managing forest ecosystems to enhance the forest as a carbon sink.