Topography and life stage regulate species aboveground biomass distribution in combination in a tropical montane rainforest

Understanding the correlation between topography, species biomass and species life stage would allow forest managers to better foster carbon storage in forests. Using census data from a 60-ha plot in south China, we first quantified aboveground biomass (AGB) and how much it varied among different topographies. Next, the specific contribution of 42 dominant species to total aboveground biomass was analyzed for each of the different topographies. We also explored whether these species-topography associations changed, in terms of species’ AGB distribution, during each of the three life stages (sapling, juvenile, adult) for these 42 species. Our results showed that the average AGB was 368.79 Mg ha-1 and that it varied noticeably among the four topographies (Low valley, Slope, High valley and Ridge, which were classified by using fuzzy C-mean clustering based on elevation, convexity, and slope). AGB was significantly lower in the two valleys than in the two other topographies. Of the 42 species, 88.1% showed topographic preferences, and 78.6% appeared to exhibit topographic preferences that changed with life stage. Our work highlights the importance of topography and life stage in species biomass distribution and suggests that different combinations of species and life stages, based on species topographic preferences across life stages, may be better suited in different tropical rainforest topographies to maximize carbon storage overall.


Introduction
Globally, forests are a major reservoir of carbon, and they can act either as a carbon source when logged or disturbed or as a carbon sink when recovering post-disturbance.When forests accumulate carbon, they mitigate the global carbon emissions that cause an imbalance in the atmosphere by accumulating biomass through photosynthesis (Deal et al., 2012;Zhao et al., 2016).Therefore, forests play important roles in the regional and global carbon cycles as well as in the corresponding maintenance of biomass (Lorenz and Lal, 2010;Pan et al., 2011).Still, forests are sensitive to disturbance and varies greatly throughout heterogeneous habitats (Lorenz and Lal, 2010;Lin et al., 2012;Ratcliffe et al., 2017).
As key environmental factors, topographical factors such as elevation, convexity and slope strongly influence local-scale variation in microclimate, hydrological conditions and soil nutrient availability (Werner and Homeier, 2015;Jucker et al., 2018;Wang et al., 2023).As such, they constrain the conditions within which plant grows, driving environmental filtering, controlling the distribution of species, and influencing the biomass of forests (Mascaro et al., 2011;Riihimäki et al., 2017;Jin et al., 2018;van der Plas, 2019;Yuan et al., 2020).Many studies have demonstrated significant associations between trees and topography in forests (Gunatilleke et al., 2006;Fortunel et al., 2018;O'Brien and Escudero, 2022;Zuleta et al., 2020Zuleta et al., , 2022)).Some tree species have been found to be more or less frequent in swamp, valley, slope, ridge, or low-lying plateau habitats (Harms et al., 2001;Queenborough et al., 2007;Guo et al., 2015).However, species abundance does not necessarily relate to biomass.Information regarding how tree species-specific contributions to forest aboveground biomass vary among different topographies is scarce, though this knowledge would allow for a novel way to assess species coexistence and would provide valuable information to forest managers looking to maximize carbon storage (Lin et al, 2012;Colgan and Asner, 2014;Smith et al., 2016).
The concept of the regeneration niche suggests that speciestopography associations change with life stage because each stage has different physiological requirements and responds to different selective pressures (Grubb, 1977;Webb and Peart, 2000;Silvertown, 2004;Cui and Zheng, 2016;Yao et al., 2020).Species-topography associations change can be quantified by counting the individuals of a species at different life stages, as measured by tree DBH (diameter at the breast height) size (Comita et al., 2007;Lai et al., 2009;Kanagaral et al., 2011).Also important is detecting how these associations change based on species' aboveground biomass during different life stages.Species aboveground biomass encompasses species abundance and tree volume with wood density, but it also expresses the outcome of myriad ecological processes (De Castilho et al., 2006;Colgan and Asner, 2014).This knowledge could further be beneficial for understanding species coexistence (Colgan and Asner, 2014;Ramachandran et al., 2023) and for promoting biomass carbon gain in controlled burns.
Tropical rainforests, with their high levels of biodiversity, are among the most productive ecosystems on the earth (Field et al., 1998).Here, we utilized the 60-ha forest dynamic plot that was established in the primary tropical montane rainforest in Hainan, China to analyze the correlations between topography, species aboveground biomass (AGB) and life stage.Our objectives were 1) to quantify AGB and its variation among topographies; 2) to determine the species-topography associations based on species-specific AGB contribution; and 3) to analyze how species-topography associations change across life stages based on species AGB distribution.

