How the Strength of Monsoon Winds Shape Forest Dynamics

: Monitoring long-term forest dynamics can improve our understanding of how tree species change under varying environmental conditions over time. The Nanjenshan Nature Reserve, located in the southernmost part of Taiwan, is covered by natural tropical forest and affected by the northeast monsoon wind in the winter. To explore how monsoon wind shapes the structure of vegetation, we used the long-term data of the Lanjenchi plot (high wind exposure forest) and Nanjenshan plot I (low wind exposure forest) located in the Nanjenshan Nature Reserve. In these two plots, all free-standing trees with diameter at breast height (DBH) ≥ 1 cm were identiﬁed, measured, and mapped ﬁve times in thirty years (1989–2019). Temporal changes of tree abundance, dominance, recruitment, and mortality were calculated for all species. We recorded 9694 ± 3924 individual ha − 1 belonging to 126 species in the high wind exposure forest and 4019 ± 1791 individual ha − 1 belonging to 105 species in the low wind exposure forest in the latest census (census 5). From 1989 to 2019, stem density in the high wind exposure forests decreased (from 12,179 ± 6689 to 9694 ± 3924 individual ha − 1 ), while stem density in the low wind exposure forests increased (from 3253 ± 1208 to 4019 ± 1791 individual ha − 1 ). The stem density of the rare species, which survived only in the high wind exposure forest, decreased over time. However, the basal areas of these two forests barely changed during the study periods. Due to environmental conditions, the two forests have developed different forest structures in spaces and time. Species stem density in the high wind exposure forest and diversity could continually decrease and affect the maintenance of species diversity in the monsoon wind-shaped forest.


Introduction
Species are the basic units of a community that occurs together in space and time [1]. A species may change its relative abundance, immigrate, or go extinct locally [2][3][4][5]. Species community dynamics may be influenced by both exogenous forces such as climate and endogenous ecological forces [2,5]. How multiple species have survived under different environments through time can provide information about their potential for facing various environmental conditions in the future [5,6].
Long-term inventory data help to better understand forest structure dynamics and ecological processes [3]. Condit et al. [7] analyzed the long-term dynamic data from Pasoh and Barro Colorado Island (BCI) forest dynamic plots and found that the different environmental conditions and climate events cause BCI to have a higher mortality and recruitment rate than the Pasoh plot. Shen et al. [8] found that intraspecific competitive exclusion and density dependence caused the high mortality of species in the Dinghushan forest dynamic 2 of 12 plot. In addition to the relationships among species, environmental conditions may alter dynamic patterns over time. Vincent et al. [9] showed that the frequency of natural disturbance might cause higher mortality in New Guinea's forests than in other tropical forests. Condit et al. [3] showed that drought resulting from El Niño, caused a higher mortality rate for the species originally associated with relatively moist soils than for other species. Overall, long-term species survival is closely related to the environment [6,10,11].
The Lanjenchi plot and Nanjenshan plot I are located within the Nanjenshan Nature Reserve, which is part of the Kenting National Park and located on the east side of the mountain range near the coast of the Pacific Ocean [12]. In winter, the northeast monsoon wind is strong. It causes differing tree morphology (e.g., leaf thickness [13], leaf area [13]) and tree physiology (e.g., evapotranspiration rate [14]) between the trees in the high wind exposure and low wind exposure plots. Moreover, the species have different habitat associations, resulting in different vegetation compositions and structures under different wind exposure levels [15][16][17]. The Lanjenchi plot is a high wind exposed plot that was dominated by Illicium arborescens and composed of two Taiwanese vegetation associations, including Lithocarpus amygdalifolius-Illicium arborescens and Helicia formosana-Schefflera octophylla types [16]. The Nanjenshan plot I is a low wind-exposed plot composed by Dendrocnide meyeniana-Ficus benjamina forest dominated by Ficus benjamina [16]. These two plots represent the heterogeneity of environmental conditions under the northeast monsoon wind stress [15,16]. However, the patterns in which the species change through time are not well studied.
This study aims to understand how the woody vegetation composition differs in the plots having high wind exposure (Lanjenchi plot) versus low wind exposure (Nanjenshan plot I). We utilized the tree survey data from Lanjenchi plot and Nanjenshan plot I, which have been surveyed five times. The questions we addressed were as follows: (1) Do trees have different change in stem densities and basal areas between the Lanjenchi plot and Nanjenshan plot I? (2) How does species composition change over time between these two plots?

