Next Article in Journal
Variation in Water Uptake Dynamics of Dominant Wood Plants of Pinus taiwanensis Hayata Communities Based on Stable Isotopes
Next Article in Special Issue
Habitat Conditions and Tree Species Shape Liana Distribution in a Subtropical Forest
Previous Article in Journal
Effect of Leptographium terebrantis on Foliage, New Root Dynamics, and Stemwood Growth in a Loblolly Pine (Pinus taeda L.) Plantation
Previous Article in Special Issue
Whole-Plant Seedling Functional Traits Suggest Lianas Also Support “Fast-Slow” Plant Economics Spectrum
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 13
Issue 8
10.3390/f13081334
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Habitat Type and Scale Dependences of Interspecific Associations in a Subtropical Evergreen Broad-Leaved Forest

by
Changchun Jiang
1,†,
Jiaqin Fu
1,†,
Yunquan Wang
1,*,
Pengtao Chai
1,
Yidan Yang
1,
Xiangcheng Mi
2,
Mingjian Yu
3,
Keping Ma
2 and
Jianhua Chen
1
1
College of Chemistry and Life Sciences, Zhejiang Normal University, Jinhua 321004, China
2
State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
3
Ministry of Education Key Laboratory of Conservation Biology for Endangered Wildlife, College of Life Sciences, Zhejiang University, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Forests 2022, 13(8), 1334; https://doi.org/10.3390/f13081334
Submission received: 6 July 2022 / Revised: 18 August 2022 / Accepted: 19 August 2022 / Published: 21 August 2022
(This article belongs to the Special Issue Long-Term Monitoring of Forest Biodiversity and Dynamics in China)
Browse Figures
Review Reports Versions Notes

Abstract

:
“Interspecific associations” refers to the interrelationship among different species in a particular spatial distribution, which plays an important role in species distribution, community assembly, and responses to environmental changes. However, the strength and/or direction of interspecific associations may vary with environmental gradients and scales. To understand the effects of habitat types and research scales on interspecific associations in subtropical forests, we modeled the interspecific associations for more than 15,000 individuals representing 74 co-occurring species from three habitat types and three scales by using the variance ratio and the Spearman rank correlation coefficient. We found that overall interspecific associations at a community level exhibited significant positive associations for most habitat types and scales. Moreover, interspecific associations of pairwise species have strong habitat dependence, and the association strengths decreased with the increase in elevation (change in habitat types). However, the scale dependence of pairwise interspecific associations varies with habitat types. The strength of interspecific associations increased with the increasing scales (10 m × 10 m, 20 m × 20 m, and 40 m × 40 m) at low-valleys and mid-hillside habitats, while the scale-dependent effect was not detected at high-ridges. In conclusion, our study highlights the importance of environmental gradients and research scales on interspecific associations in diverse subtropical forests, and environmental gradients and research scales should be considered in future studies.

1. Introduction

“Interspecific association” refers to the interactions and spatial relationships among different species occupying a habitat, which can provide prominent insights into community assembly as well as determine community structure and function [1,2,3,4]. This concept is widely used to quantify interspecific relationships and to infer ecological processes underlying community assembly [5,6,7,8,9]. Previous studies have shown that the interspecific associations of plant communities are not static and are usually constrained by a variety of influencing factors such as interspecific competition [10], environmental changes [11,12], and research scales [8]. However, little is known about how the strengths and/or directions of interspecific associations vary with changing environmental gradients and research scales.
Environmental factors are important in that they affect the interspecific associations of communities, among which factors such as elevation, convexity, and slope have important influences on the spatial pattern and species richness of communities [13,14,15]. For example, Chesson [16] suggested that the status of dominant species may be more affected by habitat changes than by other factors in multispecies-competing communities. Zhang et al. [17] found significant correlations between vegetation patterns and habitat factors such as elevation and slope in the Taihang Mountains in China. Moreover, Long and Tang suggested that the effects of topographic factors (such as elevation and slope) on spatial structure vary with the changing topographic factors in an evergreen broad-leaved forest [18]. These findings suggest that the strengths and directions of interspecific associations change along environmental gradients.
Research scale is also a critical factor affecting the strengths of the interspecific associations of communities [8,19]. Previous studies have shown that research scales should be matched to the selected plot area, plant community type, and vegetation uniformity when detecting interspecific associations [20]. For instance, Deng found that the overall interspecific associations at community level changed from negative to positive correlation accordingly when the research scales changed from the small scales of 100 m2 and 200 m2 to the large scale of 400 m2 in a mixed coniferous and broad-leaved forest [21]. Deng et al. found that with the gradual increase in the scale, the percentage of significant positive interspecific associations of a secondary forest had a clear trend of first increasing and then slowly decreasing [22]. Thus, research scales should be considered when quantifying the interspecific associations of communities.
The subtropical evergreen broad-leaved forest is a typical zonal vegetation type in China. The richness of the vegetation diversity of the subtropical evergreen broad-leaved forest is second only to that of the tropical rainforest and much higher than that of other evergreen broad-leaved forests [23,24,25]. The interspecific interactions and biodiversity maintenance mechanisms in subtropical forests have been core questions for ecologists in China [24,26,27]. This provides an ideal site for studying how interspecific associations vary with changes in habitat type and research scale.
To evaluate the effects of habitat type and research scale on interspecific associations in subtropical forests, we examined how the strength and/or direction of interspecific associations for more than 15,000 individuals representing 74 co-occurring species change with habitat type and research scale in a subtropical forest plot in China. Specifically, we focused on the following questions:
(1) How do the overall interspecific associations at community level vary with habitat type and research scale in the subtropical forest?
(2) What is the habitat dependence of the pairwise interspecific associations in the subtropical forest?
(3) What is the scale dependence of the pairwise interspecific associations in the subtropical forest?

