The influence of environmental factors on species composition and distribution in a community forest in Northern Thailand

Understanding the environmental factors that influence tree species composition is essential for successful management of biodiversity and sustainable use of community forest resources. This study aims to assess tree species composition and distribution in the deciduous Ban Mae Chiang Rai Lum Community Forest in Northern Thailand and to analyze the influence of environmental factors on tree biodiversity in the forest. We conducted a stratified systematic sampling of the forest’s total area of 3925 ha, and twenty-five 0.16 ha survey plots were established in three different stands of the deciduous forests to estimate and characterize the difference in biological diversity among the stands. Canonical correspondence analysis (CCA) was used to investigate the environment factors affecting such differences in biodiversity of the stands. The results showed a high diversity of trees in the forest as 197 species, 144 genera, and 62 plant families were recorded. The CCA ordination identified the environmental factors—the most important of which were elevation, distance to streams, soil moisture, organic matter, and distance to communities—that significantly influenced the diversity and distribution of tree species (p < 0.05) in the community forest. Our findings indicate that the implementation of drought reduction measures such as building check dams, fire protection, and monitoring community forest-product usage would be recommended to further biodiversity conservation and the sustainable use of community forest resources.


Introduction
In Southeast Asia, Community Forest Management (CFM) has been recognized as an effective approach to restore degraded forests, augment the income of the poor, and help meet the needs of local residents for food and medicine (Blair and Olpadwala 1988). CFM has been applied in pursuit of sustainability and forest biodiversity conservation in Thailand (Srisutham and Kaewjampa 2010) as the Royal Forest Department (RFD) has been promoting CFM since 1987. To date, CFM has been implemented in approximately 7% of the country's total forest area, or more than 1.2 million ha, within roughly 17,400 villages nationwide (Royal Forest Department 2014, 2020. Thailand is one of the most bio-diverse countries in Southeast Asia. Forests cover approximately 32% of Thailand's total land area (Royal Forest Department 2019) and are inhabited by roughly 8% of the world's plant species (Office and Natural Resources and Environmental Policy and Planning 2009). However, a decline of forest biodiversity has become a critical issue. There are 1442 threatened plant species in the country, 715 of which have been classified as vulnerable, 207 as endangered and 18 as critically endangered (Department of National Parks, Wildlife and Plant Conservation 2017). Participation by local people in the management of forest resources has been identified as a promising way to address the issue of biodiversity loss in the country.
Although community forest management can facilitate biodiversity conservation, there is insufficient information informing how CFM affects biodiversity and how environmental factors influence that biodiversity in community forests. There are numerous factors that may impact the distribution of trees in community forests.
Climatic factors can threaten the extinction of vulnerable species (e.g. atmospheric pressure, temperature, precipitation, solar radiation, etc.) (Körner 2007;Diem et al. 2018). Climatic factors such as temperature and precipitation are key drivers to control species distribution directly when they exceed the ecophysiological tolerances of species and affect photosynthetic activity and biological processes directly (Rowe 2009). Topographic features such as slope, aspect, and elevation can impact local climate as well as soil conditions that in turn have varied effects on vegetation structure (Zhang et al. 2006;Zhang and Zhang 2007;Zhang et al. 2013Zhang et al. , 2016. The relative distance from a water source can also affect the composition and distribution of woody vegetation because of the resulting varying amount of water available for growth (Marod et al. 1999;Scalley et al. 2009;Tavili et al. 2009;Sarvade et al. 2016;Asanok et al. 2017). Physical and chemical soil properties can inform vegetation patterns on a local scale (Oliveira-Filho et al. 1998;Hejcmanovā-Nežerková and Hejcman 2006;Han et al. 2011;Zhao et al. 2015). For example, higher soil sand ratio lessens water-holding capacity that can lead to water stress on trees (Aranguren et al. 1982;Toledo et al. 2012;Zhang et al. 2013), and acidity levels affect the distribution of species and is linked to slope and elevation in lowland tropical forests (Nguyen et al. 2015;Vahdati et al. 2017). Soil moisture also significantly changes the growth patterns of trees in drought areas (Fu et al. 2004;Yoshifuji et al. 2006;Asanok and Marod 2016;Tilk et al. 2017).
In addition, the organic matter in the soil is relevant to an analysis of environmental factors and plant communities in the forest (Sarker et al. 2014;An et al. 2015;Zhang et al. 2016). Soil nutrients such as nitrogen, potassium, phosphorus, calcium, and magnesium are correlated with the richness and distribution of plant species in tropical forests (Paulo et al. 2007;Zhang et al. 2013;Tilk et al. 2017). Further, the human impact of resource utilization, road construction, and residential development significantly affects plant diversity (Chen et al. 2014;Måren and Sharma 2018;Eghdami et al. 2019).
Overall, the environment's impact on the biodiversity in community forests is complex as it involves a multitude of factors and considerations that provide limited answers, presents additional clues and raise new questions that are crucial to a more thorough and useful understanding of the relationship. This study sought to increase awareness of CFM's role in the furtherance of biodiversity conservation and to seek additional knowledge to support these efforts. To do so, we identified ecological characteristics, ascertained species diversity indices, and investigated pertinent environmental factors impacting three unique forest stands in the community forest. We focused on the following questions: (1) what is the diversity and composition of tree species in the community forest?
(2) what are key environmental factors influencing species distribution in the different stands? and (3) what resulting information could prove to be useful in CFM to protect species diversity and its benefits to livelihoods?

Study area
The Pa Mae Phrik National Forest Reserve is located in Thailand's northern province of Lampang. Situated in the southern region of this Forest Reserve is the Ban Mae Chiang Rai Lum Community Forest (17°22′48″ N-17°27′47″ N and 99°00′47″ E-99°05′48″ E), where in this study was conducted ( Fig. 1).
This community forest has a total area of 3925 ha, is at an elevational range of 140-660 m and is composed of mixed deciduous and dry dipterocarp forest subtypes.
Deciduous forests occur in areas with 5-6 month dry seasons, with a mean annual rainfall of 1000-1500 mm in dry dipterocarp forests and 1000-1800 mm in mixed deciduous forest (Bunyavejchewin 1983;Santisuk 1988). Mixed deciduous forests typically grow on moderate fertile loam soil while dry dipterocarp forests have more sandy lateritic soil (Bunyavejchewin 1983(Bunyavejchewin , 1985Santisuk 1988). Forest fires are a regular feature of dry season in deciduous forests (Sukwong and Dhamanittakul 1977;Bunyavejchewin 1983). A previous study identified Xylia xylocarpa, Schleichera oleosa, Sindora siamensis, Shorea obtusa, Terminalia mucronata as the dominant species in the study area (Royal Forest Department 2017). This community forest has been locally managed in collaboration with the RFD under a community forest project since 2008 and plays a vital role as a source of biodiversity in providing non-timber forest products (NTFPs) while contributing to ecosystem services beneficial to livelihoods. However, unchecked tapping into the forest resources is widespread and creating a potential existential threat to the forest.
Northern Thailand in general, and Lampang in specific, have unique climatic features (Diem et al. 2018). The climate around the study area features a wet season from April to October and a dry season from November to March. Historically, this region also experiences drought conditions. Temperatures ranged from a minimum of 31.7 °C (January and November) to a maximum of 37.1 °C (March). The mean annual temperature, relative humidity, and rainfall were 33.6 °C, 76.1%, 1129.4 mm, respectively (Thai Meteorological Department 2018).

Field survey
A field survey was conducted in the forest from July to October 2018. Results of a 2016 study conducted by the RFD in the same area were used to calculate the sampling intensity with a confidence interval of 95%. An estimate of the sampling was obtained by using the standard deviation from previous surveys (Avery and Burkhart 1983;Asrat and Tesfaye 2013). The formula was expressed as: where n is the sampling intensity; Z is the z-value for the confidence interval of 95%; σ is the density of tree species, and E is the percent standard deviation of the precision required.
A total of 25 sampling plots of 40 m × 40 m (0.16 ha) were surveyed through a systematic sampling method (Asia Network for Sustainable Agriculture and Bioresources 2010). In each plot, trees with a diameter at breast height (DBH) ≥ 4.5 cm were identified and measured in every 10 m × 10 m sub-quadrat. Within 10 m sub-quadrats, samplings with DBH < 4.5 cm and height ˃ 1.30 m were recorded in 4 m × 4 m sub-quadrats, whereas seedlings were documented in 1 m × 1 m sub-quadrats within each 4 m subquadrat (Fig. 2).

Data collection on environmental factors
Climate data were generated by Hobo U23-001 and UA-002-64 which operated in six stations located in plots 3, 10, 13, 15, 17, and 24. The monitoring equipment was spread out and placed in these representative plots to provide the optimal coverage to accumulate climate data from throughout the study area. Historically, relevant climatic data did not (1) n = (Z ∕E) 2 vary significantly during the dry season. In addition, the dry season often featured forest fires. As such, data on temperature, relative humidity and light intensity were collected for 92 days during the wet season from July 11 to October 11, 2018.
Data regarding 15 environmental factors were collected from the 25 sampling plots. Those factors were broadly classified into topographic, edaphic, and anthropogenic. The topographic factors were elevation (EL), slope (SL), and aspect (AS), and the data were extracted out of 1:50,000 topographic maps produced by 20 m contour lines using the ArcGIS. The edaphic factors were investigated from samples collected at a depth of 15-20 cm in three random locations within each 40 m × 40 m plot during the rainy season from June to July 2018. Acidity (pH), organic matter (OM), soil texture [sand (SA), silt (SI), clay (CL)], available phosphorus (P), exchange potassium (K), exchange calcium (Ca), and soil moisture (SM) in the samples were analyzed. Soil pH was measured by 1:1 soil/water mol −1 KCL (National Soil Survey Center 1996). OM was estimated by wet digestion and titration using Walkley-Black method (Nelson and Sommers 1996). Soil texture was determined by Pipette method (Gee and Bauder 1986). P was measured by Bray II (Bray and Kurtz 1945), while K and Ca were extracted with 1 M NH4OAc at pH 7.0 (Thomas 1982), respectively. The soil moisture analysis was performed by oven drying at 105 °C (National Soil Survey Center 1996). ArcGIS's Euclidean distance tool was used to calculate the distance to streams (DS), distance to roads (DR), and distance to communities (DC) from the sampling plot.

Data analysis
The data analysis consisted of examining the biological diversity of the community forest and exploring The forest stands in the community forest were identified by cluster analysis employing importance value index (IVI) matrices in each sampling plot. A cluster analysis requires pruning of the dendrogram at a level representing a compromise between the group and the number of groups. Optimum pruning for the dendrogram was selected by applying the Relative Sorensen Distance and Ward's Linkage Method (McCune and Grace 2002).
To analyze the ecological characteristics of the three stands, we identified mature trees (DBH ≥ 4.5 cm) and calculated the stand density and basal area of individual tree species in each stand type.
Tree species diversity in each individual stand was analyzed, and a comparison of stand-to stand diversity was examined. A rarefaction approach was implemented in order to compare the diversity from the three different data-source stands (Colwell et al. 2012). Individual based interpolation and extrapolation of Hill numbers (q = 1, 2, 3) (Hill 1973) were carried out for each stand and the species diversity among the three forest stands was estimated as follows: species richness (q = 0), Shannon diversity (the exponential of Shannon entropy, q = 1), and Simpson diversity (the inverse of Simpson concentration, q = 2).
The ecological importance of the tree species in each stand was quantified through the IVI using the equation: where R.D is the relative density of the tree species; R.F is their relative frequency, and R.D o is their relative dominance. They were calculated as R.D = number of individuals of the species × 100/total number of quadrate studies, R.F = number of quadrates in which species occurred × 100/ total number of quadrate studies, and R.D o = total basal area of species × 100/total basal area of all the species (Curtis and McIntosh 1951).
To investigate the environmental factors that affect tree species in the community forest, a dataset of 15 factors was tested for normality using the Shapiro-Wilk test for an alpha level of 0.05. The abnormal data was transformed to a normal distribution to improve the linearity. In this study, AS, SI, SM, DR, and DC were normally distributed (p > 0.05). Ten factors had values less than the chosen alpha level; log10(y) was used to transform SL, pH, OM, CL, K, and Ca, and (y)^3, 1/(y), sqrt(y), and − 1/sqrt(y) were used to transform SA, EL, DS, and P, respectively.
The differences in stand characteristics, such as stand density, basal area, and the environmental factors among the three stand types were subjected to a one-way analysis of variance (ANOVA) with a post hoc Tukey's honestly significant difference (HSD) test (Zhang et al. 2016). (2) Species compositions in the three forest stands were compared using the Jaccard similarity index (Magurran 1988). The correlations between environmental factors were determined through Pearson correlation analysis. In case of a high correlation in which a matrix of two independent variables has correlation coefficients (r) > 0.7 across all environmental factors, only one of the factors was selected for the ordination analysis (Sarker et al. 2014).
The relationship between tree species distribution and environmental factors was finally assessed using the canonical correspondence analysis (CCA) (McCune and Grace 2002). The matrices between the IVI of mature trees (DBH ≥ 4.5 cm) in each stand and the environmental factors were analyzed in the CCA ordination. The correlation significance between the matrices was determined using the Monte Carlo test with 999 permutations. Only the environmental factors significantly related (Monte Carlo permutation test, p < 0.05) to the mature tree species were retained. All the statistical calculations were performed using PC. ORD version 5.10 (McCune and Mefford 2006) and the R program version 3.5.1 for Windows software (R Development Core Team 2019).

Climate data
Relative humidity and light intensity data did not vary significantly plot-to-plot, but temperature data was significantly different (p < 0.001). The average temperature during the 92-day collection period varied between 25.09 °C (plot 17) and 27.37 °C (plot 10). The average relative humidity ranged from 79.19% (plot 15) to 85.10% (plot 17), and the average light intensity ranged from 3011.67 l × (plot 15) to 4698.48 l × (plot 17). The overall average temperature, relative humidity, and light intensity throughout the study area were 26.85 °C, 82.23%, and 3549.09 l ×, respectively (Fig. 4).

Species composition and diversity
The field survey recorded 18,567 trees comprising 197 species (129 tree, 99 sapling, and 141 seedling species), 144 genera, and 62 plant families. The Jaccard index showed a 21.29% similarity (33 species) between the MLS and SSS, a 36.30% similarity (53 species) between the SSS and SOS, and a 17.78% similarity (24 species) between the MLS and SOS. The forest had an average density of 966 trees h −1 and an average basal area of 16.74 m 2 h −1 (Table S1). The one-way ANOVA test showed that the average density were significantly different stand-to-stand (p < 0.01), though average basal area did not vary significantly.
The five species with the highest IVI in the MLS were M. leucantha, L. duperreana, M. brandisiana, Antidesma sootepense, and Pterocarpus macrocarpus. The most dominant tree species in terms of its IVI in the SSS was S. siamensis, followed by S. obtusa, X. xylocarpa, T. mucronata, and Dipterocarpus tuberculatus. S. obtusa was the most common species in the SOS, and it included S. siamensis, X. xylocarpa, Canarium subulatum, and Ellipanthus tomentosus (Table S2).
The rarefaction-produced tree species diversity values in different sample sizes of three forest stands are shown in Table S3. A comparison of extrapolated species diversity shows that species richness was estimated to be the highest in the MLS, followed by the SSS and the SOS, respectively. The Shannon diversity and Simpson diversity also estimates are highest in the MLS, the SSS, and the SOS, in descending order (Fig. 5).

Relationship between environmental factors and species diversity
The mean values and standard deviation of the 15 environmental factors analyzed in the Ban Mae Chiang Rai Lum Community Forest are shown in Table S4. The degree of soil acidity (pH), organic matter (OM), silt (SI), clay (CL), exchange potassium (K), soil moisture (SM), distance to roads (DR), and distance to communities (DC) were greater in the MLS than in the SSS and SOS. The values of elevation (EL), slope (SL), available phosphorus (P), and exchange calcium (Ca) were highest in the SSS, and in decreasing order in the MLS and SOS. The value of aspect (AS), distance to streams (DS), and sand ratio (SA) were greater in the SOS than in the SSS and MLS. Five factors differed significantly among the three stand types; SL, OM, SA, CL, and K differed from one stand type to another (p < 0.05). There was no significant difference in EL, AS, DS, pH, SI, P, Ca, SM, DR, and DC among the stands. The post hoc Tukey's HSD test identified that the CL value was different between the MLS and the SOS, between the SOS and the SSS, and between the SSS and the MLS. OM and SA differed substantially between the MLS and the SOS, but they did not differ in the SSS. Similarly, SL differed between the SOS and the SSS, but it did not differ in the MLS. K varied when comparing the MLS with the SOS and the SSS with the SOS.
A matrix of two independent environmental factors showed that there were high correlation coefficients (r > 0.7) between the factors (Table S5). OM was positively correlated with K (r = 0.932). SA was positively associated with EL (r = 0.704), but negatively associated with SI (r = − 0.890), CL (r = − 0.929), DR (r = − 0.746), and K (r = − 0.724). CL positively correlated to DR (r = 0.710), but negatively correlated to EL (r = − 0.707). EL negatively related to SL (r = − 0.740). DC was positively related with DR (r = 0.853). Considering those correlations between the environmental factors, the ten variables, EL, AS, DS, pH, OM, P, Ca, SM, SI, and DC, were selected for further analysis using CCA. In the preliminary CCA, two poorly correlated factors, Ca and SI, were eliminated; the remaining eight factors were used in the final CCA.  (Table 1). The eigenvalues obtained for the first and second axes were 0.531 and 0.363, respectively. The Pearson correlation showed that the environmental factors had a high correlation to the species for axis 1 (0.924) and axis 2 (0.917). Therefore, these axes were considered as good indicators of relationships between the environmental factors and the species distribution in the community forest. The Monte Carlo permutation test confirmed that the first axis explained the significant variation in the community forest (p < 0.05). The significance of the second and third axes was not reported for the reason that a simple randomization test for these axes could bias the p values (McCune and Mefford 2006).
The CCA results demonstrate the relationships between environmental factors and species distribution (Figs. 6 and  7). In Grouping 1, the EL was an important factor in species distribution in plots 2, 9, and 11 of the SSS and in plots 1, 7, 8, and 12 of the SOS. The abundant species in this grouping were Atalantia monophylla (ATMO), Antidesma ghaesembilla (

Discussion
The Ban Mae Chiang Rai Lum Community Forest was divided into three forest stands: the mixed deciduous MLS and the dry dipterocarp SSS and SOS (Fig. 3). A rarefaction prediction of species diversity (Fig. 5 and Table S3) and the wide array of species, genera, and families (Table S1) recorded via the field survey indicate remarkable diversity, similar to other forests of Northern Thailand (Popradit et al. 2015;Khamyong et al. 2018). There were no noteworthy differences in the average basal areas, but average density varied significantly among the three stand types. There was a degree of species similarity among the stands whereas the forest structures and species compositions were distinct (Table S2), In addition, the importance of the environment factors varied stand-to-stand (Table S4) as did their degree of correlation (Table S5). Our results represent that the heterogeneity of species distribution patterns among the different forest stand types could be a response to a variation in the environmental factors; distinct environmental conditions in forest stands can be effective indicators of species distribution in the Ban Mae Chiang Rai Community Forest in Northern Thailand.
The results of the CCA demonstrated that species composition and distribution of the mature trees were related to topographic, edaphic, and anthropogenic factors ( Table 1). The environmental factors showed both positive and negative impacts on species diversity. Specifically, EL, DS, SM, OM, and DC were the most evidently linked to species distribution (Table 2). Moreover, environmental factors had different effects on species distribution between stand types (Figs. 6 and 7). This indicates that species composition and distribution can be limited by a combination of environmental factors across distinct forests.
Among the five significant factors (i.e. EL, DS, SM, OM, and DC), EL appeared to be the most influential on species distribution. The IVI of mature trees was positively related to elevation in the SSS and SOS; elevational gradients promote not only species diversification but also soil differentiation (Zhang and Zhang 2007;Zhang et al. 2013Zhang et al. , 2016. Though, climatic factors can vary with elevation (Zhang et al. 2006;Körner 2007), the elevational features of dry dipterocarp in Thailand occur mainly in elevations up to 900 m (Bunyavejchewin et al. 2011). At this elevation, there is no drastic change in lowland vegetation similar to that found in most tropical zones at elevations above 1000 m, including in the North of Thailand (Küchler and Sawyer 1967). Thus, Grouping 1 Fig. 6 The CCA ordination diagram representing the relationship between the distribution of tree species and the environmental factors. Full names of tree species by abbreviation are listed in Table S6 the various soil properties unique to elevation could be key to the composition and distribution of tree species in the dipterocarp forest stands this study area (Bunyavejchewin et al. 2011).
Other tree species associated with dominant species were also found in the dry dipterocarp stand. This suggests that higher elevation promotes species diversity in dry dipterocarp forests. Teejuntuk et al. (2003) reported that the diversity of tree species in Doi Inthanon National Park, Northern Thailand increased up to an elevation of about 1800 m. Also, it has been found that tree species in dry dipterocarp forests can move into higher land to grow with other tree species in an ecotone zone (Kutintara 1975;Marod et al. 2015) which may be related to a more favorable environment for growth at slightly higher elevation (Ivanov et al. 2008).
Elevational gradients often limit the distribution of tree species by decreasing soil moisture and nutrient levels (Bridge and Johnson 2000). Several studies have reported that an increase in species diversity may have an elevational limit above which diversity would become negatively correlated (Liu et al. 2018). For instance, the resources necessary for plant growth are limited at extremely high elevations because of strong winds and shallow soils (Zhang et al. 2013).
In this study area, elevation ranged between 140 and 660 m. Higher elevation and SI, CL, and K were negatively correlated whereas SA was positively correlated (Table S5). Poor sandy soil is an important driver for tree species in dry dipterocarp forests (Bunyavejchewin 1983(Bunyavejchewin , 1985Santisuk 1988). Therefore, some tree species in the SSS and the SOS may be able to grow at higher elevation because of more favorable soil conditions without the limiting climatic conditions found at extremely high elevation.
In this study, it appears that DS and SM influenced species diversity and distribution. The positive correlation between species diversity in the SSS and SOS with the DS revealed that diversity is related to water availability. This indicates that dry conditions may negatively impact species diversity and distribution. Several species in the dry dipterocarp forests are strongly associated with proximity to water sources. This suggests that tree species adapt to drought areas that have low soil moisture, high aridity, and fire disturbance. These are important factors that determine the species occurrence in dry dipterocarp forests in Thailand (Rundel and Boonprakob 1995;Marod et al. 2002;Bunyavejchewin et al. 2011).
Similarly, the negative correlation of species diversity in the MLS and SSS with SM substantiated the occurrence of these forests in dry regions (Marod et al. 1999). Tree species such as P. macrocarpus (PTMA), P. cerasoides (POCE), Z. dongnaiensis (ZODO), H. orixense (HYOR), and L. coromandelica (LACO) were closely associated with SM, which suggests that they were restricted by moisture levels and sensitive to drought conditions. Generally, these species had a higher water demand than other tree species in deciduous forests that could adapt and survive in arid areas (Marod et al. 1999). Previous findings indicate that SM decreased as the distance from a reservoir increased (Sarvade et al. 2016), suggesting that the species in MLS and SSS were also related to distance from a water source.
Deciduous forests often occur in arid areas with sandy soil; 53% of the soils in mixed deciduous forests and 70% of those in deciduous dipterocarp forests consist of sand (Myo et al. 2016). Sandy soils are more susceptible to leaching which can lead to drought (Aranguren et al. 1982;Zhang et al. 2013). Soils with higher sand content have lower waterholding capacity, thereby causing plants to compete more for soil moisture (Toledo et al. 2012). Contrarily, trees growing in clay soil have more efficient root system and are less susceptible to drought stress because of the greater capacity of the soil for water retention. Trees with deeper roots in sandy soil have a greater chance of surviving as they adapt to reduce drought stress (Wessel 1971). Hence, a deep root system enables plants to adapt to and survive in extremely dry conditions, such as being distant from a water supply or being in an area with limited soil moisture. In this study, our results imply that trees survive and grow well despite the dry conditions in the SSS and SOS; the sufficient depth of root system of trees in the both stands has enabled them to adapt to drought stress. However, the existence some species in the MLS and SSS may be limited by the level of moisture in soil for their growth.
The negative correlation of the dominant species in the MLS and SSS with OM and DC, explains that both environmental factors constrain tree growth and diversity in their corresponding forest stands. The growth, species distribution, and composition of trees in tropical areas are influenced by soil nutrients (John et al. 2007;Santiago et al. 2012). Many studies have shown that OM is a key environmental factor in plant communities (Zhang and Zhang 2007;Slik et al. 2009;Sarker et al. 2014), and that it contributes to the availability of nutrients and water by improving soil structure and physical conditions, increasing the water holding capacity of the soil, and providing a habitat for plant roots and soil organisms (Carter 2002;Meng et al. 2014). Thus, soils with more OM are more fertile and more favorable for the optimal growth of trees (Vahdati et al. 2017).
In the study forest, however, there were tree species in the MLS and SSS that were inversely associated with OM, indicating that this environmental factor was restricting their growth. As explained by Grime (1977), if nutrients are scarcely available in the soil, then plants will have to adapt to low-resource and nutrient-poor conditions, making natural selection predominant in the process.
Plant productivity is positively linked with OM, as reported by Bauer and Black (1994), thus, more OM should be made available for plants to increase their productivity. Our findings also revealed that OM had a positive correlation with pH and P (Table S5), indicating that trees are not only affected by OM content but that pH and P are major soil properties involved in plant growth. In a previous study conducted in a tropical forest in Ben En National Park, Vietnam, it was reported that an increase in OM could augment P and pH, thereby enhancing species composition and distribution of plants in the area (Hoang et al. 2011).
The dominant species were closely correlated with DC, supporting a previous report that human disturbance is a threat to plant diversity loss (Millennium Ecosystem Assessment 2005). Not only were dominant species found distributed in the stands, but shrub species such as W. tinctorial (WETI), M. rotundifolia (MIRO), B. retusa and Z. oenoplia (ZIOE) were also distributed near the communities. These species are often used for firewood and construction. As such, there is a high risk of them being harvested thereby affecting their distribution. Furthermore, they are pioneer species that play a vital role in ecological succession after disturbances (Department of National Parks, Wildlife and Plant Conservation 2007; Asanok et al. 2020). Utilization by communities may have an impact on the natural restoration of forest ecosystems in the community forest.
Many studies have shown that the impact on tree species in forests could be worse if they were situated closer to a community (Thapa and Chapman 2010;Hoang et al. 2011;Chen et al. 2014;Asanok et al. 2017;Martínez-Camilo et al. 2018). These studies reported that trees would be more abundant in less disturbed plots located far from human settlements. Furthermore, Teejuntuk et al. (2003) discussed that at an elevation less than 1800 m in Northern Thailand, the occurrence of ecologically similar species increased at higher elevation since forests at higher elevation are relatively far and inaccessible to man. Lowland forests, conversely, have become significantly more fragmented and degraded because of deforestation through illegal logging, burning, and grazing. Human disturbance is more rampant in areas with easy access, poor monitoring and weak enforcement of regulations (Måren and Sharma 2018). The apathy of locals to forest management also aggravates the problem and contributes to the decline of forest biodiversity.

Conclusions
Our study suggests that elevation, distance to streams, soil moisture, organic matter, and distance to communities are the most important factors influencing species composition and distribution in the deciduous forests of the Ban Mae Chiang Rai Lum Community Forest in Northern Thailand.
Even through tree species in deciduous forests are able to shift to higher elevation to find more favorable conditions and adapt themselves to grow in drought areas, the limitation of organic matter, water availability, and soil moisture can still impact tree species diversity. This could prove to be important and useful information for community forest management. We recommend that management practices for drought reduction such as building check dams and fire protection be implemented to protect biodiversity, especially in the MLS and SSS.
In addition, biodiversity conservation would be better served by monitoring the utilization of forest resources by nearby communities. Species propagation, silviculture, and transplantation should be initiated and/or expanded, particularly in neighboring villages or private farms, to enhance the productivity and viability of the community forest. However, future studies are needed to investigate the threat that community utilization of trees poses to tree species' very survival as this was not an area of investigation in this study.