Study Site
The study site of this research is located at the Nanjenshan Nature Reserve, situated at the southernmost part of Taiwan ( Figure 1). It is a lowland tropical primary forest, without significant human activities. Based on the nearest automatic weather station (Bing-Lang station, elevation: 250 m; 120 • 50 12 E, 22 • 04 34 N) (1.4 km from Nanjenshan plot I, and 3.4 km from Lanjenchi plot) the mean annual precipitation is 3793 mm and mean annual temperature is 22.3 • C in the study region [18]. Inside the Nanjenshan Nature Reserve, are five forest dynamics plots, including the Lanjenchi (5.88-ha), Nanjenshan plot I (2.1-ha), Nanjenshan plot II (0.64-ha), Transect (0.76-ha), and Nanjen Lake (1.86-ha) [16]. We selected the two largest plots with the same number of survey intervals to represent plots affected by monsoon wind in a high wind exposure forest (Lanjenchi plot) and low wind exposure forest (Nanjenshan plot I) [9] (Figure 1). The average wind gust speed is significantly higher in the Lanjenchi plot than in Nanjenshan plot I [19]. The maximum wind velocity could be up to 20 m s −1 during the northeast monsoon wind season in the Lanjenchi plot [19,20]. The location and topographical environment of these two plots are likely to be the main cause of the different winter wind speed patterns [21]. The Lanjenchi plot (220 m by 240 to 300 m; 120 • 51 38 E, 22 • 03 23 N) is located on the coastal eastern upper slope of Mt. Wanlite. Elevation ranges from 280 to 340 m [16,18] and the slope is 17 • -27 • [16,18]. Its vegetation is classified as an Illicium-Cyclobalanopsis tropical winter monsoon forest [17]. The Nanjenshan plot I (150 m by 140 m; 120 • 50 51 E, 22 • 04 54 N) is located in the valley at the northwest foot of Mt. Nanjen. Its elevation ranges from 224 to 275 m [16,22] and its slope is 15 • -24 • [16,22]. The vegetation type is classified as a Dysoxylum-Machilus tropical mountain zonal foothill evergreen broad-leaved forest [17].  [16,22] and its slope is 15°-24° [16,22]. The vegetation type is classified as a Dys Machilus tropical mountain zonal foothill evergreen broad-leaved forest [17].

