Variation in Niche and Interspecific Associations across Elevations in Subtropical Forest Communities of the Wuyi Mountains, Southeastern China

: Elucidating changes in the structure and function of plant communities along different elevation ranges will help researchers to analyze the strategies plant communities use in environments and processes influencing niche differentiation. The aims of this paper are to reveal the underlying mechanisms and ecological processes governing the development of subtropical forest ecosystem plant communities. This paper analyzes the forest vegetation of the Wuyi Mountains across the following three elevation ranges: low elevation, mid elevation, and high elevation, spanning from 560 to 2150 m. Twenty and twenty-three dominant tree and shrub layer species, respectively, were identified based on their importance values, and their niches and species associations were further analyzed based on the elevation range. The results showed interspecific associations between tree and shrub species, with the strongest associations observed at mid-elevations. The analysis of niche width and overlap showed that the number of pairs of species with a higher degree of niche overlap decreased with increasing elevation, suggesting that resource use varied at different elevations for both tree and shrub layer species, which may be related to the adaptive capacity of plants at different elevations to the environment and resource use strategies. These findings should contribute to a deeper understanding of the ecological functioning and structural framework of plant communities on Wuyi Mountain.


Introduction
During forest community assembly, different plant species show various kinds of interdependence and associations; their abundance and diversity change significantly over time and space.These changes affect the stability of forest ecosystems [1][2][3].The niche of a species reflects its ability to use environmental resources and interact with other species [4].Niche is therefore an important driver of species coexistence and assembly change in plant communities [5].Interspecific association provides a static description of the interrelationships between the species that make up a plant community over a certain period of time and represents a qualitative association based on binary data that tests whether two species are related.The convergence of species to the environment and their interaction during community assembly frequently give rise to interspecific associations.This statement suggests that these entities play a crucial role not only in the establishment and progression of a community but also in demarcating the said community and describing how different species interact with each other [6].Research on plant niches and interspecific associations can aid in identifying associations between plant species, predicting population changes, and revealing mechanisms of plant substitution during community assembly.This research holds valuable theoretical and practical applications, including applications related to vegetation recovery and restoration, environmental and biodiversity conservation, preventing the creation of degraded ecosystems and managing those that are created, and the implementation of related policies and regulations.
Various environmental factors play crucial roles in shaping plant niches and interspecific associations [7,8].For example, in the Amazon region, species niche breadth increases proportionally with their niche position across different soil nutrient gradients [9].Demographic trade-offs may emerge along edaphic gradients and consequently aid in species sorting among soil-related habitats in Bornean rainforests [10].Additionally, soil nutrients and water quality exert a significant impact on interspecies associations among mangroves in East Lampung, Indonesia [11].These studies illustrate that niches and associations between species undergo changes along environmental ranges.In mountainous ecosystems, environmental factors such as vegetation, climate, and soil undergo significant changes with increasing elevation, presenting an ideal platform for investigating the distribution patterns of biodiversity along elevation ranges and their responses to changes in the environment [12,13].By examining the alterations in plant ecological niches and interspecific associations across elevation ranges, we can gain insights into the mechanisms by which species distributions and community assemblages respond to environmental changes.This, in turn, allows us to develop more effective conservation and management strategies [14].
Subtropical mountain forests boast rich plant diversity [15], which contributes to the formation of several global biodiversity hotspots [16].These areas are ideal for studying biodiversity patterns and the biotic and abiotic drivers of community assembly [17].Current research has examined how some subtropical plant communities respond to environmental ranges, but the underlying mechanisms and theories remain controversial.Kröber [18] highlighted a stronger correlation between the environment and species in Chinese subtropical forests than between the species themselves.Pei [19] also asserted that niches play a significant role in plant community aggregation across different habitats within the subtropical region.However, Umaña [20] found that some growth-related traits, such as leaf area or wood density, exhibited minimal variation along elevation ranges in the subtropics.This finding may suggest that niches and interspecific associations of certain central subtropical plants do not change significantly with elevation.Nonetheless, further exploration and research are necessary to gain a better understanding of how plant communities change with elevation in a subtropical region.
The Wuyi Mountains preserved relatively complete primitive forests [21].This forest ecosystem exhibits extensive plant stratification, with the dominant species at each level being prominent and changing significantly with an increase in elevation.This type of environment can be studied to effectively characterize the response of plants to environmental factors, including hydrothermal gradients.Additionally, the mountains offer a suitable location to investigate plant communities and conduct related research in the vertical zone spectrum of the middle subtropic region [22].This setting provides a natural site to analyze the vertical zonation spectrum of plant communities in a typical central subtropical region.Currently, there is extensive research on vegetation productivity [23][24][25], forest stoichiometry characteristics [26][27][28], soil microbial community structure [29], and plant economic spectrum [30] across different elevations in the Wuyi Mountains.Nevertheless, studies on interspecific associations within plant communities and their responses to environmental changes along different elevation ranges are relatively scarce.
The present study involved conducting a species survey of the forest ecosystem on Wuyi Mountain, which was originally a natural reserve area.The interspecific associations and niches of tree and shrub layer plants at various elevation ranges were analyzed.This study aimed to answer the following questions: (1) What are the characteristics of niches and interspecific associations in plant communities on Wuyi Mountain at different elevations?(2) How do changes in niches and interspecific associations in this region respond to ecosystem responses to environmental change?The final objective of this study is to establish a preliminary understanding of the intricate associations between ecosystems and the environment.The insights obtained will provide a solid research foundation for the development of scientifically informed conservation strategies for Wuyi Mountain.

