Changes in the Ecological Parameters of Mixed Forests of Sal ( Shorea robusta Gaertn . ) Are a Function of Distance from the Human Settlements

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, MusharniMai, and Janajagaran of ParsaWildlife 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 300m 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.


Introduction
Anthropogenic activities influence the ecological processes [1][2][3] and have extensive impacts on the forest ecosystems [4][5][6][7].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 [9][10][11][12].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 4 International Journal of Ecology  1: Location of studied BZCFs and PWR forest in Bara District, Central Terai, Nepal.In these BZCFs, the forest user groups are different and also vary in a number of households using forest resources.These BZCFs are the southernmost part of a large forest landscape also known as "Charkose Jhadi".People enter from the southern border of the BZCFs to collect/harvest forest resources.
(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 [17][18][19][20].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 [16][17][18][19], 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.

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.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 speciesarea curves in the forests (Appendix B, Figure 13).In all International Journal of Ecology  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.

Experimental Design and
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".

Data Analysis.
Basal area of live woody species (dominance) was calculated as Σr 2 , where r is the radius (DBH/2) for all live woody species with DBH ≥ 1.5 cm.The density of a International Journal of Ecology  Table 3: Observed richness, total density, and total basal area along the sites of Janajagaran BZCF.

Site
Observed richness (ha −1 ) D e n s i t y± SD (individuals ha −1 ) Basal area ± SD (m 2 ha −1 ) V a r i a n c eo ff r e q u e n c y I 1    Table 5: Observed richness, total density, and total basal area along the sites of Radha Krishna BZCF.

Site
Observed richness (ha −1 ) D e n s i t y± SD (individuals ha −1 ) Basal area ± SD (m 2 ha −1 ) V a r i a n c eo ff r e q u e n c y I 8  Table 6: Observed richness, total density, and total basal area along the sites of Parsa Wildlife Reserve.

Site
Observed richness (ha −1 ) D e n s i t y± SD (individuals ha −1 ) Basal area ± SD (m 2 ha −1 ) V a r i a n c eo ff r e q u e n c y I  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 [26][27][28]: where IVI j is the importance value of the jth species, n j is the frequency of jth species, N is the sum of frequencies (N = ∑ n j ), d j 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 = ∑ d j ), x j is the sum of dominance for the jth species, and X is the total dominance across all species X = ∑ x j .Since relative frequency, relative density, and relative dominance are expressed in percentage, the maximum value of IVI is equal to 300.
International Journal of Ecology Species diversity (alpha diversity) was calculated by using Shannon-Wiener index [29] and it is defined as where 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 -Σpi 2 [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 (DIS jk ) between sites j and k using BCI can be expressed as follows [35]: where y ij and y ik 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 DIS jk and the similarity (S jk = 1 -DIS jk ) 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].
International Journal of Ecology

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.81stems 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 3-6).The median  3-6).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 15 Semecarpus anacardium L.f. 21 Terminalia alata Heyne ex Roth  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 7-10).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 11-14).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.

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.

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).
International Journal of Ecology 25 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 15-18).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

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 [8][9][10][11]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 [42][43][44].
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  [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.

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.

Figure 4 :
Figure 4: Density (stems ha −1 ) of Shorea robusta (categorized as groups under saplings, poles, and trees) plotted against the distance from the village boundary lines.(a) Janajagaran BZCF, (b) Musharni Mai BZCF, (c) Radha Krishna BZCF, and (d) Parsa Wildlife Reserve.Solid lines in the figure indicate the density of the trees.Salmon colored dotted lines show the density of poles, and purple colored dotted lines show the density of saplings.The vertical lines are standard error bars.The horizontal lines are the reference density lines for the maximum number of mature trees (85 stems ha −1 ) that 1 ha plot can withstand.

0
Distance from the boundary lines of villages (m) Number of species (Ｂ； −1
Distance from the boundary lines of villages (m)Distance from the boundary lines of villages (m)

Figure 6 :
Figure 6: Dominance (basal area) of species as a function of the distance from the village boundaries.In the figure, solid circles with solid lines, hollow circles with solid lines, and hollow circle with dashed lines indicate total dominance, dominance of Shorea robusta, and dominance of other species except for Shorea robusta, respectively.(a) Janajagaran BZCF: significant (p < 0.05), (b) Musharni Mai BZCF: nonsignificant (p > 0.05), (c) Radha Krishna BZCF: significant (p < 0.001), and (d) Parsa Wildlife Reserve: nonsignificant (p > 0.05).

Figure 13 :
Figure 13: Species accumulation curves for different sites of studied forests plotted using Kindt's exact method.The curves are for Sites I, Sites III, and Sites V for Janajagaran (J), Musharni Mai (M), Radha Krishna (RK), and Parsa Wildlife Reserve (PWR) forests.The species accumulation curves reflected increase in number of species in the beginning with addition of sampling units.The vertical lines represent the standard deviations.

Table 2 :
Geographical location of sampled sites * in the studied forests.

Table 4 :
Observed richness, total density, and total basal area along the sites of Musharni Mai BZCF.

Table 9 :
IVI of species in Radha Krishna BZCF.

Table 10 :
IVI of species in Parsa Wildlife Reserve.

Table 11 :
IVI of predominant species and Morisita Index of Dispersion in Janajagaran BZCF.

Table 12 :
IVI of predominant species and Morisita Index of Dispersion in Musharni Mai BZCF.

Table 13 :
IVI of predominant species and Morisita Index of Dispersion in Radha Krishna BZCF.

Table 14 :
IVI of predominant species and Morisita Index of Dispersion in Parsa Wildlife Reserve.

Table 15 :
Index of similarity and dissimilarity among different sites for Janajagaran BZCF.

Table 18 :
Index of similarity and dissimilarity among different sites for Parsa Wildlife Reserve. ) 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 disturbances