Data from the Study Site
This study used the long-term survey data of woody species from the Lanjen and Nanjenshan plot I. The Lanjenchi plot was originally a 3-ha rectangular plo lished in 1989 (the first survey was finished in 1991), but in 1997 was extended to ha, nearly rectangular plot. The plot was represented by surveying 588 10 by 10-m rats, and five surveys were completed in 1991 (census 1, only 3-ha), 1997 (census (census 3), 2013 (census 4), and 2019 (census 5). The Nanjenshan plot I was a 2.1tangular plot composed of 210 10 by 10-m quadrats. The Nanjenshan plot I was lished in 1993 and the five surveys were completed in 1993 (census 1), 2000 (cen 2008 (census 3), 2013 (census 4), and 2020 (census 5). All trees within a given quadr ing a diameter ≥ 1 cm were measured at a height of 130 cm, tagged, mapped, and fied to species level. We also investigated each 10 by 10-m quadrat's environmen tors, including elevation, slope, convexity, aspect, and exposure value. Chao et al. [ Ku et al. [18] described the collecting methods of each environmental factor.

Abundance and Dominance
Abundance (individual ha −1 ) was determined using the stem density of each tree. Dominance (m 2 ha −1 ) was determined using the basal area per hectare. IV (imp value, %) was calculated for relative dominance and abundance for each species census year, which indicated the relative importance of the species in the forest (I (relative abundance (%) + relative dominance (%))/2).
Generalized estimating equations were used to test the significance between fi sus years for each species. Tuckey'HSD was use for post hoc comparisons in over sus year. Statistical analysis was performed with R statistical software R v4.0.3 [23 the geeglm function in the geepack package [24]. Botanical nomenclature followed

Data from the Study Site
This study used the long-term survey data of woody species from the Lanjenchi plot and Nanjenshan plot I. The Lanjenchi plot was originally a 3-ha rectangular plot established in 1989 (the first survey was finished in 1991), but in 1997 was extended to a 5.88-ha, nearly rectangular plot. The plot was represented by surveying 588 10 by 10-m quadrats, and five surveys were completed in 1991 (census 1, only 3-ha), 1997 (census 2), 2005 (census 3), 2013 (census 4), and 2019 (census 5). The Nanjenshan plot I was a 2.1-ha rectangular plot composed of 210 10 by 10-m quadrats. The Nanjenshan plot I was established in 1993 and the five surveys were completed in 1993 (census 1), 2000 (census 2), 2008 (census 3), 2013 (census 4), and 2020 (census 5). All trees within a given quadrat having a diameter ≥ 1 cm were measured at a height of 130 cm, tagged, mapped, and identified to species level. We also investigated each 10 by 10-m quadrat's environmental factors, including elevation, slope, convexity, aspect, and exposure value. Chao et al. [16] and Ku et al. [18] described the collecting methods of each environmental factor.

Abundance and Dominance
Abundance (individual ha −1 ) was determined using the stem density of each tagged tree. Dominance (m 2 ha −1 ) was determined using the basal area per hectare. IV (importance value, %) was calculated for relative dominance and abundance for each species in each census year, which indicated the relative importance of the species in the forest (IV (%) = (relative abundance (%) + relative dominance (%))/2).
Generalized estimating equations were used to test the significance between five census years for each species. Tuckey'HSD was use for post hoc comparisons in overall census year. Statistical analysis was performed with R statistical software R v4.0.3 (R Core Team, Vienna, Austria) [23], using the geeglm function in the geepack package [24]. Botanical nomenclature followed Flora of Taiwan, 2nd edition [25].

DCA Analysis and the Environmental Factors
We used Detrended Correspondence Analysis (DCA, Hill and Gauch [26]) to understand how species were distributed in the Lanjenchi plot and Nanjenshan plot I. Hellinger transformation was conducted on the abundance data for each species in each quadrat. To better interpret the effect of the species distribution between the two forests, each quadrat's environmental factors were correlated with each DCA axis by using Pearson's correlation.
The DCA analysis and the correlation with environmental factors were conducted with R statistical software R v4.0.3 [23] by using the decorana function (for DCA) and envfit function (correlated with environmental factors) in the vegan package [27].

Diversity Index
We calculated species richness (S), Simpson's diversity index (D), and Shannon's diversity index (H') for each census year data in the Lanjenchi plot and Nanjenshan plot I. Species richness is simply the number of species in both plots. The Simpson's diversity index was calculated using the following equation: where S is the species richness, S i the ith species, and P i the proportion of the ith species in each community. Shannon's diversity index was calculated as the same parameter but with a different formula and weighting the abundance of rare species by natural logarithm: All diversity indexes were conducted with R statistical software R v4.0.3 [23] using the diversity function in the vegan package [27].

Dynamic Index
Demographic rates are defined as the rate constants of population growth, calculated for mortality from annual mortality rate (m, % yr −1 ) [7]: The recruitment rate (r, % yr −1 ) was calculated as [7]: Population growth (λ, % yr −1 ) is the difference between the r and m, the equation [7] is from: N 0 is the number of individuals in the previous census, N t the total number of individuals at the latest census of measurement interval t, and N st the number of individuals surviving from the previous census (t−1) to the latest census (t). The individuals that survived but did not meet the diameter at 1.3 m height threshold during the latest census were treated as the surviving individuals.

Forest Structure and Forest Composition in the Lanjenchi and Nanjenshan Plots
Forest structure and species composition differed between the two plots. The record in census 5 showed a stem density of 9694 ± 3924 individual ha −1 and basal area of  46.0 ± 26.3 m 2 ha −1 (average DBH of each tree is 6.91 cm) in the Lanjenchi plot and a stem density of 4019 ± 1791 individual ha −1 and basal area of 46.6 ± 43.0 m 2 ha −1 (average DBH of each tree is 10.76 cm) in the Nanjenshan plot I. A total of 160 species were recorded in the two plots in the census 5. There were 126 species in the Lanjenchi plot and 105 species in the Nanjenshan plot I, of which 71 species were shared in both plots. Species diversity in the Lanjenchi plot was higher than in the Nanjenshan plot I ( Table 1). The most important species (31.24%) in the Lanjenchi plot were Illicium arborescens, Castanopsis cuspidata var. carlesii, Schefflera octophylla, Ilex cochinchinensis, and Psychotria rubra. In the Nanjenshan plot I, the most important species were Psychotria rubra, Bischofia javanica, Dysoxylum hongkongense, Ficus benjamina, and Lasianthus obliquinervisare. DCA results showed that all 160 species have markedly different distributions within the Lanjenchi plot and the Nanjenshan plot I ( Figure 2). The first DCA axis explained 12.5% of the total variance in the data and the second DCA axis explained 4.87% of the variance. The left side of the first DCA axis was occupied by the species that were only distributed in the Lanjenchi plot (first DCA axis < −0.78), and the right side of the first DCA axis was occupied by species that were only found the in Nanjenshan plot I (first DCA axis > 0.925). The middle of the first DCA axis was occupied by the shared species of the Lanjenchi plot and Nanjenshan plot I. The shared species to the left side of the first DCA axis represent species having greater abundance in the Lanjenchi plot, and vice versa to the right side of the axis ( Figure 2). (average DBH of each tree is 10.76 cm) in the Nanjenshan plot I. A total of 160 species were recorded in the two plots in the census 5. There were 126 species in the Lanjenchi plot and 105 species in the Nanjenshan plot I, of which 71 species were shared in both plots. Species diversity in the Lanjenchi plot was higher than in the Nanjenshan plot I ( Table 1). The most important species (31.24%) in the Lanjenchi plot were Illicium arborescens, Castanopsis cuspidata var. carlesii, Schefflera octophylla, Ilex cochinchinensis, and Psychotria rubra. In the Nanjenshan plot I, the most important species were Psychotria rubra, Bischofia javanica, Dysoxylum hongkongense, Ficus benjamina, and Lasianthus obliquinervisare. The census 1 of Lanjenchi plot size was 3-ha, but was enlarged to 5.88-ha from census 2; Ϯ2. Shannon diversity index; Ϯ3. Simpson's diversity index; The small letter (a-d) indicates a significant difference between census by using Tuckey'HSD post hoc comparison.
DCA results showed that all 160 species have markedly different distributions within the Lanjenchi plot and the Nanjenshan plot I ( Figure 2). The first DCA axis explained 12.5% of the total variance in the data and the second DCA axis explained 4.87% of the variance. The left side of the first DCA axis was occupied by the species that were only distributed in the Lanjenchi plot (first DCA axis < −0.78), and the right side of the first DCA axis was occupied by species that were only found the in Nanjenshan plot I (first DCA axis > 0.925). The middle of the first DCA axis was occupied by the shared species of the Lanjenchi plot and Nanjenshan plot I. The shared species to the left side of the first DCA axis represent species having greater abundance in the Lanjenchi plot, and vice versa to the right side of the axis ( Figure 2).  DCA results showed that all 160 specie the Lanjenchi plot and the Nanjenshan pl 12.5% of the total variance in the data and variance. The left side of the first DCA axi distributed in the Lanjenchi plot (first DCA DCA axis was occupied by species that we DCA axis > 0.925). The middle of the first D of the Lanjenchi plot and Nanjenshan plot I DCA axis represent species having greater a to the right side of the axis (Figure 2).  Simpson's divers ference between census by using Tuckey DCA results showed that all 160 the Lanjenchi plot and the Nanjens 12.5% of the total variance in the da variance. The left side of the first D distributed in the Lanjenchi plot (fi DCA axis was occupied by species DCA axis > 0.925). The middle of th of the Lanjenchi plot and Nanjensha DCA axis represent species having g to the right side of the axis (Figure 2  DCA results showed that all 160 species have markedly different distributions within the Lanjenchi plot and the Nanjenshan plot I (Figure 2). The first DCA axis explained 12.5% of the total variance in the data and the second DCA axis explained 4.87% of the variance. The left side of the first DCA axis was occupied by the species that were only distributed in the Lanjenchi plot (first DCA axis < −0.78), and the right side of the first DCA axis was occupied by species that were only found the in Nanjenshan plot I (first DCA axis > 0.925). The middle of the first DCA axis was occupied by the shared species of the Lanjenchi plot and Nanjenshan plot I. The shared species to the left side of the first DCA axis represent species having greater abundance in the Lanjenchi plot, and vice versa to the right side of the axis (Figure 2).   The census 1 of Lanjenchi plot size was 3-ha, but was enlarged to 5.88-ha from c non diversity index; Ϯ3. Simpson's diversity index; The small letter (a-d) indicates a ference between census by using Tuckey'HSD post hoc comparison.
DCA results showed that all 160 species have markedly different distr the Lanjenchi plot and the Nanjenshan plot I (Figure 2). The first DCA 12.5% of the total variance in the data and the second DCA axis explaine variance. The left side of the first DCA axis was occupied by the species distributed in the Lanjenchi plot (first DCA axis < −0.78), and the right DCA axis was occupied by species that were only found the in Nanjens DCA axis > 0.925). The middle of the first DCA axis was occupied by the of the Lanjenchi plot and Nanjenshan plot I. The shared species to the left DCA axis represent species having greater abundance in the Lanjenchi plot to the right side of the axis (Figure 2).  The census 1 of Lanjenchi plot size was 3-ha, but was enlarged to 5.88-ha non diversity index; Ϯ3.
Simpson's diversity index; The small letter (a-d) indic ference between census by using Tuckey'HSD post hoc comparison.
DCA results showed that all 160 species have markedly different the Lanjenchi plot and the Nanjenshan plot I (Figure 2). The first 12.5% of the total variance in the data and the second DCA axis ex variance. The left side of the first DCA axis was occupied by the sp distributed in the Lanjenchi plot (first DCA axis < −0.78), and the r DCA axis was occupied by species that were only found the in Na DCA axis > 0.925). The middle of the first DCA axis was occupied b of the Lanjenchi plot and Nanjenshan plot I. The shared species to th DCA axis represent species having greater abundance in the Lanjench to the right side of the axis (Figure 2).

Simpson's diversity index; The small letter (a-d) indicates a significant difference between census by using Tuckey'HSD post hoc comparison.
DCA results showed that all 160 species have markedly different distributions within the Lanjenchi plot and the Nanjenshan plot I (Figure 2). The first DCA axis explained 12.5% of the total variance in the data and the second DCA axis explained 4.87% of the variance. The left side of the first DCA axis was occupied by the species that were only distributed in the Lanjenchi plot (first DCA axis < −0.78), and the right side of the first DCA axis was occupied by species that were only found the in Nanjenshan plot I (first DCA axis > 0.925). The middle of the first DCA axis was occupied by the shared species of the Lanjenchi plot and Nanjenshan plot I. The shared species to the left side of the first DCA axis represent species having greater abundance in the Lanjenchi plot, and vice versa to the right side of the axis (Figure 2).  DCA results showed that all 160 species have markedly different distributions within the Lanjenchi plot and the Nanjenshan plot I (Figure 2). The first DCA axis explained 12.5% of the total variance in the data and the second DCA axis explained 4.87% of the variance. The left side of the first DCA axis was occupied by the species that were only distributed in the Lanjenchi plot (first DCA axis < −0.78), and the right side of the first DCA axis was occupied by species that were only found the in Nanjenshan plot I (first DCA axis > 0.925). The middle of the first DCA axis was occupied by the shared species of the Lanjenchi plot and Nanjenshan plot I. The shared species to the left side of the first DCA axis represent species having greater abundance in the Lanjenchi plot, and vice versa to the right side of the axis (Figure 2).

Overall Dynamic Changes in the Plots
The overall stem density in the Lanjenchi plot decreased significantly from 12,179 ± 6689 individual ha −1 in census 1 to 9695 ± 3924 individual ha −1 in census 5 (Table 1). In the Nanjenshan plot I, stem density increased from 3253 ± 1208 individual ha −1 in census 1 to 4019 ± 1791 individual ha −1 in census 5. The overall basal area remained relatively stable in the Lanjenchi plot and the Nanjenshan plot I through time ( Table 1).
The recruitment rate was higher in the Nanjenshan plot I than in the Lanjenchi plot, except in census 5 when the recruitment rate of the Nanjenshan plot I was 2.31% yr −1 , while the Lanjenchi plot was 3.14% yr −1 ( Table 2). Table 2. Dynamic changes in the 5.88-ha Lanjenchi plot and 2.1-ha Nanjenshan plot I.

Census
Lanjenchi Plot Nanjenshan Plot I DCA results showed that all 160 species have markedly differen the Lanjenchi plot and the Nanjenshan plot I (Figure 2). The first 12.5% of the total variance in the data and the second DCA axis ex variance. The left side of the first DCA axis was occupied by the s distributed in the Lanjenchi plot (first DCA axis < −0.78), and the DCA axis was occupied by species that were only found the in Na DCA axis > 0.925). The middle of the first DCA axis was occupied of the Lanjenchi plot and Nanjenshan plot I. The shared species to t DCA axis represent species having greater abundance in the Lanjenc to the right side of the axis (Figure 2).  (averag tree is 10.76 cm) in the Nanjenshan plot I. A total of 160 species were reco plots in the census 5. There were 126 species in the Lanjenchi plot and 10 Nanjenshan plot I, of which 71 species were shared in both plots. Species Lanjenchi plot was higher than in the Nanjenshan plot I ( Table 1). The m species (31.24%) in the Lanjenchi plot were Illicium arborescens, Castanops carlesii, Schefflera octophylla, Ilex cochinchinensis, and Psychotria rubra. In t plot I, the most important species were Psychotria rubra, Bischofia javan hongkongense, Ficus benjamina, and Lasianthus obliquinervisare. The census 1 of Lanjenchi plot size was 3-ha, but was enlarged to 5.88-ha from c non diversity index; Ϯ3. Simpson's diversity index; The small letter (a-d) indicates a ference between census by using Tuckey'HSD post hoc comparison.
DCA results showed that all 160 species have markedly different distr the Lanjenchi plot and the Nanjenshan plot I (Figure 2). The first DCA 12.5% of the total variance in the data and the second DCA axis explaine variance. The left side of the first DCA axis was occupied by the species distributed in the Lanjenchi plot (first DCA axis < −0.78), and the right DCA axis was occupied by species that were only found the in Nanjens DCA axis > 0.925). The middle of the first DCA axis was occupied by the of the Lanjenchi plot and Nanjenshan plot I. The shared species to the left DCA axis represent species having greater abundance in the Lanjenchi plot to the right side of the axis (Figure 2). Figure 2. Distribution of 160 species between the first and second DCA axes. Solid the species that were only distributed in the Lanjenchi plot; hollow circles represen 1. The census 1 of Lanjenchi plot size was 3-ha, but was enlarged to 5.88-ha from census 2; r: recruitment rate; m: mortality rate; λ: population growth.
Mortality in the Lanjenchi plot was higher than in Nanjenshan plot I, except at census 4 ( Table 2), when mortality in the Nanjenshan plot I reached its highest level of 2.83% yr −1 . The highest mortality rate in the Lanjenchi plot was 2.86% yr −1 in census 5, followed by census 2, and census 4 had the lowest mortality rate (Table 2).
Population growth in the two plots showed different trends ( Table 2). The Lanjenchi plot showed the highest population growth in census 5 and the lowest population growth in census 3. In the Nanjenshan plot I, the highest population growth occurred in census 3 and the lowest population growth occurred in census 4.