The Study Area
The Wuyi Mountain National Nature Reserve is located in northwestern Fujian Province at 27 The reserve covers a total area of 100,000 hectares, making it the largest and most complete subtropical primary forest ecosystem within the subtropical latitude zone globally.The region features abundant natural resources and is a crucial area for global biodiversity conservation [31].The typical subtropical monsoon climate includes abundant rainfall, reaching up to 2152 mm annually.The annual average evaporation is 998 mm, while the relative humidity ranges from 78% to 84%.The temperature ranges from −15 • C to 18 • C. The frost-free period lasts 252-257 days per year.Additionally, the Huanggang Mountain, the highest peak in the area, is not only the main peak of the Wuyi Mountain range but is also the highest peak in Southeastern China, with a vertical increase in elevation of 450-2158 m above mean sea level.The southeast slope is relatively gentle, while the northwest slope is steep.The evergreen broad-leaved forest in the area is rich in plant species, with notable dominant species.The differences in the dominance of the dominant species in different parts of this region are significant, while the distribution of the vertical vegetation belt spectrum is complete.The primary forest types along the elevation range from low to high elevation are evergreen broad-leaved forest, coniferous forest, and subalpine dwarf forest.The main dominant tree species include Castanopsis carlesii (Hemsl.)Hayata., Quercus glauca Thunb., Pinus massoniana Lamb., Cunninghamia lanceolata (Lamb.)Hook., Pinus taiwanensis Hayata, and Taxus wallichiana var.mairei (Lemée and H. Lév.) L. K. Fu and Nan Li, which are warm or temperate tree species [32,33].