2. Materials and Methods

2.1. Study Site and Tree Census

The study was carried out in the predominantly old-growth subtropical forest of the 5-ha forest dynamic plot (FDP; 29.25° N, 118.12° E) in the Qianjiangyuan National Park, China. This region is located in the mid-subtropical humid monsoon climate zone in China, which has rich rainfall in summer, warm and humid weather in winter, and abundant microclimate habitats [28]. The mean annual temperature in this region is 15.3 °C, and the mean annual precipitation is 1964 mm [12]. The main soil types in this region are red soil, yellow-red soil, red-yellow soil, and, in some areas, swamp soil [28]. Qianjiangyuan National Park is rich in forest vegetation types, with evergreen broad-leaved forests, evergreen and deciduous broad-leaved mixed forests, coniferous broad-leaved mixed forests, coniferous forests, and subalpine wetland [29]. Among these, the mid-subtropical evergreen broad-leaved forest is the main vegetation in the park, a typical mid-subtropical zonal vegetation that is widely distributed in low-elevation areas below 800 m above sea level.
The 5-ha forest dynamic plot was established in a large ditch between the two main hills in accordance with CTFS-ForestGEO protocol [30]. The plot is 200 m long from east to west and 250 m wide from north to south (Figure 1). All freestanding woody stems with DBH (diameter at breast height, 1.3 m) ≥ 1 cm in the plot were tagged, mapped, and identified to species during the summer of 2002 [31]. Survivorship was assessed and new recruits were tagged, measured, mapped, and identified for each 5-year interval [12]. The 5-ha plot was dominated by evergreen tree species, accounting for 52.45 % of the total species, and the most dominant species in this plot was Castanopsis eyrei Champ. ex Benth. Hutch., 1905 [31].
In the present study, we used the 2012 tree census data to assess the habitat type and scale dependence of interspecific association in the 5-ha subtropical forest. In total, there were 17,398 woody plants belonging to 158 species in the 5-ha forest plot. In order to study the habitat dependence and scale dependence of the interspecific associations of the pairwise species, we analyzed all the co-occurring species in all the habitat types and scale cells. Here we modeled the interspecific associations of more than 15,000 individuals (>86% of total individuals) belonging to 74 wood species in the 5-ha subtropical forest (Table 1).

2.2. Environmental Variables and Habitat Categories

In this study, we defined environmental variables in terms of topography. Four topographic factors were identified: elevation, convexity, slope, and aspect for each 20 × 20 m cell, following Harms et al. [32] and Valencia et al. [33]. To test how interspecific associations vary with habitat categories, these four topographic variables were selected by Multivariate regression tree (MRT) procedure by using mvpart package [34], and a tree with three habitat categories was selected as the best one after 1000 cross-validation trials following Chen et al. [35]. Based on these four variables, the whole 5-ha forest plot was classified into three habitat categories with 20 m × 20 m cells (Figure 2): low-valleys (H1, 50 plots), mid-hillsides (H2, 45 plots), and high-ridges (H3, 35 plots). Although four habitat factors—including elevation, convexity, slope, and slope direction—were used as independent variables in the analysis, the multivariate regression tree was successfully constructed only by the factor of elevation. This means that elevation played a dominant role in the habitat classification of the present study and that elevation increases with the increasing number of habitat types.
In the present study, to test how interspecific associations vary with different scales, the 5-ha forest plot was grouped into cells 10 m × 10 m, 20 m × 20 m, and 40 m × 40 m in size. There were only 8 cells 80 m × 80 m, too few for statistical analysis. Based on these three scales, the entire 5-ha forest plot was divided into different scale cells (See Figures S1 and S2 for more details). These cells were used to quantify the scale dependence of interspecific associations in the 5-ha forest plot.

2.3. Statistical Analysis