Study site
This study was conducted in the 60-ha (1 000 m × 600 m) Forest Dynamics Plot (FDP) established between 2009 and 2010 in Jianfengling National Nature Reserve (JNNR: 18.33°~18.95°N,108.68°~109.20°E),Ledong County, Hainan Province, China.The site is characterized by a tropical monsoon climate, with an average annual precipitation of 2449 mm, and a mean annual temperature of 19.8° (Zhou et al., 2009).The main soil type is lateritic yellow soil with high soil moisture and humus content (Li et al., 2002).

Data collection
The Jianfengling FDP was divided into 1500 20 m × 20 m quadrats using 2464 posts.Plot altitude varies from 866.3 m to 1016.7 m, convexity varies from -7.7 to 7.2, and slope ranges from 1.7°to 49.3°.Between December 2010 and December 2011, all woody plants with DBH ≥ 1.0 cm were mapped, tagged, measured, and identified to species.Overall, a total number of 439,676 stems from 62 families, 155 genera, and 290 species were recorded (Xu et al., 2015a).The four most dominant families were Lauraceae, Fagaceae, Rubiaceae and Arecaceae, and the importance values for these families were all greater than five.The four most dominant species were Livistona saribus, Gironniera subaequalis, Cryptocarya chinensis, and Alseodaphne hainanensis.

Topography classification
We classified the Jianfengling FDP into four topographies (Figure 1; Low valley, Steep slope, Upper valley, Ridge) by using fuzzy C-mean clustering based on elevation, convexity, and slope in each of the 20 m×20 m quadrats (Wang et al., 2018; Table 1).Fuzzy C-means clustering obtains the membership degree of each sample to all cluster centers by optimizing the objective function, and then can automatically classify sample data.Among fuzzy algorithms, it is the most widely used and successful because it can reflect the real world objectively (Gong et al., 2007;Bezdek, 2013).

Aboveground biomass calculation
In this study, the AGB of each tree was estimated using the allometric model for the tropical montane rainforest in Jianfengling (Equation 1) (Zeng et al., 1997;Wang, 2017).With this equation, we calculated AGB for each tree in the Jianfengling FDP.AGB = 0:100194058 Â DBH 2:463324542 ; R 2 = 0:94 (1) Zeng et al. (1997) fitted an allometric growth equation of more than 70 tree species in Jianfengling.On this basis, Equation 1 was fitted by adding more than two hundred analytic trees.This mixedspecies allometric regression model has a high accuracy due to the large number of species considered and the large number of samples, which offset the inter-species differences.In addition, in tropical forests, the dominant species are not obvious.If the regression model of each tree species is fitted, the workload is huge.Given its convenience and high accuracy, we used Equation 1to calculate AGB.

Life stage classification
D95 0.1 represents the 95th percentile of DBH of all individuals in a species≥0.1Dmax , where D max is the largest measured DBH in a species (King et al., 2006).Previous studies have verified that D95 0.1 can be widely applied to the division of the adult stage of trees in old-growth forests (King et al., 2006;Lan et al., 2011).Therefore, individuals of each of 42 canopy species were divided into three life stages according to D95 0.1 : saplings [ (1,5) cm DBH], juveniles [(5, D95 0.1 ) cm DBH], and adults (>D95 0.1 cm DBH).

Relative contribution of dominant species to AGB in different topographies
Our analysis focused on the dominant species ranking of the top 30 species in aboveground biomass in each of the four topographies.A total of 42 canopy species were included in this analysis because some species were present in multiple topographies.For each of the 42 species, the relative contribution to total AGB was calculated at the whole plot and at the topography scale.We examined the associations between topography and the AGB distribution of these 42 species (i.e. for each of these 42 species, whether a given species contributed significantly more or less to one or more topographies) using torus-translation tests (Harms et al., 2001;Lin et al., 2012), which consider the spatial autocorrelation in both the topography and the plant distributions.This method generates a null model in which the species are distributed randomly with respect to topography, and then it tests whether the observed relative contribution of a focal species in each topography is significantly higher or lower than would be expected at random.Further details of this method can be found in Harms et al. (2001).
Torus-translation tests were also conducted for the 42 species to examine species-topography associations across different life stages (Lai et al., 2009;Oktavia and Jin, 2019).To better understand how species-topography associations changed during three life stages (sapling, juvenile, adult), we analyzed it in two aspects.On the one hand, we explored how many species showed inconsistent associations with a specific topography across three life stages and how many showed consistent associations.On the other hand, we compared the differences in the number of significant associations with each topography among the three life stages.
In addition, one-way analysis of variance model (ANOVA) and multiple comparisons were used to test differences in AGB in different topographies.To further test the effect of topography, life stages and species on AGB distribution, two-way and three-way ANOVA were also used.All the analyses in this study were conducted in R version 2.15.1 (Oksanen, 2011).