Plot Establishment and Data Collection
From May to August 2015, 34 plots measuring 20 m × 30 m each were established in the Wuyi Mountain range and were located at elevations ranging from 560 m to 2150 m in areas where the forest structure was intact (Figure 1).The vertical spacing between plots ranged from 20 to 100 m (mostly 50 m).In addition, five 2 m × 2 m shrub layer subplots were established at the center and corners of each plot.A total of 34 quadrats in the tree layer (20 m × 30 m) and 170 quadrats in the shrub layer (2 m × 2 m) were established.The selected data included species name, number of individuals in a plot, coverage, diameter at breast height, tree height, crown height, crown width, and base diameter.The woody plants with a diameter at breast height (DBH) less than or equal to 5 cm at 1.3 m in the shrub sample plots were identified as shrub layer plants, and the woody plants with a DBH greater than 5 cm were identified as tree layer plants [34,35].The shrub layer characteristics that were recorded included species name, number of individuals in a plot, coverage, height, crown width, and basal diameter (cluster diameter).Moreover, the geographical coordinates, elevation, slope, aspect, forest type, and other environmental characteristics of each plot were recorded.Sample plots 3 and 4 were part of a subalpine meadow that consisted of shrubs and had no trees.
In this study, a total of 206 tree species from 103 genera and 56 families were surveyed.Dominant families in the tree layer included Theaceae, Fagaceae, and Lauraceae.Additionally, 163 shrub species from 92 genera and 49 families were surveyed, with the dominant families being Theaceae, Rosaceae, and Ericaceae.Among these, 20 tree species and 23 shrub species had importance values exceeding 1%.The combined importance values of these 20 tree species and 23 shrub species accounted for 81.07%and 53.81% of the total importance values in the tree and shrub layers, respectively.Consequently, these species were selected for subsequent niche and association analysis; the survey results are presented in Table 1.In this study, a total of 206 tree species from 103 genera and 56 families were surveyed.Dominant families in the tree layer included Theaceae, Fagaceae, and Lauraceae.Additionally, 163 shrub species from 92 genera and 49 families were surveyed, with the dominant families being Theaceae, Rosaceae, and Ericaceae.Among these, 20 tree species  The plots were classified into three types (Table 2) using multiple regression trees, with elevation as independent variables and the species and quantity of trees and shrubs within each quadrate as dependent variables [36][37][38].These plots were divided into three elevation ranges, and they were named H1 (low-elevation range), H2 (medium-elevation range), and H3 (high-elevation range).The 10 H1 plots were located in low-elevation habitats, each featuring evergreen broad-leaved forests.The 10 H2 plots at the middle elevation had broad-leaved mixed forests, whereas the 14 H3 plots were located at high elevations with coniferous forests.The main species of each plot were determined according to the species importance values, as follows [39]: where I V , R a , R f , and R d represent the importance value, relative density, relative frequency, and relative dominance of a specific species, respectively.a i is the number of individuals in population i, f i is the number of quadrats in which the population i appears, d i is the basal area at the height of 1.3 m of population i, and S is the total number of species.

Determination of Interspecific Association
The quantitative study of interspecific associations originated from the pioneering of the discipline of probability theory and mathematical statistics at the beginning of the 20th century.Different calculation methods were proposed successively, but all of them used the 2 × 2 association table (Table 3) to calculate the occurrence of species in a specific plot [40].The association coefficient (AC) was used as the basis for the χ 2 statistic to determine interspecific associations [41,42].First, the overall connectivity of tree and shrub species within the plant community of each plot in both the tree and shrub layers was tested based on the variance ratio (VR).The significance of VR was assessed using the W statistic test, as follows:

Table 3. Measurement of the interspecific associations between species A and B.
where VR represents the overall species association, N is the total number of sample plots, T j is the number of species in a specific plot j, t is the average species count per plot, S is the total species count, and P i is the number of plots where species i is present.Under the null hypothesis of independence (VR = 1), all species are considered unrelated.Therefore, VR > 1 indicates a positive association, while VR < 1 indicates a negative association.In addition, W represents the W statistic, VR is the total species association number, and N is the total number of sample plots.If χ 2 0.95(N) < W < χ 2 0.05(N) , this indicates that the overall species association is insignificant (p > 0.05).Otherwise, it signifies a significant association (p < 0.05) [43].
A χ 2 test was used to determine whether any two species were related.Since the sampling method was non-continuous and the original data consisted of binary values indicating the presence or absence of events, the χ 2 value for non-continuous data was calculated using Yates' correction formula, as follows [44]: where n is the number of sampling species.When χ 2 < 3.841, the interspecific association is independent; when χ 2 > 6.635, a significant ecological association exists between any two species; and when 3.841 < χ 2 < 6.635, there is a certain ecological association between species.Table 3 [44].Following qualitative interspecific associativity using the χ 2 test, AC is commonly used to further determine the degree of associativity of species pairs [45,46], which are calculated as follows: If .The closer the value of AC is to 1, the greater the positive connection between any two species is.Conversely, the closer the AC value is to −1, the stronger the negative association between any two species is.When the value of AC is 0, the species are entirely independent [47,48].

