Composition of Species and Spatial Patterns of Phyllostachys edulis Mixed Forests in a Succession after Enclosure

: This paper investigates the changes in species diversity and the spatial pattern of the population of a moso bamboo ( Phyllostachys edulis ) forest community in the Wuyi Mountain Nature Reserve in


Introduction
Community succession is a natural phenomenon that results from the continuous replacement of community composition types and changes in the habitat caused by the transformation of communities in the environment [1][2][3][4].Community succession in enclosed forests is called closed succession.The results of this study showed that the vegetation cover of degraded forest land can be improved through enclosure.Numerous studies have shown that enclosure can increase the diversity of vegetation and biomass of the community [5,6] and can also improve the soil moisture [7].Studies on species diversity and spatial patterns provide important information to understand the mechanisms used by species to coexist [8,9].The study of species diversity and spatial patterns has been favored by researchers in China and throughout the world [10][11][12][13].Changes in the diversity of community species are a central objective of community succession [14,15].Such changes reflect the diversity and spatial heterogeneity of the community structure [16,17].The pattern of spatial distribution can be an important indicator of community health and stability [18].Many researchers in China and throughout the world have conducted extensive studies on the composition of species and the spatial pattern of communities at different phases of succession.Shan [19] and other researchers have studied the species diversity of grassland communities at different stages of succession.Isermann [20] studied the changes in species diversity during the succession of coastal dune communities, while Liu et al. [21] studied the spatial pattern of the distribution of Artemisia ordosica Krasch.populations in response to the cordoning off of their site.In addition, Meiners [22] studied the population distribution dynamics of plant species during different succession processes.However, little research has been reported on the changes in species diversity and spatial patterns of the plant populations of Phyllostachys edulis forest communities at different stages of succession.
The quality, production value, and economic value of Phyllostachys edulis in the Wuyi Mountain Nature Reserve, Fujian Province, China, are among the best in the country.In recent years, Wuyi Mountain City has adopted a policy of enclosing the forest, which encourages the Phyllostachys edulis forest community to use its ability for self-renewal to restore and improve the forest environment.In this study, we selected mixed Phyllostachys edulis forests on Wuyi Mountain with different years of enclosure as the research object and comprehensively investigated and studied the mixed Phyllostachys edulis forest communities with different years of enclosure using established standard plots.In this study, Phyllostachys edulis mixed forests in the Wuyi Mountain Nature Reserve in Fujian Province, China were used the experimental object, their enclosure times were 3, 5, 10, 15, 40 years and were compared with the conventional operation of the Phyllostachys edulis forest.The results of this study provide a reference for the successional progress of enclosed mangosteen forests and their management.

Study Site
Wuyi Mountain National Park is located in the northern part of Fujian Province (117 • 24 13 -117 • 59 19 , 27 • 31 20 -27 • 55 49 ), which is situated in the humid zone of the central subtropical monsoon, with an annual average temperature of 17.7 • C, precipitation of 1800 mm, an annual average duration of sunshine of 1910.2 h, and an average of 79% relative humidity.As the altitude of these mountains increases, there are four types of vegetation zones with different community characteristics, namely, evergreen broad-leaved forests, mixed coniferous and broad-leaved forests, dwarf forests, and mesic meadows, with clear demarcation lines and obvious plant vertical zones from the foot of the mountain upwards [23].

Study Design
After a comprehensive survey of the Wuyishan National Forest Park in 2021, mixed Phyllostachys edulis forests with 3, 5, 10, 15, and 40 years of enclosure in the Aotou, Tongmu, and Da'an villages in Wuyishan City were selected as the survey object, and conventionally operated Phyllostachys edulis forests were used as the control (0 years).The six communities were sequentially numbered as Type A (control), Type B (3 years), Type C (5 years), Type D (10 years), Type E (15 years), and Type F (40 years).The basic conditions of the sample plots and stands in the six different treatments are shown in Tables 1 and 2. Three sample plots were established in each of the communities that had been enclosed in different years.This resulted in 18 plots in total that were 25.8 m × 25.8 m.The enclosed mangosteen (Garcinia mangostana Linn) forests were not disturbed by anthropogenic activities, such as digging for bamboo shoots and cleaning old bamboo, and were in a completely enclosed state.Regularly operated Phyllostachys edulis forests are normally utilized for a variety of business activities.Sample plots selected during the survey were fully representative, did not cross streams, roads, or logged survey lines, and were far away from the forest edge.The sample plots were established in places where the tree species and stand densities were evenly distributed.Three different replicate sample plots were established in the same stand and did not cross the stand.The different replicate sample plots were spaced not less than 25 m and not more than 50 m apart (Figure 1).
x FOR PEER REVIEW 3 of 21 enclosed mangosteen (Garcinia mangostana Linn) forests were not disturbed by anthropogenic activities, such as digging for bamboo shoots and cleaning old bamboo, and were in a completely enclosed state.Regularly operated Phyllostachys edulis forests are normally utilized for a variety of business activities.Sample plots selected during the survey were fully representative, did not cross streams, roads, or logged survey lines, and were far away from the forest edge.The sample plots were established in places where the tree species and stand densities were evenly distributed.Three different replicate sample plots were established in the same stand and did not cross the stand.The different replicate sample plots were spaced not less than 25 m and not more than 50 m apart (Figure 1).