Aboveground biomass and its variation
In the Jianfengling FDP, the average basal area and average aboveground biomass were 56.04 m 2 •ha -1 and 368.79 Mg•ha -1 , respectively (Table 2).The top three species in AGB across the entire plot were Livistona saribus, Lithocarpus fenzelianus, and Gironniera subaequalis (Table 3; Figure 2).Just 42 of 290 species contributed to 80.3% of the total aboveground biomass in the Jianfengling FDP (Table 3).Average basal area and average aboveground biomass varied among the four topographies (Table 2).Steep slope and ridge had higher aboveground biomass than the low valley and upper valley did (P<0.05).In terms of basal area, the order was ridge > steep slope > high valley > low valley (P<0.05).

Topographic differences in speciesspecific contributions
42 dominant species contributed 77.5%-87.3% of total aboveground biomass in each topography (Supplementary Table S1; Figure 2).Relative contributions of species varied among different topographies (Supplementary Tables S1, S2; Figure 2).37 of the 42 species were significantly associated with one or more topographies (Supplementary Table S2; Figure 2).There were respectively 11, 4, 10, and 10 species positively associated with low valley, steep slope, upper valley, and ridge.Four species were positively associated with low valley and upper valley, and one species was positively associated with steep slope and ridge.In other words, 11, 4, 10, and 10 species contributed significantly more to low valley, steep slope, upper valley and ridge, respectively, than would be expected at random.Four species contributed significantly more to low valley and upper valley, and one species contributed significantly more to steep slope and ridge than would be expected at random.In addition, 14, 7, 7, and 14 species contributed significantly less to low valley, steep slope, upper valley and ridge, respectively, than would be expected at random.Four species contributed significantly less to low valley and upper valley, and four species contributed significantly less to steep slope and ridge than would be expected at random.The results of ANOVA indicated species had the largest effect on AGB distribution (Table 4).The effects of topography and their interaction effect were also significant.

Variability of species-topography associations across life stages
According to the torus translation tests based on species AGB distribution, 89 positive associations among 38 species and 101 negative associations among 39 species were observed across the three life stages (Figures 3, 4).In terms of positive associations, 16 out of 30 significant associations (D and G in Figure 3) between saplings and topography persisted into the juvenile stage in the same topography, and 12 of the associations (E and G in Figure 3) were sustained in the adult stage in the same topography.Correspondingly, 25 out of 35 negative associations between saplings and topography (Figure 4) persisted into the juvenile stage, and 10 of the associations were sustained in the adult stage.In summary, less than one-third of the positive associations and just over a quarter of the negative associations between saplings and topography coincided with the two later life stages (F and G in Figure 3).Of the adults, for one-third of positive associations and for less than one-fifth of negative associations, there were no corresponding associations between species and topography with the two prior life stages.Just 7 positive and 8 negative associations between species and topography (G in Figures 3, 4) were consistent across the three life stages.In total, 33 species showed inconsistent associations with a specific topography across three life stages, and 8 species showed consistent associations.In addition, the results of ANOVA indicated life stages had the largest effect on AGB distribution, followed by species (Tables 5, 6).The effects of topography and its interaction effect with species and life stages were also significant.
Based on species AGB distribution, the number of associations with each topography differed among the three life stages (Figures 5, 6).At the sapling and juvenile stages, more positive associations were found in the ridge than in the other three topographies, but at the adult stage, the most positive associations were found in the low valley.For    the saplings and juveniles, more negative associations were detected in the low valley than in the other three topographies, but for the adults more than one-half of all negative associations were detected in the ridge.