Calculation of Niche Width and Niche Overlap
The niche width of tree and shrub species in the sample plots was calculated using the Levins' index, as follows [49]: where B i represents the niche width value of species i, which falls within the range of [0-1], with larger values indicating a broader niche for the species; P 2 ij represents the proportion of the i th species' resource use relative to all resources of type j, and r represents the number of resource types, i.e., the total number of sample plots.
The niche overlap values of tree and shrub species in the sample plots were calculated using the niche overlap index, as follows [50]: where O ik is the niche overlap index between species i and k; and P ij and P kj are the proportion of resources j used by species i or k, respectively.In this study, niche and interspecific associations were computed using the "spaa" and "corrplot" packages in R 4.3.1.The plots were generated using R 4.3.1 and OriginPro 2024 [51].

Variations in Interspecific Associations with the Elevation Gradient
The overall interspecific associations between the tree and shrub layers within the Wuyi Mountains are detailed in Table 4.Most plot communities exhibited VR values greater than 1.00, with the exception of the shrub layer community in the H3 plot.The tree and shrub layers in H1 showed significant positive associations, as did the shrub layer in H2 and the tree layer in H3.Meanwhile, the tree layer in H2 showed an insignificant positive association, and the shrub layer in H3 exhibited an insignificant negative association.The results of the Chi-squared and the W tests are in agreement, indicating a lack of significant associations among the majority of species.In the tree layer, the locations of the prevalence of species pairs exhibiting significant associations could be ranked as follows: high-elevation > low-elevation > mid-elevation ranges.In contrast, the proportion of species pairs showing a significant association declined with increasing elevation within the shrub layer.The peaks of species pairs exhibiting highly significant associations in both the tree and shrub layers occurred in the high-elevation range at 16% and 4%, respectively.
The interspecific associations between the tree and shrub species in each plot indicated by the AC values (Figure 2) showed that overall, more species pairs in the tree layer exhibited positive associations than negative associations.The proportion of species pairs in the tree layer exhibiting significant positive associations at all three elevations was relatively high.In the mid-elevation range, the proportions of species pairs with both significant positive and negative associations in the tree layer were the highest, accounting for 35.71% and 33.77% of the total species pairs at that elevation, respectively.However, the proportions of species pairs with significant positive and negative associations in the shrub layer increased as the elevation rose.H2, and (f) H3 of the shrub layer.▲ extreme significant positive association; • significant positive association; □ + non-significant positive association; □ non-significant negative association; ○ significant negative association; ▽ extreme significant negative association.The scientific names of trees and shrubs are shown in Table 1.

Niche Width Analysis
The niche width of various species calculated using the Levins's index (Figure 3) showed that tree species with relatively larger niche widths included the following: in the low-elevation range, Castanopsis eyrei, Quercus glauca, Vaccinium sprengelii, and Rhododendron ovatum; in the mid-elevation range, Rhododendron latoucheae, Pinus taiwanensis, Vaccinium sprengelii, and Schima superba; at high-elevation range, Pinus taiwanensis, Enkian-  H2, and (f) H3 of the shrub layer.▲ extreme significant positive association; • significant positive association; □ + non-significant positive association; □ non-significant negative association; ○ significant negative association; ▽ extreme significant negative association.The scientific names of trees and shrubs are shown in Table 1.

Niche Width Analysis
The niche width of various species calculated using the Levins's index (Figure 3) showed that tree species with relatively larger niche widths included the following: in the low-elevation range, Castanopsis eyrei, Quercus glauca, Vaccinium sprengelii, and Rhododendron ovatum; in the mid-elevation range, Rhododendron latoucheae, Pinus taiwanensis, Vaccinium sprengelii, and Schima superba; at high-elevation range, Pinus taiwanensis, Enkianthus quinqueflorus, Rhododendron latoucheae, and Quercus multinervis.Detailed data are shown in Appendix A.
non-significant positive association; □ non-significant negative association; significant negative association; ▽ extreme significant negative association.The scientific names of trees and shrubs are shown in Table 1.

Niche Width Analysis
The niche width of various species calculated using the Levins's index (Figure 3) showed that tree species with relatively larger niche widths included the following: in the low-elevation range, Castanopsis eyrei, Quercus glauca, Vaccinium sprengelii, and Rhododendron ovatum; in the mid-elevation range, Rhododendron latoucheae, Pinus taiwanensis, Vaccinium sprengelii, and Schima superba; at high-elevation range, Pinus taiwanensis, Enkianthus quinqueflorus, Rhododendron latoucheae, and Quercus multinervis.Detailed data are shown in Appendix A.  In general, the width of the niche of the Ericaceae was the broadest in both the tree and shrub layers.The Fagaceae demonstrated greater dominance in the tree layer at lower elevations, while the Pinaceae and Ericaceae gradually increased their dominance as the elevation increased.Within the shrub layer, species composition was more complex at lower elevations, with larger niche widths for several families, while Ericaceae species had particularly complex species composition at middle and high elevations.
The niche overlap index analysis (Figure 4) showed that in the tree layer, the niche overlap values of Rhododendron ovatum and Vaccinium sprengelii with several dominant species at low elevations reached 0.9, which may suggest intense competition with other species.In the mid-elevation range, Rhododendron latoucheae and Pinus taiwanensis, as  1.Each color block in the outer ring represents a species or an elevation range.The direction of the lines illustrates the dominant species presented in different elevations in these samples, and the thickness of the bar is proportional to the importance value of the dominant species of each sample.
In general, the width of the niche of the Ericaceae was the broadest in both the tree and shrub layers.The Fagaceae demonstrated greater dominance in the tree layer at lower elevations, while the Pinaceae and Ericaceae gradually increased their dominance as the elevation increased.Within the shrub layer, species composition was more complex at lower elevations, with larger niche widths for several families, while Ericaceae species had particularly complex species composition at middle and high elevations.
The niche overlap index analysis (Figure 4) showed that in the tree layer, the niche overlap values of Rhododendron ovatum and Vaccinium sprengelii with several dominant species at low elevations reached 0.9, which may suggest intense competition with other species.In the mid-elevation range, Rhododendron latoucheae and Pinus taiwanensis, as well as in the high-elevation range, Pinus taiwanensis and Enkianthus quinqueflorus, exhibited higher overlap values with other dominant species.1.
In the shrub layer, Begonia acetosella competes with most other dominant species at both low and mid elevations.In the high-elevation range, most shrub layer-dominant species show lower niche overlap values with each other, suggesting less intense competition among shrub species in high-elevation ranges.1.
In the shrub layer, Begonia acetosella competes with most other dominant species at both low and mid elevations.In the high-elevation range, most shrub layer-dominant species show lower niche overlap values with each other, suggesting less intense competition among shrub species in high-elevation ranges.