Changes in Abundant and Dominant Species in the Lanjenchi and Nanjenshan Plots I
The species that only occurred in the Lanjenchi plot decreased the most in species composition. The importance values of Illicium arborescens, Eurya nitida var. nanjenshanensis, Ilex lonicerifolia var. matsudae, Syzygium buxifolium, and Symplocos theophrastifolia decreased the most (Figure 3). Litsea acutivena, and Elaeocarpus sylvestris increased in importance values and stem density (Figure 3).
The species that were only distributed in the Nanjenshan plot I showed a minor change in their importance values. The importance values of Aglaia elliptifolia increased, and Dendrocnide meyeniana and Dysoxylum hongkongense decreased the most (Figure 3). Few species had great change in their importance values in the Nanjenshan plot I (Table S2).
The shared species that survived in both Lanjenchi plot and Nanjenshan plot I showed more dynamics than the other two species distribution groups. Ficus benjamina, Schefflera octophylla, Turpinia ternate, and Lasianthus obliquinervis var. obliquinervis increased their importance values (Figure 3). These species are also more dominant in Nanjenshan plot I than in the Lanjenchi plot. Daphniphyllum glaucescens subsp. oldhamii, and Osmanthus marginatus had decreased importance values, these species were more abundant in the Lanjenchi plot than in the Nanjenshan plot I (Figure 3). Diversity 2022, 14, x FOR PEER REVIEW 7 of 12 Figure 3. Changes in stem density, basal area, and importance value of 30 species having the greatest increases and decreases by importance values from census 1 to census 5. Red bar indicates the species that only exist in the high wind exposure forest (Lanjenchi plot), the dark blue bar indicates the species that only exist in the low wind exposure forest (Nanjenshan plot I), and light purple bar indicates species that existing in both forests. *, **, and *** beside the bar representing the significant difference under the generalized estimating equations of α < 0.05, α < 0.01, α < 0.001, respectively from the preceding census.

