Topography and Soil Properties Determine Biomass and Productivity Indirectly via Community Structural and Species Diversity in Karst Forest, Southwest China

: The forest ecosystem is an important part of the terrestrial ecosystem carbon sink, and its rate of biomass accumulation inﬂuences its carbon sink potential. Therefore, it is particularly important to understand the biomass and productivity of forest ecosystems, and their driving factors, especially in karst areas with a fragile ecological environment. We established a 2 ha plot in karst forest in southwest China, and investigated species composition, community structure, topography and soil nutrients in the years 2007 and 2017. In this analysis, the correlations between tree diversity and each factor were evaluated using a Pearson correlation analysis. In addition, the relationships between soil nutrients and topographies and their effects on productivity and biomass were further evaluated, either directly or indirectly, through species and structural diversity by using a structural equation model (SEM). The results showed that the number of individuals in each species decreased, and productivity was 1.76 Mg ha − 1 yr − 1 , from 2007 to 2017. Species diversity was negatively correlated with biomass and positively correlated with productivity; structural diversity was negatively correlated with biomass and productivity, while structural diversity was negatively correlated with biomass and positively correlated with productivity. In addition, the effects of soil factors on biomass and productivity were signiﬁcantly different: TN had a signiﬁcant positive effect on productivity, while all soil factors except total nitrogen (TN) had signiﬁcant positive effects on biomass. The structural equation results also showed that topographic and soil factors can directly affect productivity; structural diversity has a strong direct negative impact on biomass, while species diversity, structural diversity and biomass have similar direct positive impacts on productivity. Structural diversity was better than species diversity when explaining biomass accumulation. In conclusion, biotic and abiotic factors both inﬂuence forest productivity in karst forests in southwest China, and improving species diversity and community structure complexity is of great signiﬁcance for forest management and productivity promotion. The research further improve the understanding of biomass and productivity in karst forest ecosystems, and their driving factors, which will provide relevant theoretical support for sustainable forest development in southwest karst. H.D. and G.H.; L.Z., H.D. and T.S.; and draft L.Z.; writing—review and editing, T.S. and G.H.; funding H.D., F.Z.