Regulation Mechanisms of Interspecific Associations Along Environmental Ranges
Interspecific associations are an important part of the structural and quantitative characteristics of a community [52].Studying plant interspecific associations allows researchers to understand the interactions between different species and the way various plants adapt to the environment [53].The results of the present study show that the number of significantly positively correlated species pairs among dominant species in tree layers at all three different elevations exceeded that of negatively correlated species pairs.However, most dominant species lacked associations with each other.Consistent findings were obtained in some previous studies, as follows: Jin [53] found that interspecific associations among plants were not apparent in the region of the Taihang Mountains.In comparison to the low-elevation range, the number of species pairs exhibiting positive or negative associations decreased at middle and high-elevation ranges.Dai [54] observed that interspecific associations among different tree species at varying elevations in the forest-grassland ecotone were not significant.They pointed out that there was no evidence of exclusion caused by interspecific competition within the local plant community.Nguyen [55] also suggested that the association patterns among dominant tree species were largely independent.This indicates that the tree communities at different elevations in the Wuyi Mountains are in relatively stable successional stages, with community structures approaching completeness, yet maintaining certain mutualistic associations [2].
The total number of negative connectivity pairs among dominant species in the shrub layer surpassed that of positive connectivity pairs in the present study.Additionally, as the elevation increased, the proportion of shrub species pairs exhibiting significantly positive and negative associations gradually increased.This trend may be caused by the influence of tree layer plants on shrub layer plants.The competition for resources between the shrub and tree layer plants was asymmetric, with the tree layer exerting pressure on the understory through factors such as sunlight, rainfall, and soil nutrients, which ultimately led to a reduction in competitive traits in the understory [45].Our study observed a decrease in the number of tall broadleaf trees and an increase in the number of coniferous trees at medium elevations, while tall trees were sparser at high elevations than at low elevations.This indicates a reduction in competitive pressures from middle-to high-elevation tree layers, thereby enabling shrub layer plants to gain access to a greater abundance of resources, such as sunlight and rainfall.This could intensify competition between species in the shrub layer, resulting in more species pairs showing negative associations.
Compared to low-and high-elevation sites, the mid-elevation range stands out due to its more pronounced characteristics.In this zone, a higher proportion of species pairs was observed that show significant positive and negative associations among tree layer plant communities than at other elevations.Additionally, both tree and shrub layer species richness and diversity peaked in the mid-elevation range.These findings align with previous research; for instance, in the southern part of the Taihang Mountains [2], at lower and higher elevations, both tree and shrub layer species richness and diversity exhibited an overall non-significant negative interspecific association, while at moderate elevations, this type of species richness and diversity demonstrated a non-significant positive association.According to Tetemke [56] and Bruun [57], the association between species richness and elevation does not follow a linear pattern but rather peaks at moderate elevations.One possible explanation for the coexistence of a higher number of species at the mid-elevation is the greater availability of resources [58,59].High-elevation ranges are distinguished by a number of factors that are generally considered to be unsuitable for plant growth.These include low temperatures, short growing seasons, unstable substrates, intense radiation, and strong winds [60].Consequently, in alpine habitats with more extreme environments, plant communities are highly susceptible to any environmental changes, rendering this environment one of the least suitable habitats for plant growth and development [61].Conversely, lower elevation ranges are characterized by warmer and moister conditions, which are conducive to plant survival.However, these ranges are dominated by species such as Quercus glauca and Castanopsis eyrei, which exert strong limitations on the growth of other species [57].This results in lower species richness and interspecific associations compared to middle elevation ranges.This phenomenon is reflected by the indicators of species diversity and higher biodiversity in the mid-elevation range, which point to a more stable and balanced community structure in the range [62].Different species maximize resource use through niche differentiation [63,64], leading to stronger overall interspecific associations in the community [65].

Response of Niche Overlap to Environmental Ranges
Niche width primarily reflects the extent to which species use resources, with a broader niche indicating a higher degree of resource use [66,67].It is generally believed that plants capable of using a variety of resources have a greater likelihood of surviving across different environmental gradients and coexisting with a greater number of species [68].The results showed that several dominant species consistently exhibited wide niches at various elevations within the study sites.Examples include Rhododendron latoucheae, Vaccinium sprengelii, Schima superba, and Rhododendron ovatum in the tree layer, as well as Begonia acetosella, Tauntonia obovatifoliolata, Itea omeiensis, and Rhododendron simsii in the shrub layer.These plants were well-suited for growth in environments with moderate sunlight exposure while avoiding direct sunlight.Consequently, they could effectively use spatial resources beneath the canopy of evergreen broad-leaved forests or mixed needle-leaf and broadleaved forests.Moreover, their smaller size resulted in lower competition pressure from neighboring species, allowing them to coexist with numerous species while possessing strong habitat utilization capabilities [69].This may explain their broader niche width at different elevations.Additionally, it is noteworthy that the niche width of the most dominant species within the study sites peaked at mid-elevation (11 tree species and 12 shrub species), confirming the previous idea that mid-elevation sites can accommodate more species.
Niche overlap serves as a reflection of the similarity in resource use and the degree of spatial position occupancy among different populations [70].A higher degree of niche overlap indicates a greater similarity in the use of natural resources among species, potentially leading to more intense interspecific competition [71].As the elevation increases, the proportion of species pairs in the Wuyi Mountain study sites in both the tree and shrub layers, with niche overlap values greater than 0.5, gradually decreases.This reduction was particularly pronounced in the high-elevation ranges, indicating significant differences in the environmental resource requirements among species at higher elevations compared to mid and low elevations.
Previous studies have suggested that the effects of competitive actions are diminished at high elevations or in harsh environments [72].However, recent experimental results have challenged this perspective, indicating that competition still plays a crucial role in species coexistence in high-elevation environments [73].The present study showed that tree layer plants with broader niche widths at various elevations generally exhibited varying degrees of niche overlap with the majority of dominant species in that range.Moreover, noticeable changes in niche overlap values were observed between the same species pairs as elevation increased, which aligns with their adaptability to diverse environments.Simultaneously, the divergent trends in the changes in niche overlap values among different species pairs imply that competition plays a role in the development of plant communities at different elevations.
It is important to note that niche width is not necessarily correlated with niche overlap.This is a result of the spatial heterogeneity of resources used by different species and their varying abilities to use them.Even species with narrow niche widths may exhibit a high level of niche overlap [74].Our results support this viewpoint, as we observed instances of high niche overlap among species with narrow niche widths at various elevations.Species with narrow niches may share similar resource and habitat requirements, potentially leading to intense competition.