Discussion AGB and its variation among topographies
There was an average AGB of 368.8 Mg ha -1 in the Jianfengling FDP.Compared with other tropical regions, the AGB of Jianfengling tropical mountain rainforest fell within the same range of biomass variation as the mature tropical lowland rainforests of Peninsular Malaysia (365~440 Mg ha -1 ) (Hoshizaki et al., 2004).The Jianfengling biomass was significantly greater than that of the French Guiana neotropical rainforest (309 Mg ha -1 ), the Panamanian tropical rainforest (281 Mg ha -1 ) (Chave et al., 2001(Chave et al., , 2003)), and also than the average biomass of the Brazilian Amazonian neo-tropical rainforest (356 Mg ha-1), the Cambodian tropical deciduous forest (189 Mg ha -1 ), and the semi-deciduous forest (244 Mg ha -1 ) and the limestone seasonal rainforest in Xishuangbanna, China (Laurence et al., 1999;Top et al., 2004;Qi and Tang, 2008).These results emphasize that Jianfengling's primary tropical montane rainforest has high carbon storage capacity in its aboveground biomass.
In this study, AGB varied noticeably among the four topographies.Steep slope and ridge supported more AGB, while the two valleys supported less.High variation in AGB at the local scale has been demonstrated due to the effects of topography (De Castilho et al., 2006;Sattler et al., 2014;Riihimäki et al., 2017).Researchers have shown that topography can influence AGB distribution in multiple ways.Topography may affect AGB distribution by influencing local conditions, such as soil, hydrology, and light conditions (Mascaro et al., 2011;Liu et al., 2014), but also by influencing disturbance intensity (Lin et al., 2012).In this plot, the two valleys had many gullies due to high rainfall, thus there was less space available for Relative contribution to total aboveground biomass (AGB) per topography type for the ten most important species contributing to the total AGB in Jianfengling FDP.Species codes are defined in Table 3. Bars show the relative contribution of the species to the aboveground biomass in each topography type, while squares indicate the relative contribution of the species to the aboveground biomass of the entire plot; thus differences between bars and squares represent the influence of topography associations.'+' and '++' respectively indicate significant positive associations at levels of 0.05 and 0.01, and '-' and '-' respectively indicate significant negative associations at levels of 0.05 and 0.01.
water-intolerant plant growth (Punchi-Manage et al., 2013).In addition, light availability was lower in the valleys than in the steep slope and ridge topographies (Körner, 2007;Huang et al., 2014Huang et al., , 2016)).This may explain why the two valleys supported less AGB and why many light-demanding species were more abundant in the steep slope and ridge topographies.For example, Madhuca hainanensis contributed significantly more biomass to the total AGB of steep slope and ridge than to that of the two valleys (Kruskal-Wallis test, P<0.01).This suggests that available light may be one key environmental factor that affects species composition and plant AGB distribution in the Jianfengling plot.

Topographic preferences of species
Our results clearly showed that different species differed in their topographic preferences, as demonstrated by their AGB contributions.This was consistent with previous speciestopography studies in other forests based on species abundance (Lin et al., 2012;Smith et al., 2016;Zuleta et al., 2020).Species differences in ecological strategies and different resource availability throughout the topographies may contribute to this pattern of biomass distribution (Yamada et al., 2007;Wang et al., 2018).Of the 42 species, 88.1% were positively or negatively associated with one or more of the different topographies.Compared with other forest dynamic plots that have been examined using similar methods, the Jianfengling FDP has a higher percentage of species demonstrating significant topography associations (Harms et al., 2001;Gunatilleke et al., 2006;Wang et al., 2011;Kanagaral et al., 2011) and this can be attributed to two main reasons: its topographic complexities and its forest structure (Lai, 2008).On the one hand, topographic complexities enhanced the environmental filtering effect, which led to a higher percentage.
On the other hand, in contrast to a secondary forest, after a long    period of environmental adaptation, the topographic niche differentiation of species in the Jianfengling primary forest is more obvious (Letcher et al., 2012;Jiao et al., 2012).
The niche complementarity hypothesis suggests that a greater number of species take greater advantage of the niche opportunities that are available in an environment, thus resulting in better use of available resources (Cardinale, 2011).In this plot, we found that different species' topographic preferences varied based on biomass contribution, which indicates that species may complement each other in space, thus allowing for high diversity and high biomass/ carbon storage at the whole plot scale.This may be a form of niche complementarity (Lin et al., 2012), which contributes to a high aboveground in Jianfengling FDP.Combined with ANOVA analysis results, our study suggests that choosing different species combinations for different topographies would overall help maximize biomass carbon storage in tropical forests.
Variation in effects of topography among life stages 33 species (34 out of 38 positively associated and 31 out of 39 negatively associated species) showed inconsistent associations with a specific topography across the three life stages.This revealed that for most species, topographic preference and limitation were not consistent across multiple life stages.This also has been demonstrated by many previous studies based on species abundance (Comita et al., 2007;Lai et al., 2009;Punchi-Manage et al., 2013).More than half of the significant positive associations between saplings and a particular topography persisted into the juvenile stage, but few of these associations were sustained into the adult stage.These results indicate that species' ecological habitat preferences differ by developmental stage.Similarly, most of the negative associations between saplings and a particular topography were maintained in the juvenile stage, but just 23% of these sapling associations were further maintained into the adult stage.This indicates that, for most species, topographic limitation shows consistency from the sapling to the juvenile stage, but it tends to differ from the adult stage.Some species were positively or negatively associated with the same topography across all life stages at this site.For example, the heliophilous species Madhuca hainanensis associated positively with the ridge and negatively with the low valley and the high valley across all three life stages.Hydrophilic species Beilschmiedia laevis, Machilus cicatricose, and Polyosma cambodiana associated positively with either the low valley or the high valley across all three stages.These above species are especially light-demanding or waterdemanding, and consequently, the saplings, juveniles and adults Number of species positively associated with each topography type at the sapling, juvenile and adult stages in the Jianfengling FDP.remained spatially associated.Alternatively, this pattern may be attributed to seed dispersal limitation, in which saplings primarily establish near to the parent trees and thus form the same topographic associations (Webb and Peart, 2000;Lai et al., 2009).The role that these species, consistent in their topographic preferences across life stages, play in community assembly and forest carbon management requires further study.
In terms of the number of associations with each topography, more negative associations were detected at the sapling and juvenile stages in the low valley than in the other three topographies.This was probably due to lower light availability in the low valley, resulting from the shelter of the surrounding mountains and the greater density of large trees (DBH≥50 cm).On the contrary, more positive associations were found at the adult stage in the low valley than in the other three topographies, which may have resulted from the higher soil water content and the higher soil nitrogen, phosphorus and potassium (Xu et al., 2015b;Wang et al., 2018).In addition, the study site was in an area that is frequently hit by typhoons.In this case, direct damage mainly influences the tree layer, and large-diameter individuals are the most severely damaged (Xu et al., 2008).Because the ridge is often more susceptible to typhoons, the notably higher number of negative associations and the notably lower number of positive associations of species with the ridge were observed in the adult stage.
Furthermore, given the associations between species and topography across different life stages, our analysis showed that in the low valley, more adults than other life stages of Alseodaphne hainanensis and Canarium album led to a higher AGB, while more saplings and juveniles than adults of Alseodaphne hainanensis led to a lower AGB.In the ridge, more juveniles than other life stages of Dacrydium pectinatum and Lithocarpus fenzelianus contributed to a higher AGB, while more adults than other life stages did not relate to greater AGB.Combined with the results of ANOVA analysis, our results suggest that choosing different combinations of different species and life stages for different topographies is important for maximizing overall carbon storage when managing tropical montane rainforests.

Conclusion
Working in a highly diverse primary montane rainforest, we found that the average aboveground biomass was 368.79 Mg ha -1 , but it varied noticeably among the four topographies.AGB was significantly lower in the two valleys than in the two other topographies.Most of the 42 dominant species varied in their topographic preferences and further exhibited inconsistent topographic preferences across life stages.Our results strongly demonstrate the importance of topography and life stage in species biomass distribution, and we suggest that land managers consider both topography and life stage to best maximize carbon storage in forests.

FIGURE 1
FIGURE 1Map of the distribution of the four topographies of the Jianfengling 60 ha plot, as classified by fuzzy C-mean clustering in 20 m×20 m cells.

FIGURE 3
FIGURE 3 Venn diagram of positive associations of species at each life stage with all four topographies in the Jianfengling FDP.Here shown is the number of positive associations instead of the number of species.(A-C) show number of only one positive association at any life stage (P<0.05).(D-F) represent number of consistent positive associations at two stages (P<0.05),and (G) is number of consistent positive associations across all three life stages (P<0.05).

FIGURE 4
FIGURE 4Venn diagrams of negative associations of species at each life stage with all four topographies in the Jianfengling FDP.Here shown is the number of negative associations instead of the number of species.(A-C) show number of only one negative association at any life stage (P<0.05).(D-F) represent number of consistent negative associations at two stages (P<0.05),and (G) is number of consistent negative associations across all three life stages (P<0.05).

TABLE 1
The characteristic and quadrats number of 4 topographies.
The values of elevation, slope, and convexity in the 2 to 4 columns were the value of cluster centers clustered by fuzzy C-mean clustering.Mean value= Mean value ± SD (a standardized deviation).

TABLE 2
Basal area and aboveground biomass in each of the four topographies and the entire plot.
AGB, aboveground biomass.The values in the third and fourth columns are the means ± one standard error of means.The different lowercase letters indicate significant differences among the four topographies.

TABLE 3
Maximum DBH, AGB and tree number of 42 dominant species in Jianfengling FDP.

TABLE 4
Variation analysis of the effects of topography and species on AGB distribution.

TABLE 5
Variation analysis of the effects of life stages and species on AGB distribution.
Df, degree of freedom.*** represents significant effects at a level of 0.001.

TABLE 6
Variation analysis of the effects of topography, life stages and species on AGB distribution.