Forests in Nepal are extremely important for supporting the livelihood of millions of people who collect forest products for their subsistence use and partly for income generation. Such inherent dependence is expected to cause disturbance in the forest ecosystem. We investigated changes in the structural assemblages caused by the interaction between anthropogenic disturbances and forest management activities in the mixed forests of Sal (Shorea robusta Gaertn.) of Terai, Central Nepal. We evaluated three buffer zone community forests (BZCFs), namely, Radha Krishna, Musharni Mai, and Janajagaran of Parsa Wildlife Reserve (PWR); the forest inside PWR was taken as a control. A transect of 2 km length was laid in each forest, and six plots, each of 1 ha size, were established at a successive interval of 300 m along the edge to the interior of the forests to count and record the diameter at breast height (DBH) of the studied plants. We observed that the species diversity increased linearly (p < 0.05) towards the forest interior in the BZCFs. Species other than S. robusta had significantly higher (p < 0.05) dominance and Importance Value Indices in the interior sites. We did not observe such trends in the control forest. Multivariate analysis showed that the sites of BZCFs had higher structural dissimilarity, but the control forest sites were closer to each other in composition. The forest sites near the settlements had undergone biotic homogenization (S. robusta mixed forest changed to S. robusta forest) due to the interaction between anthropogenic disturbances and forest management activities. On the basis of vegetation density, the edges of BZCFs appeared to be protected, but on the basis of diversity failed to do so. Future management strategies should be directed towards enhancing the diversity, heterogeneity, and forest quality, especially near the forest edges.

1. Introduction

Anthropogenic activities influence the ecological processes [13] and have extensive impacts on the forest ecosystems [47]. The fragmentation of natural forests creates a pattern in the landscapes and is regarded as one of the most important causes for the recent decline of forest species [8], and the human activities increase the risk of further degradation in the forests sharing borders with anthropogenic uses [912]. However, with an increasing concern towards forest degradation, management interventions are also underway with a focus on the restoration of degraded areas [13]. Human pressures from the nearby settlements and forest management activities interact to alter the forest structural and functional attributes [14]. Our concern was to determine the extent of such changes in the mixed forest of Sal (Shorea robusta Gaertn.; family: Dipterocarpaceae) of Terai, Central Nepal.

Terai forests of Nepal, dominated by S. robusta, are among the most disturbed ecosystems [15, 16]. Commonly occurring forests fires every year during dry months also kill young regenerations and burn the deciduous litter layer. In addition, these ecosystems are important for millions of people who depend on forest products for their subsistence (Table 1). Repeated extraction of forest resources to meet livelihood needs of adjoining human settlements can cause significant impacts on the forest regeneration, structure, and diversity, and these impacts often have been associated with the distance to the settlements [1720]. In our studied buffer zone community forests (BZCFs), besides harvesting of the forest resources, the forests are also under regular management activities carried by Buffer Zone Management Committees, and the interaction between the disturbances and management interventions might have altered the ecological parameters differently. To our understanding, there are ample studies showing the impacts of forest disturbance on forest structure [1619], but limited research exists on how management and disturbance interaction influences the composition, structure, richness, and taxonomic diversity of the forests. What ecological changes take place due to such interactions in the forests of Nepal are not known [21].

We attempted to study ecological changes in the Sal forest of Nepal as a result of interaction between inherent disturbances and ongoing management practices (for more than 20 years) in the BZCFs of Nepal. To understand responses of the woody plant species in the forests, our study aimed at elucidating (1) present condition of natural BZCFs in terms of forest community composition, distribution, density, richness, and species diversity as a function of distance from the human settlements, (2) changes in basal area (dominance) of woody species due to management and disturbance interventions from the edge to the forest interior, and (3) effect of anthropogenic disturbances and management interventions on the similarity and dissimilarities of the sites located along the distance from the human settlements. We hypothesized that forest structure differs when management interventions are applied in conjugation with the ongoing disturbances in the forests.

2. Materials and Methods

2.1. Study Area