Relationship between Niche Overlap and Interspecific Associations
Interspecific associations and the degree of niche overlap are also related [75].Species with significantly positive associations usually have similar environmental and resource requirements, leading to a relatively high degree of niche overlap values [76].The present study suggested that most species pairs in the study sites exhibited similar trends in their interspecific associations and niche overlap values even though the elevational changed.However, not all species pairs showed a correlation between interspecific associations and the degree of niche overlap.Some species pairs in the study sites had a high degree of niche overlap values but weak interspecific associations, while others species pairs exhibited significant interspecific associations without niche overlap.This discrepancy may be related to the scale of this study.Zhang [77] discovered that the relationship between interspecific associations and the degree of niche overlap changed in the central Da Hinggan Mountains as the research scale increased.Wang [62] also observed that as the research scale increased, the niche overlap in the Bayinbuluke alpine meadow did not entirely reflect interspecific associations.Previous studies have shown that individual trees can only influence other plants within a certain range.If the distance between individual trees and other plant species exceeds that range, the influence of trees on other plants will be greatly reduced [78,79].The size of the Wuyi Mountain site may have exceeded the range of influence of some plants.Although some species pairs existed in multiple sites simultaneously, they did not have significant interactions with each other.Therefore, their niche overlap values and interspecific association trends were not the same.

Conclusions
Based on the plant communities in the Wuyi Mountains, this study explores interspecies associations and niches among dominant species in the tree and shrub layers across three elevation ranges.The results indicate that at different elevation ranges, the number of significantly positively correlated dominant species pairs is higher in the tree layer than in negatively correlated pairs.In contrast, the opposite trend is observed in the shrub layer.With increasing elevation, the proportion of significantly negatively correlated species pairs in the shrub layer gradually increases, suggesting that elevation influences interspecies associations.Additionally, the niche overlap between the tree and shrub layers decreases with elevation, with significantly lower overlap at higher elevation ranges compared to mid-to low-elevation ranges.These findings indicate significant differences in resource requirements between species at higher and mid-to low-elevation ranges.Species that demonstrate significant positive associations typically exhibit higher ecological niche overlap.However, some species with high niche overlap display weaker interspecies associations and vice versa.This may be scale-dependent and suggests limitations in the extent of plant influence on other species.
The observed variability and occasional opposition in niche width responses across vertical layers (tree and shrub) underscore the complexity of ecological interactions within subtropical mountainous regions.These contrasting dynamics reflect differential resource utilization strategies and competitive pressures within these distinct vegetation strata.These insights suggest that the mechanisms governing species coexistence and community assembly in subtropical forests are multifaceted and influenced by elevation-driven environmental ranges.The findings not only contribute to a deeper understanding of Wuyi Mountain ecological processes but also offer valuable implications for broader applications in ecological conservation and management of subtropical mountainous regions worldwide.