Community Investigation Methods
For the adult trees, 25.8 m × 25.8 m sample plots were considered to be the basic unit.A per-tree survey for all plants with a diameter at breast height (DBH) > 3 cm in the sample plot was conducted.The species name, DBH, tree height (H), and crown width of each individual were recorded.Simultaneously, the lower-left corner of each sample plot was considered to be the origin when the coordinates of each individual (x-and y-axis) were recorded.For young trees, 25.8 m × 25.8 m was considered to be the basic unit, and per-tree surveys for all plants >1 m high and which had a DBH < 3 cm in the sample area were recorded along with the species name, DBH, H, and crown width of each individual plant.For the shrub and herbaceous layers, 5 m × 5 m was the basic unit for shrubs and 1 m × 1 m was the basic unit for the herbaceous plants.The species name, number of individuals, coverage, and average height of the shrubs and herbaceous plants in each plot were recorded.The environmental characteristics of each site, including latitude, longitude, altitude, geomorphology, soil type, slope, and slope position, among others, were also recorded.

Calculation of the Species Diversity Index (SDI)
The relative density, relative frequency, and relative significance (relative cover) of each species were calculated based on the information of the sample plots and the content of the survey, and their important value (IV) was calculated.
The significant values for mature and young trees were shown as follows: The importance values of the shrub and herb layers were as follows: α-diversity, also known as the diversity within habitats, focuses on the number of species in a homogeneous environment within a local area.Based on the principle of multi-index selection, richness, diversity, and evenness indices were chosen for calculation, which can reflect the changes in α-diversity in a more comprehensive way.
(1) Richness index, described as the Margalef index: (2) Diversity indices: Simpson's index: Shannon-Wiener index: Forests 2023, 14, 2325 5 of 21 (3) Uniformity index: Pielou index: Ma is the Margalef's index; H is the Simpson's index; H is the Shannon-Wiener's index; J SW is the Pielou's index; p i is the proportion of the number of individuals of species i; n i , is the total n of individuals of all the species, i.e., p i = n i /n; i = 1, 2, 3 S; S is the average number of species in the sample plot sample; and B is the total biomass of the community [24].

Point Pattern Analysis
The point pattern analysis can be divided into univariate pattern analysis, bivariate pattern analysis, and multivariate pattern analysis, and the distribution patterns of communities and populations, the relationships between and among populations, and the reflection of their ecological processes at different scales can be explored using different null hypothesis models.The L (r) function was derived from Ripley's K function and can explain the actual spatial pattern more directly and simply than the K (r) function.Thus, the L (r) function was used to describe the spatial patterns of individuals after different periods of enclosure at various spatial scales in the Phyllostachys edulis forests.In analyzing the spatial associations among different species, the K 12 (r) function was extrapolated on the basis of Ripley's K function.Similarly, the L 12 (r) function was obtained based on the variance correction and linearization of the K 12 (r) function.The principle and calculation procedure are as described by Wiegand and Moloney [25,26].

Zero Model Selection
In this study, the univariate L (r) function was chosen to be the complete spatial randomness (CSR) model [27,28], and 100 stochastic simulations were conducted using the Monte Carlo method to test the deviation of the null model from the observed values and to construct the upper and lower two packet trajectory lines by applying the simulated maxima and minima, referred to as T. The intervals of the distributions with a confidence level of 95% were then intercepted.If the actual distributions of the L (r) values and L 12 (r) fell within the envelope, the plants were randomly distributed or the two enclosure years were independent of each other and lacked a spatial correlation.If they were above the envelope, they were clustered or there was a spatial positive correlation between the two years of enclosure.If they were below the envelope, they were uniformly distributed or there was a spatial negative correlation between the two enclosures.

Statistical Analysis
Microsoft Excel 2013 (Redmond, WA, USA) was used to integrate data, and SPSS 19.0 (IBM, Inc., Armonk, NY, USA) was used for analysis of variance (ANOVA).The Kruskal-Wallis test was used for the indicators within the moso bamboo forest communities under different periods of enclosure, while significant difference was found using the Dunn test.A point pattern analysis was performed using Programita (2014 version) software.The pairwise correlation function g (r) was chosen; the spatial scale of the sample plots was 0-40 m; the ring width was established as 3 m; the analytical step size was 0.2 m; and the 95% confidence intervals were determined through 199 simulations using the Monte Carlo method.

Plant Species and Importance Values in Communities following Different Periods of Enclosure
The primary adult trees (importance value ≥ 1) were members of 30 genera in 23 families, which resulted in a total of 41 species in the six different treatment samples (Table 3).Among them, plants in the beech family (Fagaceae) had the highest number of species, which totaled ten, followed by members of the Pentaphylacaceae family and Symplocaceae, with five and four species, respectively.In addition, there were 18 families with only one species, which comprised 78.26% of the total number of families.In Type B, the ranking of plant IVs in order was: Symplocos congesta Benth., Manglietia fordiana Oliv., Castanopsis eyrei (Champ.ex Benth.)Tutch., and Castanopsis faberi Hance.Symplocos congesta Benth.( 18.869%) and Castanopsis faberi Hance (13.310%) were the most predominant of the adult trees.In Type F, the top five adult trees ranked according to the IVs of the plant species were Castanopsis faberi Hance, Pinus massoniana Lamb., Machilus thunbergii Sieb.et Zucc., Alniphyllum fortunei (Hemsl.)Makino, and Eurya nitida Korthals in descending order.Castanopsis faberi Hance (14.431%) and Pinus massoniana Lamb.(13.511%) were significantly more important than the other adult tree plants.The dominant species of adult trees after different periods of enclosure varied.With the extension of the enclosure of the forests, the IVs of Machilus thunbergii Sieb.et Zucc., Castanopsis faberi Hance, and Pinus massoniana Lamb.increased continuously; they gradually became the dominant species of adult trees.The IVs of Cunninghamia lanceolata (Lamb.)Hook.showed a tendency to decrease, increase, and then decrease.It became the dominant species of adult tree in the community following 15 years of enclosure.In the six different treatment plots, the predominant young trees (importance value ≥ 1) were members of 32 genera in 24 families, which totaled 47 species (Table 4).The highest number of species was found in the family Fagaceae, with seven species, followed by five and four species in the families Pentaphylacaceae, Symplocaceae, and Lauraceae, and 16 families with only one species, which comprised 66.67% of the total number of families.Juvenile plants that occurred within all of the communities of the different treatment samples were Machilus thunbergii Sieb.et Zucc., Schima superba Gardn.etChamp., and Castanopsis faberi Hance.In Type B, the top five juvenile tree plants in terms of the plant species IV were Liquidambar formosana Hance, Machilus thunbergii Sieb.et Zucc., Photinia bodinieri Lévl., Camellia oleifera Abel., and Schima superba Gardn.etChamp. in order, with Liquidambar formosana Hance (39.069%) and Machilus thunbergii Sieb.et Zucc.(15.851%) comprising the predominant juvenile trees.In Type C, Symplocos pseudobarberina Gontsch., Castanopsis eyrei (Champ.ex Benth.)Tutch., Photinia prunifolia (Hook.et Arn.) Lindl., Schima superba Gardn.et Champ., and Machilus thunbergii Sieb.et Zucc.were the predominant species (in order).Symplocos pseudobarberina Gontsch.(20.767%) and Castanopsis eyrei (Champ.ex Benth.)Tutch.(17.419%) were significantly more important than the other juvenile trees.In Type D, the top five juvenile trees in terms of the plant species IVs were Schima superba Gardn.etChamp., Machilus chrysotricha H. W. Li, Photinia bodinieri Lévl., Adinandra millettii (Hook.et Arn.) Benth.et Hook.f. ex Hance, and Eurya nitida Korthals, with populations of Schima superba Gardn.etChamp.(13.449%) and Machilus chrysotricha H. W. Li (11.367%) that were significantly higher than those of the other juvenile trees.In Type E, the plants ranked according to their IVs were Adinandra glischroloma var.macrosepala (Metcalf) Kobuski, Castanopsis faberi Hance, Cyclobalanopsis glauca Thunb., Eurya muricata Dunn, and Cunninghamia lanceolata (Lamb.)Hook., with the IVs of Adinandra glischroloma var.macrosepala (Metcalf) Kobuski (16.948%) and Castanopsis faberi Hance (13.640%) being significantly higher than those of the other juvenile trees.In Type F, the top five juvenile trees in terms of the IVs of the plant species were Machilus chrysotricha H. W. Li, Symplocos stellaris Brand, Adinandra millettii (Hook.et Arn.) Benth.et Hook.f. ex Hance, Diospyros morrisiana Hance, and Engelhardia roxburghiana Wall.(in order).The dominance of Machilus chrysotricha H. W. Li (10.533%) and Symplocos stellaris Brand (4.695%) was readily apparent among the juvenile trees.With the increase in the length of time that the forests were enclosed, among the young trees, the IVs of Machilus thunbergii Sieb.et Zucc.decreased, and its dominance gradually weakened, while the IVs of Machilus chrysotricha H. W. Li gradually increased, and it became the dominant plant among the young trees in the community when the length of enclosure was 40 years.
Plants in the main shrub layer (importance value ≥ 1) were members of 32 genera in 23 families and totaled 45 species in the six different treatment plots (Table 5).Among them, plants of the Rosaceae, Theaceae, and Pentaphylacaceae families had the highest number of species, four each in total, while the Primulaceae, Ericaceae, Cyperaceae, and Symplocaceae families had the second highest number of species, each with three.A total of 13 families had only one species, which accounted for 56.52% of the total number of families.Shrub plants that occurred in all the communities in the different treatment plots were Eurya rubiginosa var.Attenuata H.T.Chang.In Type A, the top five shrub-layer plants in terms of the plant species IVs were Tashiroea sinensis Diels, Maesa japonica (Thunb.)Moritzi.ex Zoll., Coptosapelta diffusa (Champ.ex Benth.)Van    Plants in the primary herbaceous layer (IV ≥ 1) were members of 27 families and 38 genera and totaled 46 species in the six different treatment samples (Table 6).The Orchidaceae family had the highest number of species, totaling six, followed by the Rosaceae and Cyperaceae families, with four and three species, respectively.A total of 13 families had only one species, which comprised 48.15% of the total number of families.The herbaceous plants that occurred within all the communities of the different treatment sample plots were Carex cruciata Wahlenb.and Woodwardia japonica (L.F.) Sm.In Type A, the top five herbaceous-layer plants ranked by the IVs of the plant species were Trachelospermum jasmi-noides (Lindl.)Lem., Woodwardia japonica (L.F.) Sm., Carex cruciata Wahlenb., Smilax china L., and Sarcandra glabra (Thunb.)Nakai (in order).Trachelospermum jasminoides Lem.(20.765%) and Woodwardia japonica (L.F.) Sm. (16.372%) in the herbaceous layer were significantly dominant.In Type B, the dominant plants, in order, were Dianella ensifolia Redouté, Carex scaposa C. B. Clarke, Viola diffusa Ging., Woodwardia japonica (L.F.) Sm., and Carex cruciata Wahlenb., among which Dianella ensifolia Redouté (21.136%) and Carex scaposa C. B. Clarke (13.537%) had significantly higher IVs than the other plants in the herbaceous layer.In Type C, the top five herbaceous plants in terms of the importance of plant species were Ardisia japonica (Thunberg) Blume, Carex cruciata Wahlenb., Liriope spicata (Thunb.)Lour., Huperzia serrata (Thunb.ex Murray) Trevis., and Dryopteris fuscipes C. Chr.(in order).Ardisia japonica Blume (51.810%) and Carex cruciata Wahlenb.(10.070%) had significantly higher IVs than the other plants in the herbaceous layer.In Type D, the most important plants were Carex cruciata Wahlenb., Woodwardia japonica (L.F.) Sm., Diplopterygium glaucum (Thunberg ex Houttuyn) Nakai, Alpinia japonica (Thunb.)Miq., and Carex maubertiana Boott (in order).Carex cruciata Wahlenb.(28.048%) and Woodwardia japonica (L.F.) Sm. (16.599%) were the dominant herbaceous plants.In Type E, the top five plants in the herbaceous layer ranked in the order of importance of the plant species were Trachelospermum jasminoides (Lindl.)Lem., Woodwardia japonica (L.F.) Sm., Diplopterygium glaucum (Thunberg ex Houttuyn) Nakai, Carex cruciata Wahlenb., and Sarcandra glabra (Thunb.)Nakai.Trachelospermum jasminoides (Lindl.)Lem.(30.140%) and Woodwardia japonica (L.F.) Sm. (18.558%) were significantly more important than the other herbaceous-layer plants.In Type F, the most dominant plants were Carex cruciata Wahlenb., Diplopterygium glaucum (Thunberg ex Houttuyn) Nakai, Woodwardia japonica (L.F.) Sm., Osmunda japonica Thunb., and Dianella ensifolia (L.) Redouté.Carex cruciata Wahlenb.(53.395%) and Diplopterygium glaucum (Thunberg ex Houttuyn) Nakai (22.166%) were significantly predominant in the herbaceous layer.The IVs of Carex cruciata Wahlenb.generally increased with the increase in the number of years of enclosure, and it became the dominant species in the herbaceous layer of the understory after 40 years of enclosure.