Discussion
This study evaluated forest composition dynamic under different northeast monsoon wind stress levels in two nearby forest dynamic plots in the Nanjenshan Nature Reserve. Stem density decreased significantly in the high wind exposure value forest and increased in the low wind exposure forest. The results showed a more remarkable change in species composition in the high wind exposure value forest, wherein the importance of the most dominant species declined while species ranked in the mid-dominance range increased.

Different Wind Stressing Affects Forest Composition
Prevailing wind stress may cause the trees to grow into different shapes and sizes [28] (e.g., the windswept form [29]). Wind stress coupled with other environmental conditions creates a diverse habitat for different species to survive [15,16,18,19]. In the Nanjenshan Nature Reserve, forests at higher elevations and higher exposures to northeast monsoon winds have higher stem density than forests at lower elevations with lower exposure to northeast monsoon winds [16]. From Li et al. [17], the association of Lithocarpus amygdalifolius-Illicium arborescens is the result of the monsoon wind in Taiwan at the Lanjenchi plot. Stem density in the Lanjenchi plot is 2.5 times that of Nanjenshan plot I. Longterm high stem density forest may be the strategy that multiple trees can resist strong wind stress [30]. Trees have a different trade-off between morphology and physiology to resist the prevailing wind stress; trees become increasingly spindlier in the main stem with slender branches sparsely distributed to maximize light capture [30]. Moreover, crowded