Cross-sectional vegetation survey was conducted in three BZCFs of Parsa Wildlife Reserve (PWR) (PWR has been currently designated as a national park) at Gadimai Municipality, Central Terai Nepal, namely, Janajagaran BZCF (area: 371.10 ha), Musharni Mai BZCF (area: 231.31 ha), and Radha Krishna BZCF (area: 621.171 ha) (Figure 1). The forest in PWR was considered as a control. These forests are the part of the largest continuous forest landscape in Terai region, also known as “Charkose Jhadi” in Nepal. The BZCFs are situated in the southern part of PWR (Bhimeshwor BZCF is the other forest located in the southern part of PWR) on alluvial plains and composed of flood or river deposits of sand, silt, clay, coarser sediments, and conglomerates. The BZCFs are adjacent to each other with the connected boundaries between them. Since the forests are connected, the transects laid in the forests were kept at least at the distance of 1.5 km. The area harbors several protected faunal species such as Elephants (Elephas maximus), Royal Bengal Tiger (Panthera tigris), Gauri Gai (Bos gaurus), Wild Dogs (Cannes aurous), and many other common wild animals including spotted deer (Axis axis), barking deer (Muntiacus muntjac), hog deer (Axis porcinus), and wild boar (Sus scrofa). Common reptiles found in the area are golden lizard (Calotes versicolor), monitor lizard (Varanus bengalensis), krait (Bungarus careleus), branded krait (Bungarus faciatus), cobra (Naja naja), and common rat snake (Ptyas mucosus).

The climate in the area is tropical. A typical year can be divided into three main seasons: cold, hot, and rainy. April and May have the highest mean maximum temperature, while the coldest months are November, December, and January. Mean monthly minimum and maximum temperature are 18°C and 30°C, respectively. The average monthly rainfall in the area is 156 mm. Maximum rainfall occurs during the monsoon (June-September) season.

2.2. Experimental Design and Sampling

A preliminary survey for the study was carried out during July 2015. Primary data collection and detail inventory of BZCFs and the protected forest inside PWR were carried out from July 2016 to October 2016. Using ArcGIS (10.2.2) and Google Earth Images (Google Earth, an azimuth for each forest was oriented towards the forest interior, away from the corresponding settlement areas. A constant azimuth (355°, 350°, 1°, and 275° for Janajagaran BZCF, Musharni Mai BZCF, Radha Krishna BZCF, and PWR, respectively) was taken and a transect line was laid along the azimuth (Appendix A, Table 1). Similar sampling design was also adopted in other primary studies [7, 16, 17, 22].

Six plots each of 1 ha size were established on the left side of each transect at a successive interval (horizontal distance) of 300 m along the edge to the interior of each forest (150 m: Site I, 450 m: Site II, 750 m: Site III, 1050 m: Site IV, 1350 m: Site V, and 1650 m: Site VI). The farthest distance from the settlement (Site VI) consisted of relatively undisturbed stands (however, there were signs of illegal cutting of some individual stems). For the control forest, a transect was laid starting near the Central Office of PWR along the main entrance (fire line) towards the forest core. We used a handheld GPS (Garmin 60CSx) to determine the position of the plot-centers along the transect line (rope of 100 m was also stretched to ease transect layout). Each plot was then divided into 10 subplots each of 20 × 50 m for sampling [23]. The total number of sampling units (subplots) was sufficient to produce stabilizing species-area curves in the forests (Appendix B, Figure 13). In all the studied forests, surface fire occurs at least once during the dry season of the year which usually destroys smaller seedlings, so we considered the woody species ≥ 1.5 cm in diameter at breast height (DBH) which have survived in moderate fire pressure. Structural attributes were collected from a total of 60 subplots in each studied BZCF. The woody species present were identified, counted, and tallied, and DBH of each individual ≥ 1.5 cm DBH was measured at 1.3 m above the ground with the help of a DBH tape (Kinglon Diameter Tape, Japan). Based on DBH, individuals were classified as saplings (1.5 to 10 cm DBH), poles (10 to 20 cm DBH), and trees (> 20 cm DBH). For a multistemmed tree at breast height, it was treated as a single individual whereby the diameters of all stems were measured separately and averaged. In case a tree had a buttress or an abnormality at 1.3 m height, the diameter was measured just above the abnormality.

Most plants were identified in situ; plants which were not identified to species in the field were identified by the National Herbarium and Plant Laboratories (KATH Herbarium Center) at Godavari, Lalitpur, Nepal (letter reference number 198 2073/2074). The species which could not be identified were considered as “unidentified”.

2.3. Data Analysis

Basal area of live woody species (dominance) was calculated as Σπr2, where r is the radius (DBH/2) for all live woody species with DBH ≥ 1.5 cm. The density of a species was calculated as the number of individuals of species per unit area. Similarly, the estimated maximum density of mature S. robusta trees that 1 ha forest area can support was calculated with the mean value of the area covered by canopy projections of mature trees in natural growth conditions and considering that there will be no overlap of canopy between individuals of the species. The canopy projections (length) over cardinal and ordinal directions were measured based on the observation made vertically upwards and recorded as to whether or not the canopy of the trees obscures the sky [24]. The canopy cover area was calculated based on horizontal projections data and using readWKT and gArea functions in the rgeos package [25] and then was compared with 1 ha area.

The Importance Value Index (IVI) for each woody species was determined to evaluate the contribution of woody species in the forest composition, structure, and biomass. The IVI for a species was calculated by summing the relative frequency, relative density, and relative dominance [2628]:where is the importance value of the jth species, is the frequency of jth species, N is the sum of frequencies (N = ∑), is the number of individuals of the jth species present in the sample population, D is the total number of individuals in the sample population (D = ∑), is the sum of dominance for the jth species, and X is the total dominance across all species X = ∑. Since relative frequency, relative density, and relative dominance are expressed in percentage, the maximum value of IVI is equal to 300.

Species diversity (alpha diversity) was calculated by using Shannon-Wiener index [29] and it is defined aswhere pi is the proportional abundance of a species i and b is the base of the logarithm. The diversity function in the vegan package uses natural logarithms as the default while calculating Shannon-Wiener index, and 1– Simpson’s alpha (which also describes overall alpha diversity) is calculated as 1 – [30]. All analyses for determining diversity (Shannon-Wiener, Simpson, and species accumulation curves) were performed using the diversity function in the vegan package [31].

To observe the functional association between forest community structure and distance from boundary lines of villages, a linear model was fit wherever applicable with the null hypothesis that there is no linear pattern between the structural attributes of forests and the distance from human settlements. Linear model assumptions were tested using the gvlma function in the gvlma package [32]. Residuals versus fitted and Q-Q plots were visualized.

Morisita’s Iδ index of dispersion was used for determining the spatial pattern of individuals of the predominant woody species. Iδ, which is independent of sample size and diversity [33], was computed as follows [34]:where q is the total number of subplots, N is the total number of individuals of the species in q subplots, and x is the number of individuals of a species in a subplot. Species were considered predominant if IVI > 10 observed in any one of the sites of the forests. We used dispindmorisita function in the vegan package to calculate intraspecific aggregation [31].

To examine the dissimilarity and similarity in species composition between the sites in the forests, Bray-Curtis Index (BCI) was used. BCI is a rank-order index and values were used to determine a standardized species score to combine the measurements of IVI in the woody plant species. The values of BCI range from 0 to 1. The value 1 reflects no similarity whereas 0 reflects complete similarity between plots. BCI of dissimilarity was calculated between the sites (from the edge to the interior of the forests) based on the IVI values determined for each species. Vegdist and hclust functions in the vegan package were used for determining dissimilarity and clustering (Ward.D2 group linkage), respectively [31]. Dissimilarity () between sites j and k using BCI can be expressed as follows [35]:where and represent measures of species i in the sample sites j and k and abs indicates absolute value. To compare the dissimilarity in the forest community between sites, a matrix representation of the and the similarity ( = 1 – ) between the paired sites were tabulated and compared. A multivariate agglomerative clustering using Ward.D2 group linkage based on the Bray-Curtis distance of IVI values was performed.

We also determined the associated plant community composition and similarity between sites based on IVI of the species using Nonmetric Multidimensional Scaling (NMDS) ordination to collapse information from multiple species into just a few, so that they can be visualized and interpreted. In community ecology, NMDS is commonly regarded as the most robust unconstrained ordination method [36] and considered as an indirect gradient analysis approach for producing an ordination based on a distance or dissimilarity matrix. We chose Bray-Curtis distance measure and a maximum of 200 iterations for NMDS distance matrix. The analysis was formed with the metaMDS function in the vegan package [31]. To minimize the stress (usually the stress values more than 0.2 should be interpreted with caution and more than 0.3 are highly suspicious), NMDS was projected in K=3 dimensions (resulting in stress value for all our analysis less than 0.05; Shepard stress plots were visualized). All statistical analyses were performed in R [37].

3. Results

3.1. Community Structure and Distribution of Woody Plants

The studied BZCFs and PWR are S. robusta dominated mixed forests. The forests exhibit high total densities of woody plants (saplings, poles, and trees) but also high variance within forest sites (Figure 2). Total density (± SD) of woody species in Janajagaran BZCF ranged from 231±79 stems ha−1 to 658±97 stems ha−1 varying across different distances towards the forest interior from the boundary lines of the villages. Similarly, Musharni Mai BZCF and Radha Krishna BZCF also showed the wide range of variation in stem density (± SD) from 286±45 to 753±159 stems ha−1 and 323±112 stems ha−1 to 643±71.81 stems ha−1, respectively. In PWR, the range of variation of densities (± SD) was from 607±113 to 1021±311 stems ha−1 (Appendix C, Tables 36). The median stem density was more than 400 individuals ha−1 (arbitrary value) in each site of PWR, but oscillations were observed across the sites in BZCFs. The variance of frequency indicates the abundance difference of woody species, within a site; high values indicate that the abundance of species in the plot differs from each other, whereas low values indicate that all the species have a similar pattern of abundance (Appendix C, Tables 36).

Significant pattern (p > 0.05) was not observed between the total density of species and the distance from the boundary lines of the villages in all the BZCFs (Figure 3). When stem densities were separated into the density of the dominant species S. robusta and the density of other woody species, a significant increasing linear trend was obtained between the density of other species (except S. robusta) and the distance from the boundary lines of the settlement areas (Figure 3) (Janajagaran BZCF: p < 0.05, Musharni Mai BZCF: p < 0.05, and Radha Krishna BZCF: < 0.01, but no significant difference was observed for PWR: p > 0.05). Similarly, the assemblage of S. robusta, when categorized as saplings, poles, and trees, reflected that the poles and trees had a high density near the human settlement areas in BZCFs (Figure 4). The BZCFs near the human settlements were receiving intense management and protection of the poles and tree-sized individuals of S. robusta. Densification has occurred in the area with very fewer plant species. The density of trees of S. robusta in the edges was very high above the estimated maximum density (85 stems ha−1) of fully grown mature trees of S. robusta which 1 ha forest area can support. For PWR such difference in the density of S. robusta along the edge to the interior was not observed (Figure 4).

We found a total of 69 woody species belonging to 30 families in all the sample plots of the studied forests. Plant species belonging to the Fabaceae family were the most abundant (12) followed by those belonging to the families Euphorbiaceae (7), Combretaceae (4), and Rubiaceae (4). The number of species varied within the sites from the forest edges to the interior. The numbers of woody species increased significantly as a function of the distance from the village boundaries for Janajagaran (p < 0.05) and Radha Krishna (p < 0.05) BZCFs. In the case of Musharni Mai BZCF, the number of species observed was the least for the site near the edge but a nonsignificant pattern was observed (p>0.05). For the control forest, such a trend was not observed (p > 0.05) (Figure 5). In BZCFs, species like Dalbergia sissoo, Dalbergia latifolia, Wendlandia tinctoria, Bombax ceiba, Toona ciliata, Trewia nudiflora, and Premma barabata were recorded only in the innermost sites.

The cross-sectional area (basal area) per unit sampling area is an indicator of woody species dominance. Basal area correlates with maturity and age of the woody plant species. Our study showed mixed patterns between total dominance and distance from the boundary lines of villages (Figure 6). A decreased linear trend was observed in Janajagaran BZCF (p < 0.05), a linear trend was absent in Musharni Mai BZCF, and an increased trend was observed in Radha Krishna BZCF (p < 0.001). Though no linear trend was present in the Musharni Mai BZCF, a decrease in the basal area moving from the edge to the interior of the forest was observed. As expected, PWR showed the oscillating pattern of increase and decrease of the total basal area along the sites. An interesting result was observed when the total basal area was segregated as contributed by S. robusta and the basal area of other species except for S. robusta. In BZCFs edges, the majority of the basal area was contributed by S. robusta but its value was observed decreasing when moving across the interior of the forests. In contrast, the basal area of other species showed an increasing trend (p < 0.05) in BZCFs (Figures 6 and 7) except for the control forest. Based on % basal area, the forest near the settlement areas was observed to have undergone biotic homogenization (S. robusta Terai Mixed Hardwood Forest to S. robusta Forest) due to disturbance and management interactions.

We calculated the Importance Value Index (IVI) for each species in every 1 ha plot. Though the forest was dominated by S. robusta, other species were identified as predominant species based on high IVI values (IVI > 10). Janajagaran BZCF, Musharni Mai BZCF, Radha Krishna BZCF, and Parsa Wildlife Reserve were observed to have 12, 9, 12, and 17 predominant species, respectively (Appendix D, Tables 710). In all the BZCFs, our results indicated that IVI of S. robusta was linearly decreasing (significant; p < 0.05 in all BZCFs) with moving farther from the villages, and consequently other species showed increasing IVI values with the distance. But the control forest did not show such an increasing or decreasing trend (Figure 8). The high IVI of S. robusta near the human settlement areas were due to the protection of S. robusta and thus increasing all the three parameters affecting IVI. Based on the variation of Morisita’s dispersion index (Iδ), all these predominant species exhibited clumped distribution (Appendix E, Tables 1114). No effect on the distribution pattern of species according to the distance from the settlements was observed (all species exhibited Iδ > 1). The clumped distribution of the species reflects that the distance between neighboring individuals is minimized for each species.

3.2. Species Diversity along the Distance from the Human Settlements

To examine the variations in the species diversity of woody plants, Shannon-Wiener index (H’) and Simpson’s index (λ) were calculated for each plot and plotted against the distance from the village boundaries. The Shannon-Wiener index for the species diversity increased significantly (p < 0.05) for all the BZCFs (Figures 9(a), 9(b), and 9(c), respectively, for Janajagaran BZCF, Musharni Mai BZCF, and Radha Krishna BZCF) towards the forest edge to the interior. For PWR, this pattern was not observed (Figure 9(d)). At the nearer distance to the human settlements in BZCFs, all the sites showed very low diversity indicating a high dominance of S. robusta (the species receiving high management and restoration priority). For PWR, the index values were comparatively higher. Simpson’s index of diversity, which ranges from 0 (highest diversity) to 1 (no diversity), also showed the similar pattern of increase when moving farther from the village boundary lines (Figure 10, 1 – Simpson’s λ; p < 0.01 for BZCFs but insignificant for the control forest). This trend also indicates the conversion of natural mixed forests to the woodlands near the edges.

3.3. Community Similarity and Dissimilarity between the Sites

To demonstrate how the forest disturbances change the composition of the plant community as the distance gets closer to the human settlement areas, we determined similarity and dissimilarity between the sites. For Janajagaran BZCF, the maximum value of dissimilarity was 0.55 (similarity: 1 - 0.55 = 0.45) and the dissimilarity was between Site I versus Site VI and Site II versus Site VI. For Musharni Mai BZCF, the maximum value of dissimilarity was observed between Site I and Site VI (0.54) followed by Site I and Site V (0.42). We observed similar trends for Radha Krishna BZCF with the highest dissimilarity between Site I and Site VI (0.59) followed by Site II and Site VI (0.49). Values showed the higher dissimilarity between the nearer plots and the farther plots in BZCFs. For PWR, the maximum dissimilarity was observed between Site II and Site VI, but the dissimilarity values across sites did not exceed 0.34. The results indicated that dissimilarity values and their clustering heights were higher between the sites of BZCFs than the control forest sites in PWR (Figure 11, Appendix F, Tables 1518).

NMDS was conducted to ordinate the sites based on the Bray-Curtis distance of IVI values of the species. The NMDS results showed that the sites of BZCFs were ordinated in the farther distances as compared to PWR forest (Figure 12). In BZCFs, Site I and Site II ordinated in the farthest distances from Site VI reflecting a high variation in the species composition due to varying gradients of anthropogenic pressures. However, in case of PWR, the sites were clumped (ordinated in the shorter distances) in the central region of the ordinated axes. The ordination distances, however, do not correspond to the original distances among the sites.

4. Discussion

In this paper, we assessed how the changes in the ecological parameters in the S. robusta mixed forest are a function of the distance from the human settlements, and the interaction between anthropogenic disturbances and the forest management. The changes in the structure and composition were determined by comparing their values from the edge to the interior sites of the BZCFs and also with the control forest conditions in PWR.

We did not observe a significant pattern between total stem density and the distance from the settlements, but when density was segregated into species, the sum density of species other than S. robusta witnessed a significant positive pattern (Figure 3). Human activities such as collection of fuelwood (leading to an absence of coarse woody debris in BZCFs that was present in the forest stands of PWR), fodder, litter, and other forest products, as well as grazing animals, had reduced the density of woody species in the forests. But the decrease had been compensated (regardless of species-wise contribution) by the forest management interventions in the BZCFs. Therefore, in overall, we observed neither an increase nor a decrease in the total stem density while moving towards the forest interior from the edges. Changes in the spatial and temporal heterogeneity as a result of anthropogenic activities have also been assessed in several other studies [811, 38].

Considering the total basal area, we observed mixed results in the three forests while moving from the edge to the interior. Basal area was observed with a significant decreasing trend for Janajagaran BZCF, an insignificant decreasing trend for Musharni Mai BZCF, and an increasing trend for Radha Krishna BZCF (Figure 6) along the edge to the interior of the forests. Previous research also demonstrated similar results of no effect on the basal area with a distance from the human settlements [17]. In contrast, another study has shown that biomass (a parameter indicating basal area) is reduced in the forest edge relative to the interior [39]. Consistent with this result, the basal area of the species other than S. robusta increased with the distance from the forest edges (Figure 6). Based on % basal area contributed, the forest edges were observed to have undergone biotic homogenization (S. robusta mixed stand to S. robusta stand) due to the disturbance and management interactions (Figure 6). The Ministry of Forests and Soil Conservation (MFSC) of Nepal has categorized the Sal (S. robusta) forest as a forest where S. robusta comprises more than 60% of the basal area and Sal Terai Mixed Hardwood Forest (STMH) as a forest where S. robusta comprises 33-60% of the total basal area [40]. In case of BZCFs, the edges had more than 60% S. robusta basal area, a percentage not observed in the interior stands of BZCFs and the control forest; in the latter forests the percent contribution of S. robusta was between 33 and 60% (Figure 7).

In our study, both Shannon-Wiener and 1 – Simpson’s λ diversity values showed a positive trend with the distance from the boundary lines of villages in BZCFs (Figures 9 and 10). But in PWR such a trend was not observed. Some primary studies have determined a clear pattern of change in diversity, richness, density, and basal area of the tree species towards the forest interior from the edges [12, 20, 41], whereas other studies have obtained contrasting results with more species richness and abundancy of heliophilic species in the forest edge than the interior [4244].

The increasing value of species diversity with the distance from the settlements in our study strongly supports that human interventions (interaction of both management and disturbances) have fostered a decrease in the richness and evenness of plant community and have changed the S. robusta mixed stands to pure S. robusta stands near the forest edges due to overexploitation of other plant species. Consistent with observations of intense monoculture-driven forest management nearby human settlements, other researchers have also posited that intensification of silviculture threatens species diversity [45, 46] and leads to the development of a community dominated by a single species or even-aged stands [47]. Anthropogenic pressures which lead to changes in the biological diversity are also responsible for changing the composition and stand quality of the forest ecosystems, microclimate, and nutrient cycling [48]. Increase in the generalist woody species can be the site of food sources, nest sites, and roosting sites for a variety of other species [49]. Higher plant diversity enhances the stability of the physical, chemical, and biotic properties of soils [50]. The intermediate disturbance hypothesis explains that diversity is maximum at the intermediate level of disturbance [51]. In contrast, species diversity declined with increasing level of disturbance in our study (Figures 9 and 10). There are also views reflecting that stability increases with the complexity of the ecosystems [52, 53]. In addition, the similarity analysis between sites in BZCFs reflected that there is a very little similarity between the sites in the edges compared to the interior (Figure 11). The result was also supported by NMDS outcomes (Figure 12). In the case of the control forest sites, such variation in the similarity was not observed.

Our study, based on the structural attributes, showed loosely defined community assemblages in the edges of BZCFs primarily due to the anthropogenic activities. In the long run, the existing inherent disturbances and ongoing management in the Terai S. robusta mixed forests can potentially lead to a decrease in the stability of the ecosystems. Natural forest areas that have been disturbed or altered to a large spatial extent or degree may not reach the similar structural complex (richness and composition) as primary forest [54]. In this context, fostering heterogeneity of the woody plant species is essential for sustainability.

5. Conclusions

Our results clearly indicated the effects of distance from human settlements on the vegetative attributes in the BZCFs. Maximum disturbance on the structural characteristics was observed near the edges which gradually decreased along the interior of the forests. The stands near the edges were under active management; however, a decrease in the richness and species diversity of the woody species was observed when moving from the edge to the interior. The management interventions in the forest sites near the settlements have directed the stands of BZCFs towards biotic homogenization and densification. These conditions cannot support habitat requirements of multiple species, and if the homogenization continues towards the interior of the BZCFs, it ultimately will impact the species diversity of the landscape. These situations of ecological changes due to intervention between anthropogenic disturbances and management practices should be assessed in larger spatial and regional scales to frame the future forest management strategies.


A. Predominant Species and Their Local Uses, and Geographical Location of the Sampled Sites

See Tables 1 and 2.

B. Species Accumulation Curves

The species accumulation curves for the BZCFs and Parsa Wildlife Reserve (control plot). The curves depend on the accumulation rates of new species over the increase in the sampling area. Comparatively more steep beginning rise of the curves showed that new species are increasing with addition of the sampling units. The accumulation rates of the new species in all forests were horizontal after the sampling area reached 1 ha (10 subplots each of 0.1 ha). This represents that 10, 0.1 ha plots are sufficient to represent the forest community at one site. See Figure 13.

C. Community Structure of the Studied Forests

See Tables 3, 4, 5, and 6.

D. Importance Value Index (IVI) for Woody Plant Species in Each Site in the Studied Forests. The Total Value of IVI Is 300 Which Is the Sum of Relative Frequency, Relative Density, and Relative Basal Area Values

See Tables 7, 8, 9, and 10.

E. IVI of Predominant Species and Morisita Index of Dispersion

Predominant woody species and their Importance Value Index (IVI) for each site. Total of 12 species were found predominant in Janajagaran BZCF, 10 species in Musharni Mai BZCF, 11 species in Radha Krishna BZCF, and 17 species in Parsa Wildlife Reserve. Predominant species are considered high IVIs (IVI > 10). The values in parenthesis show Morisita’s index (Iδ) (a value of 1 indicates a random distribution, values >1 indicate clumped distribution, and values < 1 indicate a uniform or regular distribution). Chi-squared based probability shows that all the species showed clumped distribution (at critical level of 0.05) according to the distance from the boundary lines of villages (null hypothesis describes that the species have random distribution). See Tables 11, 12, 13, and 14.

F. Similarity and Dissimilarity among Different Sites of the Studied Forests

In the tables, the upper right section and the lower left section show similarity and dissimilarity values, respectively, between sites. See Tables 15, 16, 17, and 18.

Data Availability

The primary data used to support the findings of this study are included within the Appendix.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The authors are especially thankful to the forest user groups and administrations of Janajagaran BZCF, Musharni Mai BZCF, and Radha Krishna BZCF. The authors would also like to express their gratitude to the Department of National Parks and Wildlife Conservation and Office of Parsa Wildlife Reserve (recently designated as Parsa National Park), Nepal, for granting permission (letter reference numbers ECO335 – 2071/2072 and ECO074 – 2073/2074) to study the forests. They are thankful to National Herbarium and Plant Laboratories, Nepal, and Central Department of Botany, Tribhuvan University for the identification of plant specimens. The technical and field support from Wyoming Reclamation and Restoration Center (WRRC), Amrit Campus, and Central Department of Environmental Science, Tribhuvan University, Nepal, is also thankfully acknowledged.