Community-level interspecific association: The variance ratio test (VR) was used to determine the overall interspecific association at community level. The VR value indicates whether there is a significant relationship among multiple species in the selected area, and the significance of the VR value was tested by the statistic W [1]. The VR value was calculated as:
V R = 1 N ∑ i = 1 N ( T i − t ) 2 ∑ j = 1 S ( 1 − P j )
where VR is the variance ratio, N is the total number of plots in the community, T i is the total number of target species in plot i. t is the observed mean number of species in all the plots. S represents the total number of species in the community, P j represents the frequency of species j, and P j = n j N , n j represents the total number of plots occupied by species j [1].
Under the null hypothesis of independence, the expected value of VR is 1; that is, when VR = 1, it means that there is no connection between the species. If VR > 1, it means that there is a positive connection between the species. If VR < 1, it means that there is a negative connection between the species. The statistic W was used to verify the significant degree of VR deviation from 1, W = V R × N . If W > χ 0.05 2 ( N ) or W < χ 0.95 2 ( N ) , it means that the overall connection between the species is significant. Conversely, the overall connection between the species is nonsignificant.
Pairwise interspecific association: The Spearman rank correlation coefficient was used to test the interspecific association between pairwise species and quantitatively analyze the covariant linear relationship between species pairs following Bishara & Hittner [36]:
r ( i , j ) = 1 − 6 ∑ p = 1 N ( x i p − x ‾ i ) 2 ( x j p − x ‾ j ) 2 N 3 − N
where r ( i , j ) represents the spearman rank correlation coefficient species i and species j in plot p; N represents the total number of plots in the community; x i p and x j p represent the rank of species i and species j in plot p, respectively; and x ‾ i and x ‾ j represent the mean abundance of species i and species j in all the plots. The value range of r ( i , j ) is from −1 to 1, with positive values indicating positive correlation and negative values indicating negative correlation.
All the analyses were conducted in the R 4.1.2 statistical platform (http://www.r-project.org/, accessed on 1 November 2021). Interspecific association was calculated using the function in the spaa package [37]. The significance of the Spearman rank correlation coefficients was calculated using the rcorr() function in the Hmisc package [38]. The visualization of matrix saliency was implemented using the chart.Correlation() function in the Performance Analytics package [39]. Boxplots were used to show the change in interspecific association across habitat categories or scales, and the ggplot2 package [40] and ggsignif package [41] were used to draw boxplots and calculate the significance between different groups.

3. Results

3.1. Overall Interspecific Associations at Community Level

In order to study the habitat type and scale dependence of overall interspecific associations at community level, we evaluated the overall interspecific association of the community in three habitats and three scales by using the variance ratio method (Table 2, 9 groups in total). The overall interspecific association at community level for most habitat types and scales exhibited positive associations (VR > 1) while showing a negative association in habitat H3 at a scale of 20 m × 20 m (VR = 0.846 < 1). The W statistics and χ 2 tests were employed to test the significance of the deviation of the VR values from 1. The overall interspecific associations at community level showed significant associations for most of the habitat types at different scales (Table 2, W > χ 0.05 2 ( N ) or W < χ 0.95 2 ( N ) ). However, overall community-interspecific associations were not significant for habitat H3 at either 20 m × 20 m or 40 m × 40 m scales.

3.2. Habitat Dependence of Interspecific Associations for Pairwise Species

To understand the habitat type–dependence of interspecific associations for pairwise species in this forest community, we calculated the pairwise species–interspecific associations of 74 co-occurring species from three habitat types at three scales by using the Spearman rank correlation coefficient (Figure 3). The interspecific association of the pairwise species exhibited a positive association at three scales for two habitat types (Figure 3, H1 and H2), while the interspecific associations were negative (tended to be non-associated) in habitat H3 at a scale of 10 m × 10 m and 40 m × 40 m (Figure 3). The strength of the pairwise interspecific associations decreased with the increase in elevation (from H1 to H3). Moreover, the strength of the pairwise interspecific associations varied significantly among most of the habitat types at three scales, indicating that habitat type significantly affects pairwise interspecific associations in this forest. However, there were no significant associations between habitat H2 and H3 at a 40 m × 40 m scale (Figure 3c).

3.3. Scale Dependence of Interspecific Associations for Pairwise Species

To explore the scale dependence of interspecific associations for pairwise species in the forest community, we evaluated the pairwise species–interspecific associations of 74 co-occurrence species at 3 scale levels for each of 3 habitats types (Figure 4). The interspecific association of the pairwise species for most of the scales exhibited positive associations in all the habitat types (Figure 4). However, the interspecific associations were negative (tended to be non-associated) at a scale of 10 m × 10 m and 40 m × 40 m in habitat H3. The strength of the pairwise interspecific associations increased with the increase in research scales in habitat H1 and H2 (Figure 4). In addition, the strength of the pairwise interspecific associations changed significantly among all the scales in habitat H1 and H2 (Figure 4a,b). For habitat H3, however, a significant difference was detected only between the scales of 10 m × 10 m and 40 m × 40 m (Figure 4c). These results indicated that there were significant scale dependences in the subtropical forest in Qianjiangyuan National Park, yet the scale dependence of the pairwise interspecific associations varied with habitat types.

4. Discussions

Interspecific associations are known to influence species distribution, community assembly, and responses to environmental changes [4,8,42,43]. The importance of interspecific associations on the forest community has been studied at length, but how interspecific associations depend on environmental gradients and research scales and the nature of their interactions remain unclear. In this study, we revealed how interspecific associations depend on environmental gradients and research scales in a subtropical evergreen broad-leaved forest by using the variance ratio and the Spearman rank correlation coefficient.

4.1. Overall Interspecific Associations in the Gutianshan Subtropical Forest

The overall interspecific associations are known to reflect the interspecific competition among species and the stability of species composition and community structure in the community [9,44,45,46]. We found that the overall interspecific associations were significantly positive in the old-growth forest (Table 2). The results indicated that the Gutianshan forest was in a relatively stable succession stage in which community structure and species composition tend to be complete and dynamically balanced [3,45]. Our findings highlight the importance of environmental filtering: Species in this forest tend to share similar environmental space and resources, resulting in overall positive associations and ultimately achieving stable coexistence [3,47]. Moreover, we also found negative overall interspecific associations for some specific habitats such as H3 at scale of 20 m × 20 m, but the effect was not statistically significant. This pattern could be due to the fact that species have different survival strategies in different habitats [11,48], and the final manifestation is that the interspecific associations change with environmental gradients such as elevation [9].

4.2. Habitat Dependences of Pairwise Interspecific Associations in the Subtropical Forest

Our analysis revealed significant habitat dependences of pairwise interspecific associations in the Gutianshan subtropical forest (Figure 3). The pairwise interspecific associations were positive in two of three habitat types (H1 and H2). The results of the pairwise positive associations we found here are consistent with the results of the overall interspecific associations we discussed above, indicating that environmental filtering plays an important role in pairwise species interactions in the Gutianshan forest [12,35,49]. Conversely, the pairwise interspecific associations tended to be negative- or non-associated in habitat H3. These results may due to the fact that, in addition to the effects of environmental filtering, interspecific competition for limited spaces and resources also plays an important role in determining the outcome of interspecific associations [6,11,50]. Thus, a negative relationship due to species competition can partly mask the positive relationship resulting from environmental filtering. The final detected strength of pairwise interspecific associations depends on the relative importance of the environmental filtering and the species competition.
More importantly, we found that the strength of pairwise interspecific associations varied significantly among most of the habitat types at all three scales, indicating strong habitat dependences of interspecific associations in the community. This is consistent with studies reporting that pairwise interspecific associations significantly changed with changes in altitude or environmental condition [9,43,48]. As habitat type, topography, and/or soil conditions changed greatly, species density also varied significantly across habitat types [51], finally resulting in a significant difference in species interactions across habitats types [35,43]. Interestingly, we found a decrease in pairwise interspecific associations as elevation increased (from H1 to H3), and this pattern is consistent across all three research scales. The results indicated that environmental filtering may play more important roles at a lower elevation habitat relative to the higher elevation habitat in the Gutianshan subtropical forest.

4.3. Scale Dependence of Interspecific Associations for Pairwise Species

Pairwise interspecific associations may not only vary across environmental gradients but also change along research scales, and scale is crucial when detecting and predicting interspecific associations [8,47]. It remains unclear how the strengths of pairwise interspecific associations vary from a small to a large scale [19]. The pairwise interspecific associations were positive for most scales and tend to be negative at other scales (Figure 4), consistent with the results of overall and pairwise interspecific associations we have discussed above.
The strength of pairwise interspecific associations varied significantly among all the scales in habitat H1 and H2 (Figure 4a,b), which indicated that research scale is an important driver of interspecific associations [8,47]. However, for habitat H3, the strength was significantly changed only between the scales of 10 m × 10 m and 40 m × 40 m (Figure 4c). Our results indicated that there were significant scale dependences in the Gutianshan subtropical forest, yet the scale dependences of pairwise interspecific associations vary with habitat types. Additionally, we found that the strength of positive pairwise interspecific associations increased with the increase in research scales in habitat H1 and H2. This pattern could be due to the fact that interspecific competition for limited space and resources mainly focuses on small scales, that neighborhood interactions are typically negligible beyond 20 m, and that species interactions are influenced by environmental factors beyond this scale [52].

5. Conclusions

The habitat and scale dependence of the interspecific associations in the subtropical forest reported here provide insight into how overall and pairwise interspecific associations vary with environmental gradients and research scales [20]. The interspecific associations were significantly affected by habitat type and research scale in the subtropical forest in Qianjiangyuan National Park. Community-level overall interspecific associations were significantly positive for most of the habitat types and scales. Pairwise interspecific associations have significant differences between habitat types, indicating a strong habitat type dependence of interspecific associations for pairwise species in this forest. Moreover, the interspecific association strengths decreased with the increase in elevation (from H1 to H3). For a scale-dependent effect, however, the scale dependence of pairwise interspecific associations varies with habitat types. The scale-dependent effects were detected at low-valleys (H1) and mid-hillside (H2) habitats, but not at high-ridges (H3). Our results suggest that environmental gradients and research scales significantly influence interspecific associations in diverse subtropical forests, and environmental gradients and research scales should be considered in future studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/f13081334/s1, Figure S1: Habitat categories of the 5-ha forest plot at 10 m × 10 m scale. The white color represents low-valleys (H1), the blue color represents mid-hillsides (H2), while the yellow color represents high-ridges (H3); Figure S2: Habitat categories of the 5-ha forest plot at 40 m × 40 m scale. The white color represents low-valleys (H1), the blue color represents mid-hillsides (H2), while the yellow color represents high-ridges (H3).

Author Contributions

C.J., Y.W. and J.C. conceived and designed the study; C.J. and J.F. performed statistical analyses and wrote the first draft with substantial input from Y.W. and J.C.; all the other authors provided data and contributed to the development of the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the Zhejiang Provincial Natural Science Foundation of China (LQ22C030001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available from the authors on request.

Acknowledgments

We acknowledge the Center of Ecology and Resources, Qianjiangyuan National Park for support of this study, as well as the hard work of the hundreds of people who were involved in the collection of the vast quantity of data in the Gutianshan plot over the past few years.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schluter, D. A variance test for detecting species associations, with some example applications. Ecology 1984, 65, 998–1005. [Google Scholar] [CrossRef]
  2. Wiegand, T.; Gunatilleke, S.; Gunatilleke, N. Species associations in a heterogeneous Sri Lankan dipterocarp forest. Am. Nat. 2007, 170, E77–E95. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, Z.W.; Zhu, Y.; Wang, J.J.; Ma, W.; Meng, J.H. Species association of the dominant tree species in an old-growth forest and implications for enrichment planting for the restoration of natural degraded forest in subtropical China. Forests 2019, 10, 957. [Google Scholar] [CrossRef]
  4. Keil, P.; Wiegand, T.; Tóth, A.B.; McGlinn, D.J.; Chase, J.M. Measurement and analysis of interspecific spatial associations as a facet of biodiversity. Ecol. Monogr. 2021, 91, e01452. [Google Scholar] [CrossRef]
  5. Maihaiti, M.; Zhang, W.J. A mini review on theories and measures of interspecific associations. Selforganizology 2014, 1, 206–210. [Google Scholar]
  6. Dai, J.Y.; Liu, H.Y.; Xu, C.Y.; Qi, Y.; Zhu, X.R.; Zhou, M.; Liu, B.B.; Wu, Y.H. Divergent hydraulic strategies explain the interspecific associations of co-occurring trees in forest-steppe ecotone. Forests 2020, 11, 942. [Google Scholar] [CrossRef]
  7. Fridley, J.D.; Lynn, J.S.; Grime, J.P.; Askew, A.P. Longer growing seasons shift grassland vegetation towards more-productive species. Nat. Clim. Chang. 2016, 6, 865–868. [Google Scholar] [CrossRef]
  8. Ding, Z.Q.; Ma, K.M. Identifying changing interspecific associations along gradients at multiple scales using wavelet correlation networks. Ecology 2021, 102, e03360. [Google Scholar] [CrossRef]
  9. Jin, S.S.; Zhang, Y.Y.; Zhou, M.L.; Dong, X.M.; Chang, C.H.; Wang, T.; Yan, D.F. Interspecific Association and Community Stability of Tree Species in Natural Secondary Forests at Different Altitude Gradients in the Southern Taihang Mountains. Forests 2022, 13, 373. [Google Scholar] [CrossRef]
  10. Kim, T.N.; Underwood, N.; Inouye, B.D. Insect herbivores change the outcome of plant competition through both inter-and intraspecific processes. Ecology 2013, 94, 1753–1763. [Google Scholar] [CrossRef]
  11. Yuan, Z.L.; Wei, B.L.; Chen, Y.; Jia, H.R.; Wei, Q.N.; Ye, Y.Z. How do similarities in spatial distributions and interspecific associations affect the coexistence of Quercus species in the Baotianman National Nature Reserve, Henan, China. Ecol. Evol. 2018, 8, 2580–2593. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, Y.Q.; Cadotte, M.W.; Chen, J.H.; Mi, X.C.; Ren, H.B.; Liu, X.J.; Yu, M.J.; Zhang, J.T.; Ma, K.P. Neighborhood interactions on seedling survival were greatly altered following an extreme winter storm. For. Ecol. Manag. 2020, 461, 117940. [Google Scholar] [CrossRef]
  13. Burke, I.C.; Lauenroth, W.K.; Riggle, R.; Brannen, P.; Madigan, B.; Beard, S. Spatial Variability of Soil Properties in the Shortgrass Steppe: The Relative Importance of Topography, Grazing, Microsite, and Plant Species in Controlling Spatial Patterns. Ecosystems 1999, 2, 422–438. [Google Scholar] [CrossRef]
  14. Zhang, C.S.; Xie, G.D.; Bao, W.K.; Chen, L.; Pei, S.; Fan, N. Effects of topographic factors on the plant species richness and distribution pattern of alpine meadow in source region of the Lancang River, Southwest China. Chin. J. Ecol. 2012, 31, 2767–2774. [Google Scholar]
  15. Yang, W.J.; Wang, Y.H.; Webb, A.A.; Li, Z.Y.; Tian, X.; Han, Z.T.; Wang, S.L.; Yu, P.T. Influence of climatic and geographic factors on the spatial distribution of Qinghai spruce forests in the dryland Qilian Mountains of Northwest China. Sci. Total Environ. 2018, 46, 1007–1017. [Google Scholar] [CrossRef]
  16. Chesson, P. Multispecies competition in variable environments. Theor. Popul. Biol. 1994, 45, 227–276. [Google Scholar] [CrossRef]
  17. Zhang, J.; Xi, Y.; Li, J. The relationships between environment and plant communities in the middle part of Taihang Mountain Range, North China. Community Ecol. 2006, 7, 155–163. [Google Scholar] [CrossRef]
  18. Long, J.S.; Tang, M.P. Relationship between spatial structure and terrain factors of evergreen broad-leaved forest in Mount Tianmu. J. Zhejiang AF Univ. 2021, 38, 47–57. [Google Scholar]
  19. Ovaskainen, O.; Abrego, N.; Halme, P.; Dunson, D. Using latent variable models to identify large networks of species-to-species associations at different spatial scales. Methods Ecol. Evol. 2016, 7, 549–555. [Google Scholar] [CrossRef]
  20. Xu, M.H.; Liu, M.; Zhai, D.T.; Liu, T. A review of contents and methods used to analyze various aspects of plant interspecific associations. Acta Ecol. Sin. 2016, 36, 8224–8233. [Google Scholar]
  21. Deng, H.J. Species diversity and interspecific association of Pinus tabuliformis mixed forest in Huanglong Mountain. M.D. Thesis, Northwest A&F University, Xianyang, China, 2015. [Google Scholar]
  22. Deng, L.P.; Bai, X.J.; Li, L.L.; Niu, S.S.; Han, M.N.; Qin, S.J.; Zhou, Y.B. Interspecific association and correlation among dominant woody plants of secondary forest in montane region of eastern Liaoning province, China. Chin. J. Ecol. 2015, 34, 1473–1479. [Google Scholar]
  23. Xu, X.N.; Hirata, E.; Tokashiki, Y.; Shinohara, T. Structure and species diversity of subtropical evergreen broad-leaved forest in northern Okinawa Island, Japan. J. For. Res. 2001, 6, 203–210. [Google Scholar] [CrossRef]
  24. Song, Y.C.; Chen, X.Y.; Wang, X.H. Studies on evergreen broad-leaved forests of China: A retrospect and prospect. J. East Chin. Norm. Univ. Nat. Sci. Edit. 2005, 1, 1–8. [Google Scholar]
  25. Zhu, H.; Zhou, S.S.; Yan, L.C.; Shi, J.P.; Shen, Y.X. Studies on the evergreen broad-leaved forests of Yunnan, southwestern China. Bot. Rev. 2019, 85, 131–148. [Google Scholar] [CrossRef]
  26. Li, Y.D.; Xu, H.; Chen, D.X.; Luo, T.S.; Mo, J.H.; Luo, W.; Chen, H.Q.; Jiang, Z.L. Division of ecological species groups and functional groups based on interspecific association—a case study of the tree layer in the tropical lowland rainforest of Jianfenling in Hainan Island, China. Frint. For. Chin. 2008, 3, 407–415. [Google Scholar] [CrossRef]
  27. Chen, L.; Swenson, N.G.; Ji, N.N.; Mi, X.C.; Ren, H.B.; Guo, L.D.; Ma, K.P. Differential soil fungus accumulation and density dependence of trees in a subtropical forest. Science 2019, 366, 124–128. [Google Scholar] [CrossRef]
  28. Hu, Z.H.; Yu, M.J.; Ding, B.Y.; Fang, T.; Qian, H.Y.; Chen, Q.C. Types of evergreen broad-leaved forests and their species diversity in Gutian Mountain National Nature Reserve. Chin. J. Appl. Environ. Biol. 2003, 9, 341–345. [Google Scholar]
  29. Yu, M.J.; Hu, Z.H.; Yu, J.P.; Ding, B.Y.; Fang, T. Forest vegetation types in Gutianshan Natural Reserve in Zhejiang. J. Zhejiang. Univ. Agric. Life Sci. 2001, 27, 375–380. [Google Scholar]
  30. Condit, R. Tropical Forest Census Plots: Methods and Results from Barro Colorado Island, Panama and a Comparison with Other Plots; Springer: New York, NY, USA, 1998. [Google Scholar]
  31. Wang, Y.H.; Mi, X.C.; Chen, S.W.; Li, M.H.; Yu, M.J. Regeneration dynamics of major tree species during 2002–2007 in a subtropical evergreen broad-leaved forest in Gutianshan National Nature Reserve in East China. Biodivers. Sci. 2011, 19, 178. [Google Scholar]
  32. Harms, K.E.; Condit, R.; Hubbell, S.P.; Foster, R.B. Habitat associations of trees and shrubs in a 50-ha neotropical forest plot. J. Ecol. 2001, 89, 947–959. [Google Scholar] [CrossRef]
  33. Valencia, R.; Foster, R.B.; Villa, G.; Condit, R.; Svenning, J.C.; Hernández, C.; Romoleroux, K.; Losos, E.; Magård, E.; Balslev, H. Tree species distributions and local habitat variation in the Amazon: Large forest plot in eastern ecuador. J. Ecol. 2004, 92, 214–229. [Google Scholar] [CrossRef]
  34. Therneau, T.; Atkinson, B. Mvpart: Multivariate Partitioning. R Package Version 1.6-2. 2014. Available online: https://cran.r-project.org/web/packages/mvpart/index.html (accessed on 24 June 2014).
  35. Chen, L.; Mi, X.C.; Comita, L.S.; Zhang, L.W.; Ren, H.B.; Ma, K.P. Community-level consequences of density dependence and habitat association in a subtropical broad-leaved forest. Ecol. Lett. 2010, 13, 695–704. [Google Scholar] [CrossRef] [PubMed]
  36. Bishara, A.J.; Hittner, J.B. Testing the significance of a correlation with nonnormal data: Comparison of Pearson, Spearman, transformation, and resampling approaches. Psychol. Methods 2012, 17, 399–417. [Google Scholar] [CrossRef]
  37. Zhang, J.L. Spaa: Species Association Analysis. R Package Version 0.2.2. 2016. Available online: https://cran.r-project.org/web/packages/spaa/index.html (accessed on 9 June 2016).
  38. Harrell, F.E., Jr. Hmisc: Harrell Miscellaneous. R Package Version 4.6-0. 2021. Available online: https://cran.r-project.org/web/packages/Hmisc/index.html (accessed on 13 October 2021).
  39. Peterson, B.G.; Carl, P.; Boudt, K.; Bennett, B.; Ulrich, J.; Zivo, E.; Cornilly, D.; Hung, E.; Lestel, M.; Balkissoon, K.; et al. PerformanceAnalytics: Econometric Tools for Performance and Risk Analysis. R Package Version 2.2.4. 2020. Available online: https://cran.r-project.org/web/packages/PerformanceAnalytics/index.html (accessed on 6 February 2020).
  40. Wickham, H.; Chang, W.; Henry, L.; Pedersen, T.L.; Takahashi, K.; Wilke, C.; Woo, K.; Yutani, H.; Dunnington, D. ggplot2: Create Elegant Data Visualisations Using the Grammar of Graphics. R Package Version 3.3.4. 2021. Available online: https://cran.r-project.org/web/packages/ggplot2/index.html (accessed on 16 June 2021).
  41. Ahlmann-Eltze, C.; Patil, I. Ggsignif: Significance Brackets for ‘ggplot2’. R Package Version 0.6.3. 2021. Available online: https://cran.r-project.org/web/packages/ggsignif/index.html (accessed on 9 September 2021).
  42. Gotzenberger, L.; Bello, F.D.; Bråthen, K.A.; Davison, J.; Dubuis, A.; Guisan, A.; Lepš, J.; Lindborg, R.; Moora, M.; Pärtel, M.; et al. Ecological assembly rules in plant communities-approaches, patterns and prospects. Biol. Rev. Camb. Philos. 2012, 87, 111–127. [Google Scholar] [CrossRef] [PubMed]
  43. Borah, B.; Quader, S.; Srinivasan, U. Responses of interspecific associations in mixed-species bird flocks to selective logging. J. Appl. Ecol. 2018, 55, 1637–1646. [Google Scholar] [CrossRef]
  44. Su, S.J.; Liu, J.F.; He, Z.S.; Zheng, S.Q.; Hong, W.; Xu, D.W. Ecological species groups and interspecific association of dominant tree species in Daiyun Mountain National Nature Reserve. J. Mt. Sci. 2015, 12, 637–646. [Google Scholar] [CrossRef]
  45. Gu, L.; Gong, Z.W.; Li, W.Z. Niches and interspecific associations of dominant populations in three changed stages of natural secondary forests on Loess Plateau, P.R. China. Sci. Rep. 2017, 7, 6604. [Google Scholar] [CrossRef]
  46. Haak, C.R.; Hui, F.K.C.; Cowles, G.W.; Danylchuk, A.J. Positive interspecific associations consistent with social information use shape juvenile fish assemblages. Ecology 2020, 101, e02920. [Google Scholar] [CrossRef]
  47. Hao, Z.Q.; Zhang, J.; Song, B.; Ye, J.; Li, B.H. Vertical structure and spatial associations of dominant tree species in an old-growth temperate forest. Forest Ecol. Manag. 2007, 252, 1–11. [Google Scholar] [CrossRef]
  48. Ma, Y.M.; Li, Q.H.; Pan, S.P.; Liu, C.; Han, M.S.; Brancelj, A. Niche and interspecific associations of Pseudoanabaena limnetica − Exploring the influencing factors of its succession stage. Ecol. Indic. 2022, 138, 108806. [Google Scholar] [CrossRef]
  49. Jin, Y.; Russo, S.E.; Yu, M.J. Effects of light and topography on regeneration and coexistence of evergreen and deciduous tree species in a Chinese subtropical forest. J. Ecol. 2018, 106, 1634–1645. [Google Scholar] [CrossRef]
  50. Zhou, Q.; Shi, H.; Shu, X.; Xie, F.L.; Zhang, K.R.; Zhang, Q.F.; Dang, H.S. Spatial distribution and interspecific associations in a deciduous broad-leaved forest in north-central China. J. Veg. Sci. 2019, 30, 1153–1163. [Google Scholar] [CrossRef]
  51. Guo, Y.; Wang, Y.Q.; Chen, L.; Mi, X.C.; Ren, H.B.; Chen, S.W.; Chen, J.H. Comparing tree seedling composition and distribution patterns under different sampling intensities in the 24 ha Gutianshan forest dynamics plot. Biodivers. Sci. 2016, 24, 1093. [Google Scholar] [CrossRef]
  52. Hubbell, S.P.; Ahumada, J.A.; Condit, R.; Foster, R.B. Local neighborhood effects on long-term survival of individual trees in a neotropical forest. Ecol. Res. 2001, 16, 859–875. [Google Scholar] [CrossRef]
Forests 13 01334 g001 550
Figure 1. (a) Location of Qianjiangyuan National Park, and (b) lateral view and (c) top view of the 5-ha forest plot in Qianjiangyuan National Park.
Figure 1. (a) Location of Qianjiangyuan National Park, and (b) lateral view and (c) top view of the 5-ha forest plot in Qianjiangyuan National Park.
Forests 13 01334 g001
Forests 13 01334 g002 550
Figure 2. Habitat categories of the 5-ha forest plot at 20 m × 20 m scale. The white represents low-valleys (H1), the blue represents mid-hillsides (H2), and the yellow represents high-ridges (H3).
Figure 2. Habitat categories of the 5-ha forest plot at 20 m × 20 m scale. The white represents low-valleys (H1), the blue represents mid-hillsides (H2), and the yellow represents high-ridges (H3).
Forests 13 01334 g002
Forests 13 01334 g003 550
Figure 3. Boxplots of interspecific associations of three habitats at different scales: 10 m × 10 m scale (a), 20 m × 20 m scale (b), and 40 m × 40 m scale (c). Asterisks indicate a significant difference in different habitat types (*** p < 0.001; NS, nonsignificant).
Figure 3. Boxplots of interspecific associations of three habitats at different scales: 10 m × 10 m scale (a), 20 m × 20 m scale (b), and 40 m × 40 m scale (c). Asterisks indicate a significant difference in different habitat types (*** p < 0.001; NS, nonsignificant).
Forests 13 01334 g003
Forests 13 01334 g004 550
Figure 4. Boxplots of interspecific associations of plots at three scales from different habitat types: low-valleys (a), mid-hillsides, (b) and high-ridges (c). Asterisks indicate a significant difference in different habitat types (* p < 0.05; *** p < 0.001; NS, nonsignificant).
Figure 4. Boxplots of interspecific associations of plots at three scales from different habitat types: low-valleys (a), mid-hillsides, (b) and high-ridges (c). Asterisks indicate a significant difference in different habitat types (* p < 0.05; *** p < 0.001; NS, nonsignificant).
Forests 13 01334 g004
Table
Table 1. Basic information about the habitat categories of the 5-ha forest plots at different scales.
Table 1. Basic information about the habitat categories of the 5-ha forest plots at different scales.
Research ScalesHabitat CategoriesHabitat
Types		No. of
Plots		Total Area
(ha)		No. of
Species		No. of
Co-Occurring Species		Individuals
/Plot
H1Low-valleys20021467436.58
10 m × 10 mH2Mid-hillsides1801.81347429.37
H3High-ridges1001827424.73
H1Low-valleys50214674146.30
20 m × 20 mH2Mid-hillsides451.813474117.49
H3High-ridges251827498.92
H1Low-valleys111.7614474580.36
40 m × 40 mH2Mid-hillsides101.613074454.50
H3High-ridges50.88074383.20
Notes: N: the total number of plots in the selected area community; H1: low-valleys; H2: mid-hillsides; H3: high-ridges.
Table
Table 2. The overall interspecific associations at community level of three scales and habitats.
Table 2. The overall interspecific associations at community level of three scales and habitats.
Research ScalesHabitat CategoriesVariance Ratio (VR)W Statistic χ 0.95 2 ( N ) ,   χ 0.05 2 ( N ) Test
Results
H12.278410.006149.969, 212.304significant
positive correlation
10 m × 10 mH25.017993.402166.444, 231.829significant
positive correlation
H31.478147.80477.929, 124.324significant
positive correlation
H12.506112.75730.612, 61.656significant
positive correlation
20 m × 20 mH26.378318.90334.764, 67.505significant
positive correlation
H30.84621.13914.611, 37.652Nonsignificant negative correlation
H12.64426.4393.940, 18.307significant
positive correlation
40 m × 40 mH24.79752.7724.575, 19.675significant
positive correlation
H31.1485.7411.145, 11.070Nonsignificant positive correlation
Note: N: the total number of plots in the selected area community; H1: low-valleys; H2: mid-hillsides; H3: high-ridges.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jiang, C.; Fu, J.; Wang, Y.; Chai, P.; Yang, Y.; Mi, X.; Yu, M.; Ma, K.; Chen, J. The Habitat Type and Scale Dependences of Interspecific Associations in a Subtropical Evergreen Broad-Leaved Forest. Forests 2022, 13, 1334. https://doi.org/10.3390/f13081334

AMA Style

Jiang C, Fu J, Wang Y, Chai P, Yang Y, Mi X, Yu M, Ma K, Chen J. The Habitat Type and Scale Dependences of Interspecific Associations in a Subtropical Evergreen Broad-Leaved Forest. Forests. 2022; 13(8):1334. https://doi.org/10.3390/f13081334

Chicago/Turabian Style

Jiang, Changchun, Jiaqin Fu, Yunquan Wang, Pengtao Chai, Yidan Yang, Xiangcheng Mi, Mingjian Yu, Keping Ma, and Jianhua Chen. 2022. "The Habitat Type and Scale Dependences of Interspecific Associations in a Subtropical Evergreen Broad-Leaved Forest" Forests 13, no. 8: 1334. https://doi.org/10.3390/f13081334

APA Style

Jiang, C., Fu, J., Wang, Y., Chai, P., Yang, Y., Mi, X., Yu, M., Ma, K., & Chen, J. (2022). The Habitat Type and Scale Dependences of Interspecific Associations in a Subtropical Evergreen Broad-Leaved Forest. Forests, 13(8), 1334. https://doi.org/10.3390/f13081334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop