Abstract
Forest gaps play an important role during forest succession in temperate forest ecosystems. However, the differences in spatial distribution and replacement patterns of woody plants (trees and shrubs) between primary and secondary forests remain unclear during the gap-filling processes, especially for temperate forests in Northeast China. We recorded 45,619 regenerated trees and shrubs in young gaps (<10 years), old gaps (10~20 years), and closed forest stands (i.e., filled gaps) in the primary broadleaved Korean pine (Pinus koraiensis Sieb. Rt Zucc.) forests vs. secondary forests (degraded from primary forests). The gap-filling processes along horizontal (Cartesian coordinate system) and vertical (lower layer: 0~5 m, medium layer: 5~10 m, and upper layer: >10 m) dimensions were quantified by shade tolerance groups of trees and shrubs. We found that gap age, competition between species, and pre-existing regeneration status resulted in different species replacement patterns within gaps in primary vs. secondary forests. Gap formation in both primary and secondary forests increased species richness, with 33, 38, 39, and 41 in the primary closed stands, primary forest gaps, secondary closed stands, and secondary forest gaps, respectively. However, only 35.9% of species in primary forest gaps and 34.1% in secondary forest gaps successfully reached the upper layer. Based on the importance values (IVs) of tree species across different canopy heights, light-demanding trees in the upper layer of the secondary forests were gradually replaced by intermediate and shade-tolerant trees. In the primary forests, Korean pine exhibited intermittent growth patterns at different canopy heights, while it had continuous regeneration along vertical height gradients in the secondary forests. The differences in Korean pine regeneration between the primary and secondary forests existed before gap formation and continued during the gap-filling processes. The interspecific competition among different tree species gradually decreased with increasing vertical height, and compared to the primary forests, the secondary forests showed an earlier occurrence of competition exclusion within gaps. Our findings revealed the species replacement patterns within gaps and provided a further understanding of the competition dynamics among tree species during the gap-filling processes.
1 Introduction
The escalating issue of forest degradation has emerged as a formidable global challenge confronting humanity (). How to guide degraded forests toward positive succession has become a global concern (). Degraded forests, such as secondary forests, often require forest management to accelerate their restoration to the primary forest with more stability and resilience. Strategies like near-natural forest management, which imitates forest gaps, have been proposed to enhance species diversity and guide positive succession in secondary forests (; ; ). Forest gaps, caused mainly by strong winds, lightning storms, or artificial disturbances, have the potential to alter the above- and below-ground resources in closed forests (; ; ), thereby affecting the structural and compositional progression of the forest (; ).
The intermediate disturbance hypothesis suggested that disturbances, such as forest gaps, promote the coexistence of species with different life histories (; ). The spatial distribution pattern of tree species in canopy gaps represents the competition they are experiencing within the gaps and can also help researchers predict the future development trajectory of canopy gaps (; ; ; ; ). Tree species with different shade tolerance will survive in different positions within the forest canopy (). For example, observed that intermediate shade-tolerant Betula alleghaniensis was mainly located in the south part of the nine-year-old large gaps. Meanwhile, light-demanding P. tabulaeformis preferred regenerating in the northeastern part of the gap created seven years ago (). Understanding the spatial distribution patterns of different species within forest gaps can provide important clues regarding the future development trends of the forest and serve as a scientific basis for ecological conservation and management (; ).
Moreover, the resource competition theory suggests that the horizontally and vertically spatial distribution patterns of species are also determined by resource competition (; ). The types of competition experienced by individual trees vary at different vertical heights within canopy gaps (). In the lower layers of canopy gaps, shrubs often compete for growing space and negatively impact the growth of tree seedlings (). In the middle and upper layers of canopy gaps, the competition between tree species gradually decreases, and the spatial distribution pattern of trees is controlled by resource heterogeneity at the landscape scale (; ). Only a few tree species can successfully grow from lower to middle and upper layers (). The average tree species regeneration from a horizontal perspective may lead to a limited estimation of the progress of tree species succession in forest gaps (). The vertical arrangement of temperate forests holds significant importance once the trees have been established (). According to , the vertical structure plays a crucial role in the regeneration of populations and the dynamics of communities. It also provides valuable insights into population dynamics and competition in different height layers. Therefore, it is necessary to have a clear understanding of the horizontal and vertical regeneration processes of different shade-tolerant tree species within forest gaps (; ).
Secondary forests formed due to severe disturbances in primary forests account for over 70% of the forests in Northeast China (). Unlike primary forests, secondary forests have the characteristics of destruction or loss of basic structure and inherent functions (), decreased biodiversity (), weakened stability and resilience (), and low system productivity (). For example, observed that the Shannon diversity of the primary broadleaved Korean pine forest was higher than that in the secondary forest. Previous studies indicated that it would take over 200 years for secondary forests to restore functions similar to primary forests under natural development (). The current issue for forest managers is how to guide the positive succession of secondary forests to restore the forest structure to resemble that of the primary forests. In recent decades, forest managers have developed various forest management programs, especially simulating natural gap disturbances, to accelerate the restoration of degraded forests and enhance ecosystem services (). Consequently, understanding the process and mechanisms of gap regeneration is important to achieve the management goals of promoting natural forest succession (; ).
This study aims to compare the horizontal and vertical regeneration patterns of tree species with different shade tolerances (i.e., light-demanding, intermediate, and shade-tolerant) in the primary and secondary forests during the gap-filling processes. We used point pattern analysis methods to explore the spatial distribution patterns of trees and shrubs at different vertical layers and the underlying implications for competition. We investigated all trees and shrubs within 12 natural gaps in the primary and secondary forests. The specific questions were: (1) How does the regeneration of tree species change along vertical dimensions during gap-filling in primary vs. secondary forests? (2) How do tree species of different shade tolerance fill the gaps regarding distribution patterns at different vertical layers? (3) How does the interspecific competition change at different vertical layers in primary vs. secondary forest gaps? We hypothesized that: (1) As the gap closure process progresses, the IV of trees in both forest types will increase, but species-specific differences exist. (2) Light-demanding species dominate in the early stages of gap formation, while intermediate and shade-tolerant species gradually replace part of the light-demanding species. (3) Compared to primary forests, competition for trees in secondary forests might be less intense.
2 Materials and methods
2.1 Study area
This study was carried out within the Lushuihe Forest Bureau, situated in the northwestern region of the Changbai Mountains (127°29’-128°02’ E, 42°20’-42°40’ N). The study area is characterized by mixed broadleaved Korean pine forests, and it experiences a typical temperate continental climate. The region has an average annual temperature of 2.8°C and an annual precipitation ranging from 800 to 1040 mm. The elevation in this region ranges from 600 to 800 meters above mean sea level, and the soil is typical mountainous dark brown forest soil (). Before experiencing significant disturbances, the forest type in the Lushuihe Forest Bureau was the primary broadleaved Korean pine forest. The dominant species were Korean pine and broadleaved trees such as Tilia amurensis, Acer mono, Fraxinus mandshurica, and Juglans mandshurica (). The secondary forest in our study area is derived from natural regeneration after clear-cutting and is dominated by Betula platyphylla. In 1998, the Chinese government initiated implementing Natural Forest Conservation Projects, and the Lushuihe Forestry Bureau began prohibiting timber logging ().
2.2 Gap selection
We set up six treatments to compare the gap regeneration in different development stages, including primary young gap (PYG), primary old gap (POG), primary closed stand (PCK), secondary young gap (SYG), secondary old gap (SOG) and secondary closed stand (SCK). Each treatment was repeated three times.
We adopted the belt transect method in the field investigation in 2019. A total of 198 and 30 forest gaps were pre-selected in the primary broadleaved Korean pine forest and in the secondary birch forest, respectively. The purpose of our pre-selection was to identify ideal forest gaps (and CK plots). For all pre-selected gaps, we measured gap size, recorded gap shape, and estimated gap age using tree ring analysis, and then we decided which gap could be used. The minimum distance between any two sample plots was more than 100 m, and there were no other gaps between any two selected plots. Then, we only selected forest gaps that meet the following conditions for field surveys: (1) The diameter-to-height ratio of forest gaps >0.73 (medium gap and large gap); (2) No obvious signs of human or animal disturbances; (3) The tree core shows obvious growth and release, which can accurately determine the age of forest gaps (). Consequently, six medium and large gaps were selected from the primary broadleaved Korean pine and secondary forests. Three gaps in the primary forest were 440.7 ± 153.7 m2 on average with a gap age of 6.3 ± 1.7 years as the PYG, and three gaps were 612.7 ± 115.2 m2 on average with a gap age of 15 ± 1.6 years old as the POG. Three gaps in the secondary forest were 395.6 ± 80.7 m2 on average with a gap age of 6 ± 1.4 years old as the SYG, and three gaps were 417.6 ± 78.6 m2 on average with a gap age of 14 ± 2.2 years old as the SOG. We selected three control plots (30 m × 30 m) in PCK and SCK as the plot filled and closed through forest gaps.
2.3 Gap regeneration surveys
We recorded a total of 45,196 trees and shrubs with height >0.1 m in the plots, documenting all their species, coordinates, heights, diameters at breast height (DBH) or root collar diameter (RCD), and other relevant indicators. We divided them into three height groups, including the lower layer (height: 0~5 m), medium layer (height: 5~10 m), and upper layer (height: >10 m) within the PYG, POG, PCK, SYG, SOG, and SCK plots. We established a coordinate system with the center of the forest gap as the origin (0, 0). The north-south direction of the forest gap is the Y-axis, and the east-west direction of the forest gap is the X-axis. To map the regeneration locations more accurately, we divided every plot into 100 subplots (3 m × 3 m) to avoid missing individuals during the investigation. An ultrasound-based positioning instrument (Fieldscout TDR 350, Haglöf PosTex Långsele, Sweden) was used to acquire the precise coordinates of each individual, which could improve the accuracy of spatial point pattern analysis results.
2.4 Statistical analysis
2.4.1 Regeneration within gaps
To assess the transition patterns of tree species in different stages of gap closure between primary and secondary forests, we calculated the importance values of different tree species in PYG, POG, PCK, SYG, SOG, and SCK. In order to explore the regeneration dynamics of different tree species at varying vertical heights during different stages of gap closure in primary and secondary forests, we computed the importance values of different tree species in lower layer (height: 0~5 m), medium layer (height: 5~10 m) and upper layer (height: >10 m) within the PYG, POG, PCK, SYG, SOG, and SCK plots. The importance value was calculated with the following equations:
where RD is the relative dominance, RF is the relative frequency, RC is the relative coverage, ni is the number of individuals of the ith species, ai is the basal area belonging to the ith species, fi is the number of quadrats the ith species appeared, and S is the total number of species.
To assess the regeneration process in gaps of different vertical layers, we calculated the species richness in the lower layer (height: 0~5 m), medium layer (height: 5~10 m), and upper layer (height: >10 m) of PYG, POG, PCK, SYG, SOG, and SCK plots. The species richness was calculated as the number of species in each plot.
To assess the regeneration process in gaps of different ages, we calculated the Shannon diversity index of PYG, POG, PCK, SYG, SOG, and SCK plots. The Shannon diversity index was calculated with the following equations:
where Pi is the proportion of the entire species community, and S is the total number of species in plots.
To assess the regeneration growth in gaps of different ages, we analyzed the height differences of light-demanding, intermediate, and shade-tolerant species in PYG, POG, PCK, SYG, SOG, and SCK plots. One-way ANOVA was used to test the effects of gap age on the heights of woody species. Tukey’s post hoc tests were used to test the differences between different gap ages using SPSS (SPSS software, 22nd edition, Chicago, USA).
2.4.2 Univariate and bivariate spatial patterns
Point pattern analysis is widely applied in describing gap regeneration patterns (; ). The univariate spatial patterns of trees at the lower layer (height: 0~5 m), medium layer (height: 5~10 m), and upper layer (height: >10 m) in PYG, POG, PCK, SYG, SOG, and SCK plots were analyzed using g(r). Point pattern analysis was conducted using the Programita software ().
Many previous studies used K(r) to calculate the spatial pattern, but K(r) has a nature of cumulation (). The pair correlation function g(r) is widely used because of its non-cumulative character. The pair correlation function employs annuli as distance categories rather than circles. g(r) is related to the K(r) () as follows:
where s is the area of the forest gaps, n is the number of points, m and n are all ordered point pairs, I is an indicator (1 if dmn ≤ 1), dmn is the distance between two points, and e is an edge correction. Points above the upper g(r) envelope represent clustered, and points below the g(r) envelope represent regularity. The g(r) envelopes were derived from 199 Monte Carlo simulations of CSR ().
The bivariate relationship between trees and shrubs in PYG, POG, PCK, SYG, SOG, and SCK plots was analyzed by g12(r):
If the observed g12(r) was located above or below the simulation envelope at a specific distance r, indicating the attraction or repulsed pattern of shrubs and trees, respectively. If the observed g12(r) was located in the simulation envelope, indicating the independent relation of shrubs and trees. In the context of point pattern analysis, the significance level of 199 Monte-Carlo test was 0.02.
2.4.3 Kernal density estimation for density
Kernel density estimation (KDE) was used to analyze the density distribution of light-demanding, intermediate, and shade-tolerant trees in lower, medium, and upper layers of PYG, POG, PCK, SYG, SOG, and SCK ().
The KDE analysis, introduced by , is a nonparametric statistical method. The KDE analysis utilizes the quadratic kernel function to represent a smooth surface for every data point. The KDE analysis enables the computation of the magnitude per unit area, as described by . The kernel density estimator is defined as follows:
where S is a non-negative function with an integral of 1, with an average value of zero, p is the number of the points of trees and shrubs (x), and r is the search radius (). The research radius r is 1 m in our research. We converted the initial coordinates of the trees and shrubs into a new coordinate system, with the origin coordinates set at the center of each gap. These new coordinates system effectively represented all 45,196 trees and shrubs in the plots, with the forest gap center points as a common anchor point. Trees and shrubs were first grouped by layers (lower layer, medium, and upper layer) and shade-tolerance (light-demanding, intermediate, and shade-tolerant) and were aligned based on the coordinates of the center of each gap. Then, we ran kernel density estimates for the 12 gaps and 6 closed forest stands. ArcMap 10.1 (ESRI, 2012) was used to calculate the kernel density estimates for PYG, POG, PCK, SYG, SOG, and SCK.
3 Results
3.1 General regeneration patterns
Generally, the number of species in the PCK, the primary forest gap (PYG and POG), SCK, and the secondary forest gap (SYG and SOG) were 33, 39, 38, and 41, respectively (Table S2). The Shannon diversity index of PCK, PYG, POG, SCK, SYG, and SOG were 2.50, 2.41, 2.75, 2.54, 2.08 and 2.67, respectively. However, the IV changes in tree species in the primary and secondary forests were quite different (Table S1). According to the data in Figure 1, the IV of Korean pine in the primary forest gap was far lower than that in the secondary forest gap, which was 1.4% in PYG, 2.3% in POG, 6.9% in SYG and 10.4% in SOG, respectively. Acer spp. maintained high IV in the gap plots and control plots of the primary forest or secondary forest (Figure 1).
Figure 1
3.2 Regeneration pattern in different layers
In both primary and secondary forest gaps, there was a trend of decreasing species richness with increasing vertical height (Figure 2). However, the IV of some species in primary and secondary forests showed a distinct pattern across lower, medium, and upper layers. As was apparent from Figure 2, Korean pine had the highest importance value in the medium layer of secondary forest gaps but barely had any regeneration in primary forest gaps. Acer spp. had a high proportion in three layers of primary and secondary forests. Surprisingly, the importance value of B. platyphylla in the upper layer of the SCK was 50%, while its importance value in the lower and medium layers of the SCK tended to be zero. Similar results were also found in the primary forest.
Figure 2
From the pattern analysis results in Figures 3, 4, it can be seen that trees tended to randomly distribute as the height layer increased. Trees in the upper layer were spatially random across all scales (100%) examined in PYG, POG, PCK, SYG, SOG, and SCK (Figures 3, 4). Similarly, trees in the lower layer were clustered at all distances (100%) examined in PYG, POG, PCK, SYG, and SOG, and 50% scale in SCK. However, the trees in the medium layer clustered in PYG, POG, and PCK at a small scale (Figure 3), contrasting with a random distribution at 100% distances in the medium layers of SOG and SCK (Figure 4). The density distribution of light-demanding and intermediate trees was denser than shade-tolerant trees (Figure 5). As the process of gap closure advanced, there was a gradual decline in the highest density of light-demanding trees in the lower layer appealing, with a discernible pattern of decreasing from 8-11 stems m-2 in the PYG to 3-4 stems m-2 in the POG, and ultimately reaching 1-2 stems m-2 in the PCK (Figure 5). However, the density of shade-tolerant trees in the lower layer of secondary forests remains relatively stable (Figure 5). In both primary and secondary forest gaps, there was a similar pattern of the density of intermediate trees increasing with the age of the gaps in the lower layer regeneration (Figure 5). The height of upper layer trees increased gradually with the gap closure of primary and secondary forests (Figures 6, 7). However, the average height of light-demanding, intermediate, and shade-tolerant trees remained stable in the medium layer of primary forest and secondary forest (Figures 6, 7). Competition between trees and shrubs in the lower layer showed a small-scale positive correlation and a large-scale negative correlation in PYG, SYG, and PCK (Figure 8). The mutual exclusion between trees and shrubs in primary forests intensified with forest succession. Compared with secondary forests, the competition relationship between trees and shrubs was more obvious in primary forests (Figure 8).
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
4 Discussion
4.1 General regeneration patterns
We found that forest gaps increased species richness in both primary and secondary forests. Gap formation increases light penetration, creating opportunities for plants to thrive and for new seedlings to establish (). This led to a greater variety of habitats and resources, ultimately supporting a more comprehensive range of species in the forest ecosystem (; ). In accordance with the findings of , our study also observed a similar pattern of Shannon diversity index in primary and secondary forests, namely: PYG<PCK<POG and SYG<SCK<SOG. This pattern can be explained by the early stages of forest gap formation, during which dominant species had more regeneration compared to other species, resulting in a lower Shannon diversity index in young gaps. As the forest gaps aged, the diversity of tree species increased, leading to a rise in the Shannon diversity index (). Our results showed that the variation in importance values of the same tree species differs between plots in the primary and secondary forests. The importance value of Korean pine was much lower in the primary forest gaps than in the secondary forest gaps. To investigate the reasons, we inferred the advanced regeneration of Korean pine before gap formation based on the calculated age of the Korean pine (). The advantage of Korean pine regeneration in the secondary forest was already present before gap formation (Figure S1). After gap formation, the regenerating Korean pine maintained this advantage (). Previous results showed a positive correlation between the growth of Korean pine seedlings and canopy transparency (). The light intensity in the secondary forest gap was higher than in the primary forest gap, which led to the good regeneration of Korean pine in the secondary forest gap. claimed that habitat filtering and intraspecific density constraints simultaneously affect seedling survival. The other reason for the poor regeneration of Korean pine in the primary forest might be the high proportion of Korean pine within the primary broadleaved Korean pine forest, resulting in the conspecific negative density dependence phenomenon (; ).
Conversely, the importance value of Acer spp. was notably high in the forest gaps and closed stands of the primary forest, and a similar trend was observed in the secondary forest. (). Some studies reported that the regeneration of Acer spp. was not sensitive to canopy disturbances (), which was consistent with our results. Although Acer spp. were not sensitive to the formation of canopy gaps in both primary and secondary forests, they had high IVs in all plots. This suggested that Acer spp. might become the main tree species in the mixed broadleaved Korean pine forest.
4.2 Regeneration pattern in different layers
The proportion of importance values of Korean pine at different height layers in the primary forest plots differed from those at different height layers in the secondary forest plots. It was consistent that the Korean pine’s medium layers () were barely absent in the primary forest gap and closed stand (). In our research, the importance value of Korean pine was relatively high in the medium layer of secondary forests but relatively low in the lower and upper layers. This meant Korean pine continuously filled gaps from the lower layer to the medium layer and then to the upper layer in secondary forests. The species richness of shrubs and trees in the understory of secondary forests was high. As the fill-in process of the medium and large forest gaps, only several tree species could successfully ascend from the lower layer to the medium layer and then extend from the medium layer to the upper layer, ultimately completing the filling process of the canopy. Our results showed that the IV of shade-tolerant trees in the upper layers of the SCK (1.7%) was far lower than that of shade-tolerant trees in the upper layers of the PCK (14.8%). In the lower layer of the primary forest gap, there was a decline in the density of light-demanding trees as the age of the forest gap increased. Conversely, there was an increase in the density of intermediate and shade-tolerant trees. Similar to our research findings, in canopy gaps within temperate hardwood forests in Connecticut, USA, there was also an observed phenomenon where light-demanding species initially dominated but were eventually replaced by shade-tolerant species (). This phenomenon was likely caused by the suitable light conditions after gap formation, promoting the growth of light-demanding species, which subsequently made the light conditions more favorable for the growth of shade-tolerant tree species (). However, there were also studies suggesting that the differences in growth rates between tree species with varying light requirements were insufficient in the short term to result in the complete replacement of one species by another (; ). Therefore, based on the current trend of species importance value changes, we predicted that the proportion of shade-tolerant trees in the upper layer would increase with forest succession.
The formation of gaps had expanded the scale of tree aggregation, with the aggregation scale of medium layers ranging from 0-1m in PCK,0-4 m in PYG, and 0-6 m in POG. Given these variations in scale, it is necessary to focus on the competitive individuals of the dominant tree species within these ranges and promptly implement forest management practices (). suggested that an aggregated distribution pattern implied that individual trees were still engaged in intense competition, while a random distribution pattern was the consequence of community-level competition. Our research findings suggested that there was still competition among trees with different scales in the medium layer of primary forest gaps. The intensity of this competition was greater within old forest gaps compared to young forest gaps and control plots. In the medium layer of older gaps and closed stands of secondary forests, trees had completed their competition and exhibited a random distribution pattern. Some studies suggested that the spatial pattern of trees was related to the gap size. For example, indicated that as the gap size increased (from 49.3 m2 to 185.8 m2), the aggregation of Chinese pine in forest gaps increased. Other researchers argued that the aggregated distribution was due to heterogeneity within the study area or the dispersal patterns of seeds (; ). Our findings indicated that the gap formation in both primary and secondary forests did not affect the clustered distribution in the lower layer and random distribution in the upper layer. The likely reason for the difference in distribution patterns in the medium layer of primary and secondary forest gaps was that the growing conditions in secondary forests were more suitable for tree growth, while the resources within primary forest gaps were limited.
The differences in height and density distribution across vertical layers also indicated differences in the primary and secondary forest regeneration processes (). The average height of shade-tolerant trees in the lower layer of the POG was significantly higher than that in the PYG and PCK. With the gap-filling progress, the density distribution of light-demanding trees gradually decreased in the lower layer. Similarly, the average height of intermediate trees in the lower layer of SOG was significantly higher than that in the SYG and SCK. However, the density distribution of intermediate trees in the lower layer remained unchanged. These trends indicated that light-demanding tree species were more likely to access the medium layer at an earlier stage within the lower layer of primary forest gaps. Theoretically, differences in tree height growth contributed to niche differentiation, ultimately allowing for the coexistence of multiple species (). observed a significant decrease in the density of light-demanding species in the lower layer, suggesting a lack of new members to complement the population. Consequently, it can be inferred that light-demanding species may only have a single opportunity for success, but intermediate species maintain continuous renewal in forest gaps.
4.3 Competition between trees and shrubs
Our findings indicated that the correlation between shrubs and trees varied with scale, exhibiting a trend of positive correlation at small scales in PYG, POG, PCK, SYG, SOG, and SCK and negative correlation at large scales in PYG and PCK. Previous research showed that the correlation pattern at small scales was attributed to biological characteristics (). Our research indicated that the formation of forest gaps did not change the positive or independent relationship between shrubs and trees on a small scale (0-4 m). Some studies also suggested that the presence of shrubs promoted the regeneration of trees (; ). Shrubs could relieve non-biological stress by reducing solar radiation, wind speed, and soil temperature, thereby reducing water loss caused by plant transpiration and soil moisture loss caused by evaporation (). Shrubs could also inhibit tree growth by blocking moss’s direct inhibitory effect on them (). In addition, shrubs could indirectly promote the growth of tree seedlings by protecting them from herbivorous animals (). However, the negative correlation between trees and shrubs at a large scale in the primary forest plots was opposite to the positive association at a large scale in the secondary forest plots. This could be related to the differences in environmental conditions within the primary and secondary forests. The correlation effect at large scales was due to environmental factors (). suggested that the high amount of shrub regeneration directly offsets the quality of the regeneration environment provided by forest gaps. Removing shrubs from primary and secondary forests could promote seedling growth of shade-tolerant trees (). Once a seedling’s height growth surpasses the shrub layer, it may have a higher chance of successfully entering the canopy layer, as there are fewer potential competitors in the forest gap. Understory shrubs reduced seedling recruitment by shading out light (). Therefore, in both primary and secondary forest plots, shrubs protected tree regeneration at small scales. However, at large scales, shrubs competed with trees for growth resources in newly formed forest gaps.
5 Conclusions
In this study, we used point pattern analysis to compare the gap-filling processes of tree species in primary and secondary forest gaps across horizontal and vertical dimensions. We found that the formation of canopy gaps in both primary forests and secondary forests indeed increased tree species richness. However, the main tree species in the primary broad-leaved Korean pine forest, such as Quercus mongolica, F. mandshurica, and T. amurensis, have not regenerated well in the lower layer. Q. mongolica, F. mandshurica, and T. amurensis could be planted in the forest management process. At the same time, to end the aggregation of the trees in the middle and lower layers of the forest gaps in advance, non-main tree species could be removed to promote the growth of the main trees. The regeneration advantage of Korean pine in secondary forest plots existed before and after the formation of forest gaps. In secondary forests, Korean pine successfully grew from the lower to the upper layers, whereas it failed to do so in primary forests. The competition differed between tree species in canopy gaps of primary and secondary forests manifested in the medium layer. In secondary forest gaps, shade-tolerant and intermediate species replaced the regeneration advantage of light-demanding trees. In primary forest gaps, the regeneration advantage of intermediate trees were consistently maintained. Our findings further explained the spatial pattern of trees in the gap-filling phase.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
DW: Data curation, Formal Analysis, Investigation, Writing - original draft, Writing - review & editing. DL: Data curation, Formal Analysis, Methodology, Writing - original draft, Writing - review & editing. JZ: Conceptualization, Funding acquisition, Methodology, Project administration, Writing - review & editing. XG: Data curation, Investigation, Writing - review & editing. JZ: Data curation, Investigation, Writing - review & editing. LL: Formal analysis, Writing - review & editing. XW: Writing - review & editing. HL: Data curation, Investigation, Writing - review & editing. GZ: Methodology, Writing - review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was financially supported by the National Natural Science Foundation of China (31830016 and 32001301), the Youth Innovation Promotion Association of CAS (2023205), and the CAS Foundation for Special Research Assistant (2019).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1305535/full#supplementary-material
References
1
Abd LatifZ.BlackburnG. A. (2010). The effects of gap size on some microclimate variables during late summer and autumn in a temperate broadleaved deciduous forest. Int. J. Biometeorol.54, 119–129. doi: 10.1007/s00484-009-0260-1
2
AgyemanV. K.SwaineM. D.ThompsonJ. (1999). Responses of tropical forest tree seedlings to irradiance and the derivation of a light response index. J. Ecol.87, 815–827. doi: 10.1046/j.1365-2745.1999.00400.x
3
BeckageB., ClarkClintonB. D.HainesB. L. (2000). A long-term study of tree seedling recruitment in southern appalachian forests: the effects of canopy gaps and shrub understories. Can. J. For. Res.30, 1617–1631. doi: 10.1139/x00-075
4
BowmanA. W. (1985). A comparative study of some kernel-based nonparametric density estimators. J. Stat. Comput. Simul.21, 313–327. doi: 10.1080/00949658508810822
5
DaiL.LiS.ZhouW.QiL.ZhouL.WeiY. (2018). Opportunities and challenges for the protection and ecological functions promotion of natural forests in China. For. Ecol. Manage.410, 187–192. doi: 10.1016/j.foreco.2017.09.044
6
DekkerM.Sass-KlaassenU.den OudenJ.GoedhartP. W. (2009). The effect of canopy position on growth and mortality in mixed sapling communities during self-thinning. Eur. J. For. Res.128, 455–466. doi: 10.1007/s10342-009-0293-3
7
FanX. H.HeH. J.WangJ.XuC. Y.GadowK. V. (2017). Tree species richness in northeastern China: geographical variation and minimum areas. Community Ecol.18, 244–252. doi: 10.1556/168.2017.18.3.3
8
FanX. H.XuC. Y. (2019). Formation mechanism of taxonomic structures for forest communities in different successional stages in Changbai Mountains of northeastern China. J. Beijing For. Univ41, 24–32. doi: 10.13332/j.1000–1522.20180301
9
FineganB. (1996). Pattern and process in neotropical secondary rain forests: the first 100 years of succession. Trends Ecol. Evolution.11, 119–124. doi: 10.1016/0169-5347(96)81090-1
10
GeX.ZhuJ.LuD.ZhuC.GaoP.YangX. (2022). Age calibration and growth stage division for naturally regenerated seedlings and saplings of Pinus koraiensis Sieb. et Zucc. Chin. J. Eco.41, 417–426. doi: 10.13292/j.1000-4890.202202.016
11
FAO. (2020). Global Forest Resources Assessment 2020: report. Rome: FAO. doi: 10.4060/ca9825en
12
GrubbP. J. (1977). The maintenance of species-richness in plant communities: the importance of the regeneration niche. Biol. Rev.52, 107–145. doi: 10.1111/j.1469-185X.1977.tb01347.x
13
GuariguataM. R.OstertagR. (2001). Neotropical secondary forest succession: changes in structural and functional characteristics. For. Ecol. Manage.148, 185–206. doi: 10.1016/S0378-1127(00)00535-1
14
GuldinJ. M.LorimerC. G. (1985). Crown differentiation in even-aged northern hardwood forests of the Great Lakes region, U.S.A. For. Ecol. Manage.10, 65–86. doi: 10.1016/0378-1127(85)90014-3
15
HaleS. E. (2003). The effect of thinning intensity on the below-canopy light environment in a Sitka spruce plantation. For. Ecol. Manage.179, 341–349. doi: 10.1016/S0378-1127(02)00540-6
16
HaoZ.ZhangJ.SongB.YeJ.LiB. (2007). Vertical structure and spatial associations of dominant tree species in an old-growth temperate forest. For. Ecol. Manag.252, 1–11. doi: 10.1016
17
HolmgrenM.LinC.-Y.MurilloJ. E.NieuwenhuisA.PenninkhofJ.SandersN.et al. (2015). Positive shrub-tree interactions facilitate woody encroachment in boreal peatlands. J. Ecol.103, 58–66. doi: 10.1111/1365-2745.12331
18
HuaF.BruijnzeelL. A.MeliP.MartinP. A.ZhangJ.NakagawaS.et al. (2022). The biodiversity and ecosystem service contributions and trade-offs of forest restoration approaches. Science.376, 839–844. doi: 10.1126/science.abl4649
19
HubbellS. P. (2001). The Unified Neutral Theory of Biodiversity and Biogeography (Princeton: Princeton University Press).
20
HubbellS. P.FosterR. B.O’BrienS. T.HarmsK. E.ConditR.WechslerB.et al. (1999). Light-gap disturbances, recruitment limitation, and tree diversity in a neotropical forest. Science.283, 554–557. doi: 10.1126/science.283.5401.554
21
JacquemynH.ButayeJ.HermyM. (2003). Influence of environmental and spatial variables on regional distribution of forest plant species in a fragmented and changing landscape. Ecography.26, 768–776. doi: 10.1111/j.0906-7590.2003.03620.x
22
KernC. C.BurtonJ. I.RaymondP.D’AmatoA. W.KeetonW. S.RoyoA. A.et al. (2017). Challenges facing gap-based silviculture and possible solutions for mesic northern forests in North America. Forestry.90, 4–17. doi: 10.1093/forestry/cpw024
23
KernC. C.ReichP. B.MontgomeryR. A.StrongT. F. (2012). Do deer and shrubs override canopy gap size effects on growth and survival of yellow birch, northern red oak, eastern white pine, and eastern hemlock seedlings? For. Ecol. Manage.267, 134–143. doi: 10.1016/j.foreco.2011.12.002
24
LaMannaJ. A.ManganS. A.MyersJ. A. (2021). Conspecific negative density dependence and why its study should not be abandoned. Ecosphere12(1):e03322. doi: 10.1002/ecs2.3322
25
LangfordM.UnwinD. J. (1994). Generating and mapping population density surfaces within a geographical information system. Cartogr. J.31, 21–26. doi: 10.1179/000870494787073718
26
Lara-RomeroC.de la CruzM.Escribano-ÁvilaG.García-FernándezA.IriondoJ. M. (2016). What causes conspecific plant aggregation? Disentangling the role of dispersal, habitat heterogeneity and plant-plant interactions. Oikos.125, 1304–1313. doi: 10.1111/oik.03099
27
LarocqueG. R.LuckaiN.AdhikaryS. N.GrootA.BellF. W.SharmaM. (2013). Competition theory - science and application in mixed forest stands: review of experimental and modelling methods and suggestions for future research. Environ. Rev.21, 71–84. doi: 10.1139/er-2012-0033
28
LiJ.WangY. (1986). Waves features of population changes of pinus koraiensis in natural forest. J. Ecol.5, 1–5. doi: 10.13292/j.1000-4890.1986.0078
29
LimpensJ.Van EgmondE.LiB.HolmgrenM. (2014). Do plant traits explain tree seedling survival in bogs? Funct. Ecol.28, 283–290. doi: 10.1111/1365-2435.12148
30
LiptzinD.AshtonP. M. S. (1999). Early-successional dynamics of single-aged mixed hardwood stands in a southern New England forest, USA. For. Ecol. Manage.116, 141–150. doi: 10.1016/S0378-1127(98)00448-4
31
LiuJ.HongW.XuZ.YuL.LiJ. (2003). Study on dynamic pattern of species diversity in gaps of Castanopsis Kawakamii Forest. Scientia. Silvae. Sinicae.39, 159–164.
32
LiuC.LiZ. (1987). Habitat under Birch-Aspen Forest and Natural Recruitment of Korean pine. J. Northeast. For. Univ.15, 22–28. doi: 10.13759/j.cnki.dlxb.1987.04.004
33
LortieC. J.CushmanJ. H. (2007). Effects of a directional abiotic gradient on plant community dynamics and invasion in a coastal dune system. J. Ecol.95, 468–481. doi: 10.1111/j.1365-2745.2007.01231.x
34
LuD.WangG. G.YuL.ZhangT.ZhuJ. (2018). Seedling survival within forest gaps: the effects of gap size, within-gap position and forest type on species of contrasting shade-tolerance in Northeast China. For. Int. J. For. Res.91, 470–479. doi: 10.1093/forestry/cpy007
35
LuD.ZhuJ.ZhangG.SunX.SunY.HuL.et al. (2023). Disentangling regeneration by vertical stratification: A 17-year gap-filling process in a temperate secondary forest. For. Ecol. Manage.539, 120994. doi: 10.1016/j.foreco.2023.120994
36
MartinM.WoodburyD.GlogowerY.DuguidM.FreyB.AshtonM. (2021). Within-gap position shapes fifty years of forest dynamics in a temperate hardwood forest in Connecticut, USA. For. Ecol. Manage.494, 119311. doi: 10.1016/j.foreco.2021.119311
37
McClureJ. W.LeeT. D. (1993). Small-scale disturbance in a northern hardwoods forest: effects on tree species abundance and distribution. Can. J. For. Res.23, 1347–1360. doi: 10.1139/x93-172
38
NakashizukaT.MatsumotoY. (2002). Diversity and interaction in a temperate forest community, ecological studies (Tokyo: Springer Japan). doi: 10.1007/978-4-431-67879-3
39
OmelkoA.UkhvatkinaO.ZhmerenetskyA.SibirinaL.PetrenkoT.BobrovskyM. (2018). From young to adult trees: How spatial patterns of plants with different life strategies change during age development in an old-growth Korean pine-broadleaved forest. For. Ecol. Manage.411, 46–66. doi: 10.1016/j.foreco.2018.01.023
40
OwenS. M.SiegC. H.Sánchez MeadorA. J.FuléP. Z.IniguezJ. M.BaggettL. S.et al. (2017). Spatial patterns of ponderosa pine regeneration in high-severity burn patches. For. Ecol. Manage.405, 134–149. doi: 10.1016/j.foreco.2017.09.005
41
PoznanovicS. K.PoznanovicA. J.WebsterC. R.BumpJ. K. (2014). Spatial patterning of underrepresented tree species in canopy gaps 9 years after group selection cutting. For. Ecol. Manage.331, 111. doi: 10.1016/j.foreco.2014.06.029
42
RaymondP.MunsonA. D.RuelJ. C.CoatesK. D. (2006). Spatial patterns of soil microclimate, light, regeneration, and growth within silvicultural gaps of mixed tolerant hardwood white pine stands. Can. J. For. Res.36, 639–651. doi: 10.1139/x05-269
43
RipleyB. D. (1977). Modelling spatial patterns. J. R. Stat. Soc Ser. B Methodol.39, 172–192. doi: 10.1111/j.2517-6161.1977.tb01615.x
44
RozasV. (2003). Tree age estimates in Fagus sylvatica and Quercus robur: testing previous and improved methods. Plant Ecol.167, 193–212. doi: 10.1023/A:1023969822044
45
SeidlR.ThomD.KautzM.Martin-BenitoD.PeltoniemiM.VacchianoG.et al. (2017). Forest disturbances under climate change. Nat. Clim. Change.7, 395–402. doi: 10.1038/nclimate3303
46
SilvermanB. W. (1986). Density Estimation for Statistics and Data Analysis (London: Chapman & Hall). doi: 10.1007/978-1-4899-3324-9
47
TesfayeG.TeketayD.FeteneM.BeckE. (2010). Regeneration of seven indigenous tree species in a dry Afromontane forest, southern Ethiopia. Flora.205, 135–143. doi: 10.1016/j.flora.2008.12.006
48
ThomD.SeidlR. (2016). Natural disturbance impacts on ecosystem services and biodiversity in temperate and boreal forests. Biol. Rev.91, 760–781. doi: 10.1111/brv.12193
49
TurnerI. M.WongY. K.ChewP. T.IbrahimA. B. (1997). Tree species richness in primary and old secondary tropical forest in singapore. Biodivers. Conserv.6, 537–543. doi: 10.1023/A:1018381111842
50
UrzaA. K.WeisbergP. J.ChambersJ. C.SullivanB. W. (2019). Shrub facilitation of tree establishment varies with ontogenetic stage across environmental gradients. New Phytol.223, 1795–1808. doi: 10.1111/nph.15957
51
WangS.XuJ.LiX.ZhaoS.ChenG.LiG. (1993). Analysis of life table method on the potential restoration of broad-leaved korean pine forests in the xiaoxing'an mountains. for. Sci. Technol.18, 1–4.
52
WangJ.TengD.HeX.QinL.YangX.LvG. (2021). Spatial non-stationarity effects of driving factors on soil respiration in an arid desert region. Catena.207, 105617. doi: 10.1016/j.catena.2021.105617
53
WangZ.YangH.DongB.ZhouM.MaL.JiaZ.et al. (2017). Effects of canopy gap size on growth and spatial patterns of Chinese pine (Pinus tabulaeformis) regeneration. For. Ecol. Manage.385, 46–56. doi: 10.1016/j.foreco.2016.11.022
54
WiegandT.GrabarnikP.StoyanD. (2016). Envelope tests for spatial point patterns with and without simulation. Ecosphere.7, e01365. doi: 10.1002/ecs2.1365
55
YeJ.HaoZ.WangX.BaiX.XingD.YuanZ. (2014). Local-scale drivers of multi-stemmed tree formation in Acer, in a temperate forest of Northeast China. Chin. Sci. Bull.59, 320–325. doi: 10.1007/s11434-013-0013-8
56
YinZ.ZhangC.FanX.ZhangN.ZhaoX.Von GadowK. (2023). Effects of tree density and herbaceous plants on tree seedling survival across the growing and non-growing season in a temperate forest. For. Ecol. Manage.545, 121234. doi: 10.1016/j.foreco.2023.121234
57
YuZ.HaoZ. (1998). Canopy gap characteristics and its influence on the regeneration of broad-leaved Korean pine forests in Changbai Mountain. J. For. Res.9, 160–165. doi: 10.1007/BF02910062
58
YuD.ZhouW.ZhouL.DaiL. (2019). Exploring the history of the management theory and technology of broad-leaved Korean pine (Pinus Koraiensis Sieb. et Zucc.) forest in Changbai Mountain Region, Northeast China. Chin. J. Appl. Ecol.30, 1426–1434. doi: 10.13287/j.1001-9332.201905.004
59
ZhangC.ZhuJ.YanQ.YuL.SunY.XiX. (2009). Comparison of plant species diversity and factors influencing species diversity of primeval and secondary forests on the northern slope of chanbai mountains, china. J. Northeast For. Univ.37, 52–55. doi: 10.13759/j.cnki.dlxb.2009.11.024
60
ZhangJ.GeS.LiangJ.LiJ. (2022). Population age structure and dynamics of Pinus koraiensis in a broadleaved Korean pine forest in Changbai Moutain, China. Chin. J. Plant Ecol.46, 667–677. doi: 10.17521/cjpe.2021.0498
61
ZhangM.ZhuJ.LiM.ZhangG.YanQ. (2013). Different light acclimation strategies of two coexisting tree species seedlings in a temperate secondary forest along five natural light levels. For. Ecol. Manage.306, 234–242. doi: 10.1016/j.foreco.2013.06.031
62
ZhaoF.YangJ.HeH. S.DaiL. (2013). Effects of natural and human-assisted regeneration on landscape dynamics in a korean pine forest in northeast China. PLoS One8, e82414. doi: 10.1371/journal.pone.0082414
63
ZhuJ.ZhuC.LuD.WangG. G.ZhengX.CaoJ.et al. (2021). Regeneration and succession: A 50-year gap dynamic in temperate secondary forests, Northeast China. For. Ecol. Manage.484, 118943. doi: 10.1016/j.foreco.2021.118943
Summary
Keywords
forest gap, natural regeneration, pattern and dynamic, primary forest, secondary forest
Citation
Wu D, Lu D, Zhu J, Ge X, Zhang J, Lin L, Wang X, Liu H and Zhang G (2023) Spatial and temporal regeneration patterns within gaps in the primary forests vs. secondary forests of Northeast China. Front. Plant Sci. 14:1305535. doi: 10.3389/fpls.2023.1305535
Received
02 October 2023
Accepted
13 November 2023
Published
28 November 2023
Volume
14 - 2023
Edited by
Ming Dong, Hangzhou Normal University, China
Reviewed by
Changxiao Li, Southwest University, China; David Pessanha Siqueira, New Zealand Forest Research Institute Limited (Scion), New Zealand
Updates
Copyright
© 2023 Wu, Lu, Zhu, Ge, Zhang, Lin, Wang, Liu and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jiaojun Zhu, jiaojunzhu@iae.ac.cn
†These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.