Species Diversity Indices within the Communities following Different Periods of Enclosure
In the six treatment plots, the Margalef, Simpson, and Shannon-Wiener indices of the adult trees gradually increased with the number of years that the forest had been enclosed (Tables 7 and 8).The Margalef index of Type B was significantly different from those of both Types D and F (p < 0.05), and the Margalef index of Type B was 64% and 81% higher than those of Types D and F, respectively.The Simpson index of Type B was significantly different from that of Type F (p < 0.05), while the Simpson index of Type B was 32% higher than that of Type F. The Shannon-Wiener index of Type F was significantly different (p < 0.05) from those of Types D, C, and B. The Shannon-Wiener index of Type F was 61%, 195%, and 202% higher than those of Types D, C, and B, respectively.There was no significant change in the Pielou index of the six sample plots over different periods with the increase in the number of years of enclosure.The Pielou index of Type F was significantly different from that of Type D (p < 0.05), and the Pielou index of Type F was 11% higher than that of Type D. This indicated that most of the diversity indices of the adult arboreal plants increased with the extension of the enclosure period.The Margalef and Shannon-Wiener indices of the juvenile trees increased with the extension of the periods of enclosure in the six treatment samples (Tables 9 and 10).The Margalef index of Type B was significantly different from those of both Types E and F (p < 0.05), and the Margalef index of Type B was 44% and 64% higher than those of Types E and F, respectively.The Shannon-Wiener index of Type B was significantly different from those of Types E and F (p < 0.05), and the Shannon-Wiener index was 28% and 40% higher than those of Types E and F, respectively.The Simpson index of Type B was significantly different from that of Type F (p < 0.05), and there was no significant difference between the Pielou indices of the six different treatment samples.With the increase in the number of years of enclosure, the Simpson and Pielou indices of the six different treatment samples did not change significantly, indicating that the changes in the diversity index of juvenile arborvitae were less affected by the number of years of enclosure.In the six treatment plots, the Margalef, Shannon-Wiener, and Pielou indices of the shrub layer increased and then decreased with the extension of the enclosure years (Tables 11 and 12).The Margalef index of Type E was significantly different from that of Type A (p < 0.05), and the Margalef index of Type E was 340% higher than that of Type A. The Shannon-Wiener indices of Types F and E were significantly different from that of Type A (p < 0.05), and the Shannon-Wiener indices of Types F and E were 105% and 127% higher than that of Type A, respectively.The Simpson index did not change significantly with the increase in the years of enclosure, and there was no significant difference between the Simpson and Pielou indices in the six sample plots with different treatments (p > 0.05).This indicated that, in the plant diversity index of the shrub layer, the change in enclosure years had a greater effect on the Margalef and Shannon-Wiener indices.The Margalef index and Shannon-Wiener index were significantly affected by the change in enclosure years, but the Simpson index and Pielou index were less affected.In the six treatment plots, the Margalef, Simpson, and Shannon-Wiener indices of the herbaceous layer gradually decreased with the extension of the enclosure years (Tables 13 and 14).There was a significant difference (p < 0.05) between the Margalef index of Types F and A, and the Margalef index of Type F was 83% higher than that of Type A. There were significant differences (p < 0.05) in the Shannon-Wiener index values of Types F and C and Types B and A. In addition, the Shannon-Wiener index of Type F was higher than those of Types C, B, and A by 65%, 68%, and 68%, respectively.The Pielou index did not change significantly with the increase in years of enclosure, and there was no significant difference between the Simpson and Pielou indices of the six sample plots with different treatments (p > 0.05), which indicated that, in the plant diversity indices of the herbaceous layer, the change in enclosure years had no significant effect on the Margalef and Shannon-Wiener indices and less on the Simpson and Pielou indices.

Size of the Bamboo Forest Community Structure and Spatial Distribution of Individuals in Each Period of Enclosure
A total of 421 individual tree-layer plants with DBH ≥ 1 cm were recorded in the study sample (Figure 2).The largest DBH was 69.0 cm, and the average DBH of the plant community was 15.7 cm.A total of 208, 97, 50, 29, 24, 12, and one plants had DBHs of 1-10 cm, 10-20 cm, 20-30 cm, 30-40 cm, 40-50 cm, 50-60 cm, and 60-70 cm, respectively, which accounted for 49.4%, 23.0%, 11.9%, 6.9%, 5.7%, 2.9%, and 0.2% of the total number of individuals in the tree layer, respectively.The number of individuals in the tree layer with a DBH of 1-10 cm comprised the largest percentage of the total number of individuals, and the number of individuals with a diameter in the range of 1-70 cm decreased with the increase in the diameter level.In terms of the distribution of plant size, there were fewer old trees in the tree layer in the study site, and there was a higher proportion of small plants.The distribution of tree-layer plants showed an inverted "J" shape, and the population belonged to the growth type.This indicated that the tree-layer plant population in the Phyllostachys edulis forest community had a high reserve of forest renewal, and the rate of renewal was rapid and very effective.This indicates that the plant population in the tree layer of the Phyllostachys edulis forest community has a high reserve of forest trees, the forest trees renewed quickly, the state of renewal is very good, and the community is in a relatively stable state.The years of enclosure had some effect on the spatial distribution of the six different individuals in the same site (Figure 3).When the duration of enclosure was 3 years, the tree layer was densely distributed in the upper right part of the sample site and sparsely distributed in the middle and lower parts of the sample site.When the period of enclosure was 15 years, the tree layer was more densely distributed in the lower left part of the sample plot.When the enclosure period was 40 years, the tree-layer plants were densely distributed in the upper left, lower left, and upper right of the sample plots, and they were sparsely distributed in the lower middle and lower right of the sample plots.Overall, the layers of trees within the Phyllostachys edulis forest community varied in the sparseness of their individual distribution within the sample plot.2).The largest DBH was 69.0 cm, and the average DBH of the plant community was 15.7 cm.A total of 208, 97, 50, 29, 24, 12, and one plants had DBHs of 1-10 cm, 10-20 cm, 20-30 cm, 30-40 cm, 40-50 cm, 50-60 cm, and 60-70 cm, respectively, which accounted for 49.4%, 23.0%, 11.9%, 6.9%, 5.7%, 2.9%, and 0.2% of the total number of individuals in the tree layer, respectively.The number of individuals in the tree layer with a DBH of 1-10 cm comprised the largest percentage of the total number of individuals, and the number of individuals with a diameter in the range of 1-70 cm decreased with the increase in the diameter level.In terms of the distribution of plant size, there were fewer old trees in the tree layer in the study site, and there was a higher proportion of small plants.The distribution of tree-layer plants showed an inverted "J" shape, and the population belonged to the growth type.This indicated that the tree-layer plant population in the Phyllostachys edulis forest community had a high reserve of forest renewal, and the rate of renewal was rapid and very effective.This indicates that the plant population in the tree layer of the Phyllostachys edulis forest community has a high reserve of forest trees, the forest trees renewed quickly, the state of renewal is very good, and the community is in a relatively stable state.The years of enclosure had some effect on the spatial distribution of the six different individuals in the same site (Figure 3).When the duration of enclosure was 3 years, the tree layer was densely distributed in the upper right part of the sample site and sparsely distributed in the middle and lower parts of the sample site.When the period of enclosure was 15 years, the tree layer was more densely distributed in the lower left part of the sample plot.When the enclosure period was 40 years, the treelayer plants were densely distributed in the upper left, lower left, and upper right of the sample plots, and they were sparsely distributed in the lower middle and lower right of the sample plots.Overall, the layers of trees within the Phyllostachys edulis forest community varied in the sparseness of their individual distribution within the sample plot.

Spatial Distribution Pattern of Moso Bamboo Forest Communities with Different Enclosure Periods
The year of enclosure had an effect on the pattern of spatial distribution of the populations in the six different sites (Figure 4).When the period of enclosure was 3 years, the estimates of arborvitae plants at the 0-3.5 cm scale were above the confidence interval, and the plants of the community showed a clustered distribution.When the period of enclosure was 5 years, the estimates of the arborvitae at the 0-4 cm scale were above the confidence interval, and the plants in the community showed a clustered distribution.

Spatial Distribution Pattern of Moso Bamboo Forest Communities with Different Enclosure Periods
The year of enclosure had an effect on the pattern of spatial distribution of the populations in the six different sites (Figure 4).When the period of enclosure was 3 years, the estimates of arborvitae plants at the 0-3.5 cm scale were above the confidence interval, and the plants of the community showed a clustered distribution.When the period of enclosure was 5 years, the estimates of the arborvitae at the 0-4 cm scale were above the confidence interval, and the plants in the community showed a clustered distribution.When the enclosure period was 10 years, the estimates of the arborvitae at the 0-7 cm scale were below the confidence interval, and the plants in the community were uniformly distributed.When the enclosure period was 15 years, the estimates of the arborvitae at the 0-7 cm scale were above the confidence interval, and the plants in the community showed a clustered distribution.At 40 years of enclosure, the estimates of the arborvitae at the 0-7 cm scale were above the confidence interval, and the plants in the community showed a clustered distribution.This indicates that the pattern of spatial distribution of the arborvitae in the Phyllostachys edulis forest community was generally clustered with the extension of the enclosure period.When the enclosure period was 10 years, the estimates of the arborvitae at the 0-7 cm scale were below the confidence interval, and the plants in the community were uniformly distributed.When the enclosure period was 15 years, the estimates of the arborvitae at the 0-7 cm scale were above the confidence interval, and the plants in the community showed a clustered distribution.At 40 years of enclosure, the estimates of the arborvitae at the 0-7 cm scale were above the confidence interval, and the plants in the community showed a clustered distribution.This indicates that the pattern of spatial distribution of the arborvitae in the Phyllostachys edulis forest community was generally clustered with the extension of the enclosure period.