Discussion
This study evaluated forest composition dynamic under different northeast monsoon wind stress levels in two nearby forest dynamic plots in the Nanjenshan Nature Reserve. Stem density decreased significantly in the high wind exposure value forest and increased in the low wind exposure forest. The results showed a more remarkable change in species composition in the high wind exposure value forest, wherein the importance of the most dominant species declined while species ranked in the mid-dominance range increased.

Different Wind Stressing Affects Forest Composition
Prevailing wind stress may cause the trees to grow into different shapes and sizes [28] (e.g., the windswept form [29]). Wind stress coupled with other environmental conditions creates a diverse habitat for different species to survive [15,16,18,19]. In the Nanjenshan Nature Reserve, forests at higher elevations and higher exposures to northeast monsoon winds have higher stem density than forests at lower elevations with lower exposure to northeast monsoon winds [16]. From Li et al. [17], the association of Lithocarpus amygdalifolius-Illicium arborescens is the result of the monsoon wind in Taiwan at the Lanjenchi plot. Stem density in the Lanjenchi plot is 2.5 times that of Nanjenshan plot I. Long-term high stem density forest may be the strategy that multiple trees can resist strong wind stress [30]. Trees have a different trade-off between morphology and physiology to resist the prevailing wind stress; trees become increasingly spindlier in the main stem with slender branches sparsely distributed to maximize light capture [30]. Moreover, crowded forests reduce the impacts of wind on single individual more than in a non-crowded forest [31]. Reducing tree height may also resist the effect of physical damage so that the height of the trees in a high wind exposure forest may be smaller than in a low wind exposure wind forest [28]. In the Lanjenchi plot, the exposure value of the northeast monsoon wind is the determinate factor that causes species to have different associations with different habitats. For the Nanjenshan plot I, the environment is wind-shaded and more homogeneous than the Lanjenchi plot, causing their vegetation composition to be relatively similar [15,16,22]. Overall, the northeast monsoon wind and topography cause the higher stem density and smaller trees diameter in the high wind exposure forest of the Lanjenchi plot than in the low wind exposure forest of the Nanjenshan plot I. Although the total basal area in the Lanjenchi plot was higher than in the Nanjenshan plot I, it means that numerous small trees in the Lanjenchi plot had comparable biomass with fewer larger trees in Nanjenshan plot I.

From Species to Overall Change in Dynamics
Any increase in some species must be accompanied by a matching decrease in the collective numbers of all other species in the community [32]. Species composition and structure have different distribution between the high wind exposure forest and low wind exposure forest and show the different dynamics. The following will discuss species changes over time within the three major groups: species only found in the high wind exposure forest, species only found in the low wind exposure forest, and species found in both forests.
The species that were only found in the high wind exposure forest decreased their stem density, basal area, and importance values over time, but tended to survive at higher elevations and under higher exposure values. Species including Illicium arborescens, Eurya nitida var. nanjenshanensis, Ilex lonicerifolia var. matsudae, Syzygium buxifolium, Osmanthus marginatus, Symplocos theophrastifolia, and Symplocos congesta were the most abundant species in census 1 (Supplementary Material Table S1); however, these exhibited large decrease in stem density, basal area, and importance value over time. These species, except for Symplocos congesta, and Symplocos theophrastifolia, had a strong species-habitat association with the windward habitat and were designated as windward distributed species in Ku et al. [18]. Windward distributed species specialize in small, thick, and waxy leaves to prevent rapid evapotranspiration under the strong wind stress conditions in winter [12,13]. Thus, those specialized functions are an adaptation to environmental stress, but in spite of these adaptations, their stem density, basal area, and importance values decreased over time. Moreover, eight out of 22 windward distributed species are contained in the red list of vascular plants of Taiwan [33], including critically endangered (Lithocarpus formosanus), endangered (Nageia nagi and Podocarpus macrophyllus), vulnerable (Camellia hengchunensis and Symplocos shilanensis), and near-threatened species (Eurya nitida var. nanjenshanensis, Myrsine seguinii, and Ilex lonicerifolia var. matsudae) [33]. Six of the eight endangered species showed a decrease in stem density, basal area, and importance values ( Figure 3, Table S1). Among these species, the distribution of Lithocarpus formosanus was only recorded in the Nanjenshan Nature Reserve. Currently, there is no apparent cause for those endangered species to undergo the observed decreasing patterns, but in Ku et al. [21], large trees tend to have a high growth rate, and small trees tend to self-thin. Further examination of the dynamics of the size distribution of these species is needed to better understand how the red-listed species survive over time. Nonetheless, conservation plans for these endangered species in the wind-affected forest are needed to ensure their persistence in the future. Among the species known to exist only in high wind exposure forests, Melastoma candidum increased in stem density significantly over the study period. Melastoma candidum is a fast-growing tree whose photosynthetic capacity can range between 24.0-25.9 µmol CO 2 m −2 s −1 [34]. The increase of this species could be due to the Nanjenshan Nature Reserve being disturbed by a typhoon in 2016, causing defoliation in the forest and felling Diversity 2022, 14, 169 9 of 12 numerous large trees. The created gaps could have provided a suitable regenerative habitat for this species [35,36].
The species found in both high and low wind exposure plots show a positive change in stem density, basal area, and importance value. These species are not the most dominant species in the forest, but increased their basal area over time and became the dominant species. For examples, Schefflera octophylla, Schima superba var. kankaoensis, Ilex uraiensis, and Litsea acutivena were the larger diameter species in the forest, and they accumulated more aboveground biomass through time, showing that the importance of the large tree species is crucial to maintaining biomass in the forest [37,38].
The species that were only found in the low wind exposure forest tended to be bigger than those in the high wind exposure forest. Ficus benjamina dominated the most space with its large diameter and numerous prop roots. Other species in the Nanjenshan plot I showed an increasing trend in importance values, in contrast to the dynamic of trees in the wind-stressed forest. The originally dominant species were no longer dominant over time, and the species that were originally lower in the ranking became more prominent. Legendre and Condit [5] showed that the drought from El Niño caused the stem density to decrease and changed the beta diversity difference in different habitats at Barro Colorado Island forest dynamics plot. Species from other habitats invaded BCI swamp habitat due to drought [5].