Introduction
The forest ecosystem is a natural system with a well-defined structure, function and regulation, and is the most extensive component of the terrestrial ecosystems [1,2]. Meanwhile, the forest ecosystem is also the largest carbon pool among terrestrial ecosystems, accounting for more than 3/4 of terrestrial biomass [3]. Therefore, the forest ecosystem plays an important role in regulating terrestrial carbon cycle [3,4]. Global change and human disturbance strongly affect ecosystem components, structures, functions, processes and services and their spatial patterns [5]. Forest fire disturbance is a natural ecological process of forest ecosystem, which is one of the main driving factors of ecosystem structure and function, and has long-term influence on nutrient cycling and composition of fuel types. Driven by the dual pressures of global change and human interference, biodiversity on the earth is facing serious loss, and the structure and function of forest ecosystems are seriously affected [6]. Therefore, maintaining the stability of forest ecosystem services and functions is necessary to understand the mechanisms underlying forest community dynamics [7].
The relationship between biodiversity and productivity is an important indicator to measure ecosystem dynamics, function [8,9], an important topic in sustainable forest management [10], and one of the issues of concern for global scientists [11]. A large number of studies have shown that species diversity is positively correlated with productivity [12][13][14][15]. For example, Liang et al. indicated that continued loss of biodiversity will lead to an accelerated decline in global forest productivity [13]. But negative correlation [16], a hump relationship [11], and irrelevance [17] have also been reported. Some studies have also shown that mixed forest has higher productivity than pure forest [10,12], and that increasing species diversity improves productivity [10,18]. Such results can be explained by the niche complementing effect [11]. Of course, other studies have shown that this effect cannot explain all the results: the niche complementarity effect is not only related to interspecific differences, but is also affected by a variety of biotic and abiotic factors [19]. So far, the effects of species diversity on productivity and biomass in natural forests and their maintenance mechanisms remain controversial. Therefore, it is urgent to further strengthen the research on the relationship between diversity and productivity and biomass.
Community structure diversity plays an important role in productivity [20,21], and is generally defined as one or more combinations of spatial structure distribution, species diversity, and tree size variation (such as tree height and diameter) [20,22,23]. Complex structures can cause resource heterogeneity or variable resource use efficiency, and may be the cause of different results in past studies [6,24]. Through access to environmental resources and by promoting biomass accumulation [3,4,25,26], forests can increase productivity. In adddtion, the age of a forest ecosystem greatly affects its biomass/carbon storage and production. The biomass of different forest ecosystems may reach their peaks at different ages [27]. In another word, similar size structures may have different production rates. Previous studies show that community structure diversity is positively with productivity relationship [18,21], uncorrelated or negatively correlated [16,20]. Individual differences, DBH variation index, species fragrance index, stand age density and other factors will each affect the distribution of tree species and lead to changes in forest structure [24,[28][29][30]. Overall, the effects of species diversity and community structure on biomass and productivity need further discussion and research.
In addition, the distribution pattern of plant communities is the result of climate, soil, topography and other environmental factors at different scales. The productivity of vegetation is also influenced by species composition and geographical environment [12]. Soil fertility and topography are key factors affecting plant community structure and productivity [20]. Studies have shown that high-nutrient soils can promote faster plant growth and thus have a more positive effect on tree productivity than poor soils [31], while topographic factors such as elevation, slope and aspect also influence the accumulation of biomass and changes in productivity [32,33]. Therefore, it is necessary to study the relationship between biodiversity and productivity and environmental factors in complex natural forest ecosystems.
Southwest China is the largest area of karst forests in the world, with the characteristics of weak ecosystem stability and strong heterogeneity of habitat (exposed rock, barren ground and discontinuous soil, with an extremely fragile geological background, small environmental carrying capacity, low tolerance to disturbance and establishment on rocky surfaces, and excessive soil calcium, etc.) [34,35]. Mixed evergreen-deciduous broadleaf forest is unique and representative of the karst landscape in southwest China. This forest type has a complicated community structure, rich biodiversity, and high habitat heterogeneity, making it an ideal community in which to study forest productivity.
In the early stages, large-scale productivity research is mainly carried out by remote sensing. Most small-scale biomass studies are single community surveys, and productivity research is mostly carried out in space instead of time. Many studies have explored influencing factors of karst biomass, but few studies have focused on the interaction between biomass and productivity in karst forests. In particular, biotic and abiotic factors contributing to biomass changes are less studied in this region [36][37][38]. In recent years, due to climate change and other factors, the structure of the forest ecosystem has changed, so that it is urgent to strengthen research on the relationship between biodiversity and biomass and productivity of this ecosystem. Therefore, a detailed study of karst forest productivity is of great significance for understanding the spatial distribution and characteristics of forest vegetation in China.
This study focused on karst forests in southwest China, exploring the structure and species diversity of karst forests to indirectly determine biomass and productivity through the analysis and assessment of topography and soil. The study combined multiple statistical methods, with the following objectives. (1) What are the dynamic changes of forest biomass and productivity in the primary forest 2 ha plot? (2) Species diversity or structural diversity-which has the greatest impact on biomass and productivity? (3) How do topographic and soil factors affect changes in biomass and productivity?

Study Description
The study was conducted in Mulun National Nature Reserve (hereafter, Mulun Reserve), which is among the most important refuges in China for primary and secondary karst forest flora. The reserve lies in the Huanjiang Autonomous County, northwestern Guangxi Province, China (107 • 54 01 -108 • 05 51 E, 25 • 07 01 -25 • 12 22 N), with an area of 301 km 2 ( Figure 1). It was established in 1991 to protect the subtropical mixed evergreen and deciduous broad-leaf forest ecosystem, which has developed on the limestone substrate, and contains the world's largest and best preserved karst forest.
The forest coverage rate is 94.8%, and the region has a subtropical monsoon climate with annual average temperature of 19.3 • C, average annual precipitation 1529.2 mm and frost-free period of 310 d [35]. The elevation range is 413.6 to 637.8 m. The geology is predominantly soluble and porous limestone, resulting in a karst landscape of steep hills separated by lowland depressions with numerous potholes and caves and extensive below-ground streams. Hillside rock exposure rate is high, soil is sparsely and unevenly distributed, and is mainly lime soil; vegetation coverage is very low [38].  . All the individuals of woody plant with a diameter at breast height (DBH) ≥ 1 cm were tagged, identified, measured, and georeferenced from August to October 2007, (main: tree species name, DBH, tree height, crown width, coordinates and growth status were recorded). The important forest species included Sinosideroxylon, pedunculatum, Boniodendron minus, Ligustrum quihoui, Handeliodendron bodinieri, Radermachera sinica.

Topographic Variables and Soil Variables
Topographical factors, including elevation (ELE), slope (SLO), cos (aspect) (ASP), convexity (CON), and rock exposure rate (ROR) for each quadrat were measured. ELE was obtained from a Garmin GPSmap 60CSX (Garmin Ltd., Olathe, KS, USA), ASP was obtained from an electronic tachometer (NTS-302R, SOUTH Ltd., GuangZhou, China), and slope was obtained as the mean angular deviation from horizontal of each of the four triangular planes formed by connecting three of its corners. The ROR was calculated as the ratio of rock to the total area of a grid using a 100 cm × 100 cm grid with 25 subgrids of 20 cm × 20 cm, taking the average of five points along the diagonal for each quadrat. ASP and ROR were converted to quantitative data (ASP: 0, 1, 2, 3, 4, 5, 6, 7, and 8 for a depression, northern slope, northeast slope, northwest slope, eastern slope, western slope, southeast slope, southwest slope, and southern slope, respectively) [39].
About 8-10 soil subsamples were collected at a depth of 15 cm, along an S-shaped transect, using a soil corer 3 cm in diameter. The soil subsamples were pooled into one composite sample per quadrat of 20 m × 20 m. The soil samples were air-dried, thoroughly mixed, and passed through a 2-mm sieve to remove gravel and roots [38].  All the individuals of woody plant with a diameter at breast height (DBH) ≥ 1 cm were tagged, identified, measured, and georeferenced from August to October 2007, (main: tree species name, DBH, tree height, crown width, coordinates and growth status were recorded). The important forest species included Sinosideroxylon pedunculatum, Boniodendron minus, Ligustrum quihoui, Handeliodendron bodinieri, Radermachera sinica.

Topographic Variables and Soil Variables
Topographical factors, including elevation (ELE), slope (SLO), cos (aspect) (ASP), convexity (CON), and rock exposure rate (ROR) for each quadrat were measured. ELE was obtained from a Garmin GPSmap 60CSX (Garmin Ltd., Olathe, KS, USA), ASP was obtained from an electronic tachometer (NTS-302R, SOUTH Ltd., GuangZhou, China), and slope was obtained as the mean angular deviation from horizontal of each of the four triangular planes formed by connecting three of its corners. The ROR was calculated as the ratio of rock to the total area of a grid using a 100 cm × 100 cm grid with 25 subgrids of 20 cm × 20 cm, taking the average of five points along the diagonal for each quadrat. ASP and ROR were converted to quantitative data (ASP: 0, 1, 2, 3, 4, 5, 6, 7, and 8 for a depression, northern slope, northeast slope, northwest slope, eastern slope, western slope, southeast slope, southwest slope, and southern slope, respectively) [39].
About 8-10 soil subsamples were collected at a depth of 15 cm, along an S-shaped transect, using a soil corer 3 cm in diameter. The soil subsamples were pooled into one composite sample per quadrat of 20 m × 20 m. The soil samples were air-dried, thoroughly mixed, and passed through a 2-mm sieve to remove gravel and roots [38].
Soil pH was determined in a 1:5 soil-water slurry using a combination glass electrode. Soil organic matter (SOM) was determined by the chloride potassium dichromate volumetric-external heating method. Total nitrogen (TN) was determined by the semi-micro Kjeldahl determination of nitrogen method. Total phosphorus (TP) was determined by sodium hydroxide (NaOH) fusion-ammonium molybdate spectro photometry. Total potassium (TK) was determined by NaOH fusion-flame spectrophotometry. Available nitrogen (AN) was determined by the diffusion-absorption method. Available phosphorus (AP) was determined by NaHCO 3 extraction-ammonium molybdate spectrophotometry. Available potassium (AK) was determined by (NH 4 ) 2 CO 3 extraction-flame spectrophotometry. Root content and gravel content were obtained in the process of sieving as a percentage of the total weight. To improve the precision of the results, three replicate samples were analyzed and averaged for each of the soil analyses, except for root and gravel content. Statistical data are shown in Table 1.

Species Diversity and Community Structure Diversity
Species richness (R), Shannon-Wiener diversity index (H ), Simpson diversity index (D) and Pielou evenness index (E) were selected as four commonly used indices to measure species diversity of communities. Basal area standard deviation (SD ba ), basal area coefficient of variation (CV ba ), DBH Shannon-Wiener index (H ba ), and stand density (DEN) were selected as four indicators to measure the diversity of community structure [22]. The coefficient of variation quantifies the degree of dispersion of DBH. Calculation formulas and results of each index are shown in Table 2.  N is the total number of individuals in a subplot; Ns is the total number of species in a subplot; n i is the individual number of the ith species; n j is the individual number of the jth diameter class.

Biomass and Productivity
Biomass includes aboveground (branches, leaves, fruits) and underground biomass (roots). Since there may be errors in estimating biomass, only the allometric growth equation based on DBH was used in this study. Allometric growth equations were selected as appropriate for the subtropical forests of China [40]. Biomass is reported for 2017, and the change in biomass accumulation from 2007 to 2017 was used to measure productivity [41]. The calculation is as follows: Here, BIO is biomass, PRO is productivity, DBH is diameter at breast height (cm), a and b are regression coefficients, t is the time interval (2007 to 2017), Biomass2007 and Biomass2017 are the total biomass of all individuals in the sample at the time of investigation in the respective year. Considering that different initial amounts of biomass within the quadrat may in turn affect the change of biomass within the quadrat, this study used the relative change of biomass to represent the primary productivity of the community (referred to as "productivity" in this paper). To ensure the normality of data, biomass and productivity were log transformed before further analysis.

Statistical Analysis
Before analysis, the explanatory variables of topographic factors, soil factors, species diversity and structural diversity were standardized by log function normalization method. Correlation analysis was conducted using the transformed response and explanatory variables. Then, principal components analysis (PCA) was used to analyze topographic factors, soil factors, species diversity and structural diversity. A unitary regression model was used to examine the relationships between tree diversity, biomass and productivity. Then, Pearson correlation analysis was used to explain the correlation between species diversity, structural diversity, productivity, biomass, soil nutrients and topographic factors. We randomly selected 30 samples to calculate the p-value and repeat this process 1000 times (bootstrapping). Then, we reported the average p-values in our manuscript. The statistical significance level was set as p < 0.05.
To compare the relative effects of biotic and abiotic factors on productivity and to determine the relative contributions of key factors and community dynamics, we used the structural equation model (SEM). We also used the SEM to investigate their direct and indirect effects on biomass and productivity. The guaranteed loading values of the first two PCA axes were chosen as topographic factors, soil factors, species diversity and structural diversity. SEM was performed in Amos 26 and was evaluated for goodness-of-fit using Bentler's Comparative Fit Index (CFI > 0.90), standardized root mean square residual (RMSEA ≤ 0.05), and a Chi-square (χ2) test statistic (p > 0.05) [17,42]. Statistical analysis and plotting were completed in R software version 4.12.

Changes in Community Composition, Structure and Biomass
The number of genera and species in the primary forest decreased by 5 genera and 7 species, and the number of individual trees decreased by 2680, including 1978 new trees and 4658 dead trees from 2007 to 2017. Overall, the number of species and individuals in the whole forest decreased gradually, but the average AGB increased from 53.19 Mg ha −1 to 70.47 Mg ha −1 , and the 10-year average productivity was 1.76 Mg ha −1 yr −1 (Table 3).

The Relationship between Biodiversity and Biomass, Productivity
PCA results (Table S1) show that the first two axes of topographic factors explain 61.47% of total variance, the first two axes of soil factors explain 93.01% of total variance, the first two axes of species diversity explain 96.88% of total variance, and the first two axes of structural diversity explain 76.28% of total variance. The first axis can represent the changes of soil factors and species diversity, and the first and second axes represent the changes of topographic and structural diversity.

Correlation of Topography and Soil Factors with Biodiversity, Biomass, and Productivity
Pearson correlation analysis (Figure 3) showed that productivity was significantly (b) Relationship between loading on the first axis of structural diversity and plant biomass and productivity. (c) Relationship between loading on the second axis of structural diversity and plant biomass and productivity. The grey shadow represent 95% confidence intervals.

Correlation of Topography and Soil Factors with Biodiversity, Biomass, and Productivity
Pearson correlation analysis (Figure 3) showed that productivity was significantly positively correlated with TK (p < 0.001), and significantly negatively correlated with pH, SOM, TN, AN, AP (0.05 < p < 0.01), and there was no significant correlation with TP, AK and topographic factors. Biomass was significantly positively correlated with soil factors except TK (p < 0.001), and significantly negatively correlated with ASP and ROR, and showed no significant correlation with ELE, SOL and CON. Shannon-wiener index of DBH was significantly positively correlated with ELE (p < 0.001) and significantly negatively correlated with CON, and not significantly correlated with other factors. Stand density was significantly positively correlated with soil factors except TK (p < 0.001), and significantly positively correlated with ASP and ROR, and not significantly correlated with ELE or SLO (p > 0.05).
The basal area standard deviation was significantly positively correlated with pH, TP, TK, AP and AK (p < 0.001), and significantly negatively correlated with SOM, TN and AN (p < 0.01); the basal area coefficient of variation was significantly positively correlated with pH, TP and AK (p < 0.001), and significantly negatively correlated with ASP, SOM, TN, TK, AN and AP (p < 0.01). Neither basal area standard deviation nor basal area coefficient of variation were significantly correlated with topographic factors (except ASP). Shannonwiener index of DBH was significantly positively correlated with ELE (p < 0.001) and significantly negatively correlated with CON, and not significantly correlated with other factors. Stand density was significantly positively correlated with soil factors except TK (p < 0.001), and significantly positively correlated with ASP and ROR, and not significantly correlated with ELE or SLO (p > 0.05).

Dynamic Changes of Forest Biomass and Productivity
Over the 2007-2017 study period, the number of genera and species of native forest woody plants decreased, but the biomass increased. The results showed that the average AGB (70.47 Mg ha −1 ) in 2017 was 24.52% higher than that in 2007 (53. 19 Mg ha −1 ), indicating that the forest is still growing and is rich in species composition and stable in structure [43]. However, its AGB was lower than that of the main forests surrounding the karst mountains in the southwest China (122.95-136.4 Mg ha −1 ) [36]. This shows that the primary forest is still in the period of community succession and growth, and that the biomass and productivity are gradually increasing [44]; in addition, different study areas and scales can also cause differences. We also found that the karst forest had lower AGB than typical forests in non-karst regions (152.88 Mg ha −1 ) in the same climate zone [45], which may be due to the shallow soil layer and poor water retention capacity of the karst environment, which restricts the growth of plants and accumulation of biomass [46]. This pattern was attributed to the number of individuals and different distributions of the differ- Figure 4. The impacts of topography, soil nutrients, species diversity and structural diversity on biomass and productivity. The SEM considered all plausible pathways, to increase the degrees of freedom. The numbers on the arrows indicate standardized path coefficients (* p < 0.05, ** p < 0.01, *** p < 0.001). Black arrows represent significant influence and dashed arrows represent non-significant influence. Arrow width is proportional to the strength of the path coefficients. R 2 is the proportion of variance explained.

Dynamic Changes of Forest Biomass and Productivity
Over the 2007-2017 study period, the number of genera and species of native forest woody plants decreased, but the biomass increased. The results showed that the average AGB (70.47 Mg ha −1 ) in 2017 was 24.52% higher than that in 2007 (53. 19 Mg ha −1 ), indicating that the forest is still growing and is rich in species composition and stable in structure [43]. However, its AGB was lower than that of the main forests surrounding the karst mountains in the southwest China (122.95-136.4 Mg ha −1 ) [36]. This shows that the primary forest is still in the period of community succession and growth, and that the biomass and productivity are gradually increasing [44]; in addition, different study areas and scales can also cause differences. We also found that the karst forest had lower AGB than typical forests in non-karst regions (152.88 Mg ha −1 ) in the same climate zone [45], which may be due to the shallow soil layer and poor water retention capacity of the karst environment, which restricts the growth of plants and accumulation of biomass [46]. This pattern was attributed to the number of individuals and different distributions of the different species [47]. Terrain with good growing conditions can capture enough light to help promote growth. Trees in valleys are adversely influenced by harsh conditions such as high humidity and high canopy density, leading to reduced biomass accumulation in valleys than on ridges [48].

The Influence of Species Diversity and Structural Diversity on Community Dynamic
The results of unary linear regression (Figure 2a) showed that species diversity was significantly positively correlated with productivity, consistent with previous studies [15,49]. Again, this suggests that the Karst primary forest site is in the stage of succession and growth, and that the forest structure is relatively stable and the resources can be fully utilized. It is also evident that when the individual species distribution is more uneven, productivity is higher [44]. However, species diversity was negatively correlated with biomass, which is consistent with Zhu, J. et al. [50], and contrary to Wu, C. et al. [15]. The total biomass of the community mainly depends on the species with high carbon storage. When the species with high carbon storage compete or die, the biomass accumulation will be reduced [14]. The complex karst topography and landforms lead to spatially heterogeneous elevations, slopes and niche types, which encourage a faster regeneration rate of some species. This influences the growth of species with high storage, resulting in negative correlation between species diversity and biomass.
The community structure through light complementary utilization to promote productivity, and reflects individuals the complexity of spatial structure, and its diversity reflects the resource allocation of constituent species [51]. The results showed that the structural diversity (PCA1) was negatively correlated with productivity. Productivity was positively correlated with SD ba and CV ba (Figure 2b, Table S1). This indicates that when individuals were more evenly distributed across the full range of diameters, the productivity level was higher [52], and individuals in different diameter classes, different forest layers and different habitats could promote complexity in the canopy structure [51]. Meanwhile, the structural diversity (PCA2) was positively correlated with productivity. Productivity was negatively correlated with H ba and DEN (Figure 2c, Table S1). When the karst forests are influenced by external disturbance, the spatial structure of the community will change, and the impact of community structure (DBH Shannon index) and stand density on productivity will be weakened by scarce resources and decreased absorption efficiency [10]. Likewise, communities with high structural heterogeneity can have reduced productivity [14]. Due to the influences of spatial structure, climate and geographical conditions on karst forest, some species grow slowly: in particular, young plants are subjected to extremely intense pressure of intraspecific competition, with limited resources (such as light and space), leading to high sapling mortality. Together these factors explain the negative correlation between community structure diversity and productivity.
The results showed that the structural diversity axes (PCA1, PCA2) were all negatively correlated with biomass (Figure 2b,c). With increasing community age, community structure complexity and canopy density, some positive species were at a disadvantage in the competition for environmental resources and may not have survived. These factors explain the rapid regeneration rate of the community [18], as well as the higher SD ba index, the greater dispersion, and the increase of stand density, which resulted in fierce resource competition and the decline in biomass.

The Relative Influences of Biotic and Abiotic Factors on Community Dynamics
Biotic and abiotic factors were both closely related to the structure and function of forest communities. Changes in topographic factors and soil nutrient determined species community distribution, forest carbon pool and forest function. Biodiversity is a major driver of ecosystem function. During community succession, species and environment interact and gradually form a stable community structure.
The results above showed that topography had a direct impact on species diversity, structural diversity, biomass and productivity (Figure 4), while topography had a significant negative correlation with biodiversity and productivity. There is more light in convex high elevation areas (such as topographic high points), which is beneficial to the growth of trees, while the restricted light together with shallow and barren soil in the valleys are not suitable for growth, the number of species decreases and forest structure change with increasing elevation from the depression to the slope top and cause species diversity and complexity of species structure gradually decreased from the depression to the slope top ( Figure 3) [53,54], so topography had a significant negative correlation with species diversity. Topography was positively correlated with biomass, which is consistent with Hui et al. and de Castilho et al. [55], at greater slope angles, improved light conditions are conducive to the growth of species with high carbon storage, promoting the increase of biomass. Topography affected the distribution of different diameter classes of species. Although small-diameter individuals rapidly regenerated and improved community productivity, it mainly depended on large-diameter individuals, the number of large-diameter individuals affected by topographic factors decreased, leaded to decreased productivity and had a negative correlation with topography.
Changes in topography and soil nutrients determine multiple aspects of species distribution, forest carbon pool composition and function. Species in natural communities have evolved over a long period of time, forming a specific community structure of interdependence and interaction between species and the environment. The results of showed that topographic factors had no significant effect on biomass and productivity, which was consistent with the results of Zhu Jie et al.'s study on evergreen broad-leaved forests in the subtropical region [50], and different from the results of Wu Chuping' s study on secondary forest in Dinghai, Zhejiang [15], this indicates that the structure of the community is influenced by region and forest type. Elevation was positively correlated with species diversity, while slope was negatively correlated with species diversity. The possible reason was that the topography of the study area changed from depression to steep slope, and species diversity varies with elevation and slope, leading to multiple relationships between topographic factors and species diversity. The nutrient status of the soil is an indicator to measure the nutrient limitation of the environment on the growth of vegetation, the soil nitrogen and phosphorus required by plant growth mainly come from soil. The low availability of phosphorus and high availability of calcium due to high soil weathering in karst areas of southwest China. In this study, it was found that the community structure was complex in the environment with high content of soil phosphorus and soil available potassium in karst primeval forest. While the environment with high soil organic carbon and soil nitrogen content reduced community structure, which was similar to the results of Zhu Jie et al.'s study [50]. The possible causes: soil nitrogen deposition affects the changes of species diversity and structural diversity, and the more nitrogen in soil, the lower the community species diversity and structural diversity, and the lower the biodiversity will affect the growth of forest productivity. Therefore, soil nutrients have a significant impact on vegetation biomass and productivity in karst forests.
There was a direct positive correlation between species diversity and productivity, and a direct negative correlation between species diversity and biomass, but these relationships only explained 13% of variance. Previous studies have confirmed that species richness has a positive impact on forest productivity [56,57]. Other studies have also found that species richness is positively correlated with biomass and has a significant quadratic influence on productivity [15]. It is sometimes believed that species diversity is one of the determinants of ecosystem function [9]: the higher the species richness, the stronger the aggregation traits of dominant species, which can promote productivity. Studies have also shown that biomass can instead decrease at higher functional diversity [19], because when the number of species in the community is small, the dominant species are the pioneer species with high productivity and low carbon reserves. With the increase of species richness, the number of species with high carbon storage increases in the community. Owing to the particular characteristics of the karst sample, most of the species with high carbon storage were located on the slope or on a topographic high, leading to direct and significant positive correlation between biomass and productivity ( Figure 4). There was a direct and significant negative correlation between structural diversity and biomass. Research has shown that stand density can alter biomass accumulation by changing the interaction between species [10]. Species competition increases mortality and changes seedling stand density and spatial pattern [58]. Furthermore, the amount of biomass is mainly controlled by some minority species or individuals in the community, and the sizes of different species and trees have different influences on biomass [42]-in particular, large trees have a higher capacity for photosynthesis and growth and therefore a more direct influence on biomass.

Conclusions
Our study highlights how multiple abiotic and biotic factors interact to influence the biomass and productivity dynamics in a 2 ha dynamic forest plot in Southwest China. We found that the productivity was 1.76 Mg ha −1 yr −1 between 2007 and 2017. Both species diversity and structural diversity were influenced by topographic and soil factors, but their influences on species diversity were stronger. There was no significant difference between species diversity and structural diversity in explaining productivity, while structural diversity was better than species diversity in explaining biomass accumulation. Both biotic factors and abiotic factors can directly explain biomass and productivity, but biotic factors are better than abiotic factors in explaining biomass and productivity. This study explored the dynamics of community biomass and productivity in subtropical forests, and identified key influencing factors. There are two major limitations in this study that could be address in future research. First, the area of the study sample plot is only 2 ha, and the results obtained are of reference value, but not completely representative of other forests. Second, the depression in the study sample was partially disturbed by natural disasters, which caused a little damage to the species in the depression, but the overall impact was not significant, and this is why the topographic factors have various influence relationships on species diversity and the topographic factors have no effect on species diversity. Over time, the forest restoration process provides an ideal place for subsequent sustainable development and research.

Implication
Based on the results presented above, it is reasonable to conclude that forest productivity is closely related to tree species structure, stand density, and species diversity in subtropical karst forest, and therefore, measures such as forest tending should be emphasized to increase forest productivity, so as to achieve sustainable management of forests. Therefore, the combined influences of biotic and abiotic factors should be considered in future studies on biodiversity and ecosystem function, providing a scientific basis for understanding the services and functions of karst ecosystems. These results provide novel insights into the driving mechanism of forest productivity, which is beneficial to achievie effective protection and the sustainable utilization of forest resources, maintain the stability of the forest ecosystem, as well as help policy makers in planning sustainable forest management are of great significance.
Supplementary Materials: The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/su14137644/s1, Table S1: Mean and contribution rate of topography factors, soil factors, species diversity factors and structural diversity factors variables of each 10 m × 10 m quadrat and Correlation between them with the PCA axes.