Discussion
The relative importance of plants within a community can be expressed by the IVs of plants, and the IV is also one of the important indicators to assess the amount of plant diversity within a population.In this study, with the extension of the period of enclosure, the IVs of Machilus thunbergii Sieb.et Zucc., Castanopsis faberi Hance, and Pinus massoniana Lamb.increased continuously, and they gradually became the dominant species of adult trees.The IVs of Cunninghamia lanceolata (Lamb.)Hook.tended to decrease, increase, and then decrease again, and it became the dominant plant of adult trees within the community when the forests had been enclosed for 15 years.This study showed that Cunninghamia lanceolata (Lamb.)Hook.and Pinus massoniana are the establishment species of subtropical evergreen coniferous forests in China, with a strong adaptive ability, and the vegetation gradually transitioned to mixed coniferous and broadleaf forests during the process of sequestration; they finally developed into a plant community dominated by evergreen broadleaf forests [29].In the process, Machilus thunbergii Sieb.et Zucc.and Castanopsis faberi Hance were the dominant zonal plants, and they were the top species in the succession [30].In this study, the IVs of Eurya rubiginosa var.attenuata H.T.Chang in the shrub layer showed a general trend of increasing and then decreasing, and it became the dominant species in the shrub layer under the Phyllostachys edulis forest after 10 and 15 years after enclosure, while the IVs of Carex cruciata in the herbaceous layer generally increased.It became the dominant species in the herbaceous layer under the Phyllostachys edulis forest after 40 years of enclosure.The results showed that Eurya rubiginosa var.attenuata H.T.Chang has shallow roots and usually grows in the understory and barren slope thickets on mountains and hilltops.It is a common understory shrub-layer plant in subtropical broad-leaved evergreen forests in China [31,32], and thus, it can become the dominant species of shrub layer vegetation after 10 and 15 years after enclosure.Carex cruciata Wahlenb.has a well-developed root system and a strong ability to adapt to the environment, and it can grow tenaciously in different environments; this makes it an important constituent of the herbaceous layer [33].It gradually eliminated the other species and became the absolute dominant plant in the herbaceous layer as the number of years of enclosure increased.
The restoration of species diversity and community structure within forests is an important task for modern researchers [34].In this study, the shrub layer dominance Simpson index was lower in the community after 10 years of enclosure than in the other communities, primarily owing to the presence of the absolutely dominant shrub plant, Eurya rubiginosa var.attenuata H.T.Chang, in the community after 10 years of enclosure.The Shannon-Wiener index of plant diversity in the shrub layer tended to increase and then decrease with the increasing period of enclosure.With the further increase in the period of enclosure, the diversity of shrub layer vegetation tended to decrease, primarily owing to the increase in the density of Phyllostachys edulis and trees.Consequently, there is an increasingly large degree of canopy enclosure; the penetration of sunlight is weakened, and the competition between plants is more intense, which inhibits the growth of shrubs [35,36].The shrub layer uniformity Pielou index tended to increase and then decrease with the period of enclosure, and the shrub layer uniformity Pielou index was lower in the community that had not been subjected to a full year of enclosure.This may be owing to the fact that the species richness, S, of the shrub layer was higher in the zero-year community, which increased the denominator of Jsw = H /InS and led to a decrease in the Pielou index.The vegetation evenness index and the Pielou index were positively correlated with the diversity index and the Shannon-Wiener index and were negatively correlated with species richness, S. The results of other studies were similar to those of this study.
In this study, the Margalef richness index of the herbaceous layer tended to decrease as the period of enclosure increased, and the Simpson dominance index and the Shannon-Wiener diversity index were negatively correlated with the age of enclosure.With the increase in the enclosure time, the stand enclosure gradually increased, which led to a significant decrease in the light available to the herbaceous plants.Simultaneously, the increase in stand density led to a more intensive competitive environment for the herbaceous plants, so that only the better adapted plant species could survive and become the absolute dominant species in the herbaceous layer.This pattern resulted in a gradual decrease in the diversity of the plant species and the trend of the gradual decline of the herbaceous-layer plants [37].
The distribution of plants according to the DBH or the percentage of different DBH sizes in a community is called the diameter structure of the community, and the current status of the population structure and renewal strategy can be reflected by the diameter structure, which can also be used as an effective method to predict the direction of the succession of a population or a community [38,39].In this study, the structure of the adult trees was determined by the size of the tree.For example, the average diameter class of the adult trees tended to increase in parallel with the increase in the period of enclosure, which made them part of the growing population.They showed that the number of young trees in these communities increased with the duration of enclosure; the vegetation regenerated effectively, and the number of plants in the communities remained stable or continued to increase [40].The results of this study showed that the number of young trees in these communities is increasing with the time of enclosure.
The spatial pattern of distribution of the plant populations can be represented by a quantitative expression of the distribution of populations in two dimensions [41], which is closely related to scale and changes with scale [42].In this study, the spatial pattern of distribution of adult arborvitae plants within the community was broadly shown as an aggregated distribution with the extension of the years of enclosure.This is related to the scale of the study, where the spatial distribution pattern of plant populations is primarily related to the biological characteristics of the plants within the populations [43].The pattern of spatial distribution of the plant populations at large scales is affected by different factors, and this pattern of plant populations at different scales is not the same.The research sample sites were selected in a small-scale environment, which resulted in the spatial distribution pattern of adult tree plants being mostly aggregated.With the increase in the time of enclosure, the need for nutrients, light, and other materials of plants in the community gradually increased, and the competition between different individuals was also strengthened.Thus, many individuals of different species of plants died along with those within the same species owing to the intense competition, the phenomenon of selfthinning, and other thinning.This resulted in a substantial reduction in population density, and the intensity of aggregation within the population was also weakened [41].However, this phenomenon differs from the pattern derived from the results of this study.This may be owing to the fact that there were still many Phyllostachys edulis plants in the study area.Thus, the competition within the same species and between different species did not lead to a gradual decrease in the aggregation strength of the adult tree plant populations with the development of closed succession.In contrast, with the development of closed succession, the spatial distribution pattern of adult tree plant populations was mostly aggregated, which led to the conclusion that is inconsistent with the phenomenon of self-thinning and other thinning.