How the Monsoon Wind Shaped the Forest of the Future
Fluctuations in environmental stress may influence species changes in density or biomass. Based on the stress gradient hypothesis, high environmental stress release may transfer species interaction from facilitation to competition [39]. Grassland research in southern Norway showed that global warming increases competition in dominant species [40]. Some alpine research also found that the intensity of stress may affect the relationship between biomass and density [41,42]. Crabtree and Ellis [43] revealed that decreasing wind speed induced by global warming caused the stress-tolerant lichens to compete for space and light in the Cairgorm Mountains. The stress gradient hypothesis has been tested in herb-, shrub-, and small tree-dominated ecosystems [40,44], but seldom in tropical forests. At our study site, the northeast monsoon wind is an important influential factor affecting two different forest compositions [15,18,22]. According to the research from Hung and Kao [45] and Wang and Lu [46], global warming changes the pattern of the Siberian high and causes the East Asian winter monsoon to decrease. Based on the stress gradient hypothesis, once the environmental stress decreases, species interaction in the forest may change. This hypothesis may explain why the stem density decreased in windward distributed species but increased in other species. More field survey work and experiments are needed to investigate this hypothesis.

Conclusions
This study describes species changes in two forests having different levels of wind exposure at the Nanjenshan Nature Reserve, Taiwan. High wind exposed forests had high stem density and small trees to defend against wind stress in winter. However, stem density decreased significantly over time in the most abundant and endangered species. In contrast, the low wind exposure forest had a low stem density, and its large trees increased their basal area to dominate the forest over time. The high wind exposure forest brought about more negative changes in the tree community than the low wind exposure forest. Species fluctuation in the Nanjenshan Nature Reserve caused the dominant species to no longer be dominant.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/d14030169/s1. Table S1: Change of the density, basal area, and importance value of each species under different census times for the species that only distributed in the Lanjenchi plot. Table S2. Change of the density, basal area, and importance value of each species under different census times for the species that only distributed in the Nanjenshan plot 1. Table S3. Change of the density, basal area, and importance value of each species under different census times for the species that distributed in both Lanjenchi plot and the Nanjenshan plot 1.

Institutional Review Board Statement: Not applicable.
Data Availability Statement: All of the raw data in each census are belonging to the funding unit. The data can be derived from the site principal investigator. The information of the principal investigator in Nanjenshan site can be found on the ForestGEO website (https://forestgeo.si.edu/ sites/asia/nanjenshan (accessed on 24 February 2022)).