Species
No.

Figure 1 .
Figure 1.Distribution map of the 34 research sites on Wuyi Mountain National Nature Reserve.The number next to the sample site represents the sample site number, the vertical interval of the contour line is 200m, and the number on it is the elevation of the contour line in meters.

Figure 1 .
Figure 1.Distribution map of the 34 research sites on Wuyi Mountain National Nature Reserve.The number next to the sample site represents the sample site number, the vertical interval of the contour line is 200 m, and the number on it is the elevation of the contour line in meters.
a + b + c + d Note: a represents the number of sample plots in which the paired species co-occur, d represents the number of sample plots in which neither of the paired species occurs, and b and c represent the number of sample plots in which species A and B occur, respectively.The sum of a, b, c, and d equals the total number of sample plots.

Figure 2 .
Figure 2. Half-matrix diagram using association coefficient (AC) values of interspecific associations in the tree and shrub layers: sample plots (a) H1, (b) H2, and (c) H3 of the tree layer and (d) H1,

Figure 2 .
Figure 2. Half-matrix diagram using association coefficient (AC) values of interspecific associations in the tree and shrub layers: sample plots (a) H1, (b) H2, and (c) H3 of the tree layer and (d) H1, (e)

Figure 3 .
Figure 3. Levins's niche width of dominant plant species in tree and shrub layers.Note: The scientific names of trees and shrubs are shown in Table1.Each color block in the outer ring represents a species or an elevation range.The direction of the lines illustrates the dominant species presented in different elevations in these samples, and the thickness of the bar is proportional to the importance value of the dominant species of each sample.Among shrub species, those with larger niche widths included the following: in the low-elevation range, omeiensis, Rhododendron ovatum, Begonia acetosella, and Laurocerasus spinulosa; in the mid-elevation range, Vaccinium sprengelii, Begonia acetosella, Rhododendron simiarum, and Rhododendron tsoi; in the high-elevation range, Enkianthus quinqueflorus, Litsea cubeba, Stranvaesia davidiana, and Vaccinium japonicum.

Figure 3 .
Figure 3. Levins's niche width of dominant plant species in tree and shrub layers.Note: The scientific names of trees and shrubs are shown in Table1.Each color block in the outer ring represents a species or an elevation range.The direction of the lines illustrates the dominant species presented in different elevations in these samples, and the thickness of the bar is proportional to the importance value of the dominant species of each sample.

Figure 4 .
Figure 4. Niche overlap analysis of dominant plant species in tree and shrub layers: sample plots (a) H1, (b) H2, and (c) H3 of the tree layer and (d) H1, (e) H2, and (f) H3 of the shrub layer.The size and color of the circles represent the value of the niche overlap.The legend on the right shows the colors of the niche overlap index from 0-1.The scientific names of trees and shrubs are shown in Table1.

Figure 4 .
Figure 4. Niche overlap analysis of dominant plant species in tree and shrub layers: sample plots (a) H1, (b) H2, and (c) H3 of the tree layer and (d) H1, (e) H2, and (f) H3 of the shrub layer.The size and color of the circles represent the value of the niche overlap.The legend on the right shows the colors of the niche overlap index from 0-1.The scientific names of trees and shrubs are shown in Table1.

Table 1 .
Scientific names and importance values of dominant tree and shrub layer species.

Table 2 .
Sample plot numbers, elevations, and forest community types in each of the three elevation ranges.
presents the measurements of the interspecific associations between any two species A and B. In that table, a represents the number of sample plots in which the paired species co-occur, d represents the number of sample plots in which neither of the paired species occurs, and b and c represent the number of sample plots in which species B and A occur, respectively.The sum of a, b, c, and d equals the total number of sample plots.When ad > bc, then the interspecific association is a positive connection, and if ad < bc, it is a negative connection.If χ 2 > 3.841 (0.01 < p < 0.05), this means that a strong interspecific association exists, and if χ 2 > 6.635 (p < 0.01), this indicates that an extremely significant interspecific association exists

Table 4 .
Overall associations among dominant species.