Conclusions
Comprehensive analyses showed that, with the increase in the time of enclosure, the dominance of adult trees in the community increased and gradually formed a symbiotic pattern with Phyllostachys edulis.The spatial pattern of distribution of the populations was generally in a state of aggregation, and it was apparent that the Phyllostachys edulis forest community under prolonged enclosure conditions responded in a positive direction and ultimately formed a more stable Phyllostachys edulis mixed forest community.In order to better understand the effects of Phyllostachys edulis succession on ecosystem functions, it is necessary to study the growth and environmental changes of Phyllostachys edulis.

Figure 1 .
Figure 1.Experimental scheme of the Phyllostachys edulis forest community.

Figure 1 .
Figure 1.Experimental scheme of the Phyllostachys edulis forest community.

Figure 2 .
Figure 2. DBH class distribution of the Phyllostachys edulis forest community.DBH, diameter at breast height.Type B: 3 years; Type C: 5 years; Type D: 10 years; Type E: 15 years; Type F: 40 years; different lowercase letters indicate significant differences at p < 0.05.

Figure 2 .
Figure 2. DBH class distribution of the Phyllostachys edulis forest community.DBH, diameter at breast height.Type B: 3 years; Type C: 5 years; Type D: 10 years; Type E: 15 years; Type F: 40 years; different lowercase letters indicate significant differences at p < 0.05.

Figure 3 .
Figure 3. Spatial distribution of individuals of the Phyllostachys edulis forest community.

Figure 3 .
Figure 3. Spatial distribution of individuals of the Phyllostachys edulis forest community.

Figure 4 .
Figure 4. Single variable point pattern analysis of the Phyllostachys edulis forest community after different periods of enclosure.L ^iso (r): estimated value; L ^trans (r): fitting value; L ^bord (r): upper limit of the confidence interval; Lpois (r): lower limit of the confidence interval.

Figure 4 .
Figure 4. Single variable point pattern analysis of the Phyllostachys edulis forest community after different periods of enclosure.Liso (r): estimated value; Ltrans (r): fitting value; Lbord (r): upper limit of the confidence interval; L pois (r): lower limit of the confidence interval.

Table 1 .
The characteristics of Phyllostachys edulis forests.

Table 2 .
Basic information of the Phyllostachys edulis forest stands.
Note: Tree refers to other trees in the community except Phyllostachys edulis.

Table 3 .
The importance values of mature trees after different periods of natural enclosure.

Table 4 .
The importance values of young trees after different periods of natural enclosure.

Table 5 .
The importance values of shrubs in different years of enclosure.

Table 6 .
The importance values of herbaceous plants after different periods of enclosure.

Table 7 .
Kruskal-Wallis test of mature arbors after different periods of enclosure.

Table 8 .
Diversity indices of mature arbors after different periods of enclosure.
Note: Different lowercase letters indicate significant differences at p < 0.05.

Table 9 .
Kruskal-Wallis test of the young arbors after different periods of enclosure.

Table 10 .
Diversity indices of the young arbors after different periods of enclosure.
Note: Different lowercase letters indicate significant differences at p < 0.05.

Table 11 .
Kruskal-Wallis test of shrubs after different periods of enclosure.

Table 12 .
Diversity index of shrubs after different periods of enclosure.
Note: Different lowercase letters indicate significant differences at p < 0.05.

Table 13 .
Kruskal-Wallis test of the herbaceous plants after different periods of enclosure.

Table 14 .
Diversity index of the herbaceous plants after different periods of enclosure.
Note: Different lowercase letters indicate significant differences at p < 0.05.
Size of the Bamboo Forest Community Structure and Spatial Distribution of Individuals in Each Period of Enclosure A total of 421 individual tree-layer plants with DBH ≥ 1 cm were recorded in the study sample (Figure