Abstract

Industrial activities have various effects on biodiversity, posing significant threats to forest ecosystems. The current study describes the species composition, taxonomic diversity, and stand structure at Bhawal Sal forest, Gazipur, Bangladesh, as they are affected by industrialization. To achieve the goal, 30 different categorized industries within the forest were considered sampling points and three distance gradient sites viz. Site-1 = Adjacent to industries (0 m), Site-2 = 160 m distance and Site-3 = 320 m distance from industries were designated as treatments. A total of 90 quadrate plots (10 m × 10 m) were taken randomly, of which 30 plots were from each site. Through forest inventory, 46 species (24 trees, 8 shrubs, 5 climbers, and 9 saplings) were recorded from three sites. The study revealed that the stand density and basal area of mature trees (257 stems ha−1 and 8.06 ± 0.60 m2·ha−1) at Site-1 were significantly lower due to diverse industrial operations than other sites. Statistically, all the biodiversity indices of mature trees; Shannon–Wiener’s index (1.72), Simpsons index (0.82), Margalef’s index (1.38), Pielou’s evenness Index (0.39) was found to be lower at proximity to industries. The lowest species richness (12) of all plants was recorded from Site-1. However, the diameter and height distributions of Site-1 comprised young (10–20 cm·dbh) to medium-sized (20.1–30 cm·dbh) trees, while the medium to large sized (>30 cm·dbh) trees was contained at Site-3 in this study. The population structure of tree species at Site-1 also showed a fluctuating curve. Overall, this study highlights that plant ecosystems and tree population structure have declined tremendously due to industrialization. Hence, the current research could be significant for developing the management framework for the disturbed deciduous forest.

1. Introduction

The tropical forest is crucial for maintaining the world’s most diverse plant communities and controlling the earth’s climate and biodiversity. However, tropical moist deciduous Sal forests in South Asia are the most disturbed forests because of the high demand for forest resources [1]. Bangladesh is a biodiversity-rich country with 2.52 million ha of forest land (hill forest, deciduous and tidal mangrove forest), accounting for 17.47% of the country’s total land area [2, 3]. Deciduous forests are one of the most prevalent forest types in Bangladesh, encompassing 0.12 million hectares of land and accounting for roughly 4.7% of the forest area [4]. The Madhupur tract consists of two national parks: Madhupur and Bhawal, and it is one of Bangladesh’s most important moist deciduous forest regions. Although tropical moist deciduous Sal (Shorea robusta) forests are rich in diverse plant communities, they are also vulnerable to deforestation and environmental degradation [57]. This Sal forest had 66% forest cover in 1985, but by 2009, it had shrunk to 45% forest cover [8].

In Bangladesh, Bhawal Sal forest is a historically significant and ecologically diverse tropical moist deciduous Sal ecosystem. It has a tremendous influence on the ecological and economic balance, but the forest’s biodiversity is currently endangered [9]. Industrialization, agricultural expansion, land grabbing, tree plantations (including invasive species), lack of management, and unauthorized logging contribute to its rapid extinction [10]. Bhawal Sal forest is in peril due to the industrial revolution and is gradually evolving into an industrial and recreational zone. Thousands of workers from various businesses live in and around the buffer zone, which is home to more than 150 different industries [11]. In 2013, about 354 illegal industrial set-ups in the protected areas seized approximately 255.3 ha of forest land [12]. Industrial effluents are wreaking havoc on the landscape, causing natural plant communities to be disrupted and depleted ecosystems, making plant development extremely difficult. As a result, ecosystem productivity and biological diversity are lost, and the area suffers from overall environmental degradation [13, 14].

Nowadays, the pollution from industrial activities significantly affects every ecosystem component. According to its management plan, industrialization harms the park’s flora, wildlife, hydrology, and quality of air [15]. In addition, a lack of forest employees encourages illicit land grabbing to establish industries within the forest. The forest floor is disturbed by industrial development activities such as infrastructure, road networks, labor settlement, labor walking, gathering, transportation systems, and other related activities within the forest region. The coppicing potential of Sal species is rapidly dwindling due to industrial waste discharge channels. Human disturbances influence diversity and the regrowth and dominance of tree species after they have been disturbed [16]. Long-term deposition of industrial effluents in the forest floor reduces soil productivity and fertility. However, industrialization and its impacts negatively influence forest dynamics and diversity.

Several researchers [1, 5, 7, 8, 1720] have studied anthropogenic disturbance and its consequences on plant species diversity in Bangladesh’s moist deciduous forest and other forest ecosystems in world perspective. Previous studies have elucidated how human disruptions affect species composition, structure, and regeneration. There have been several studies on the effects of coal mining on vegetation composition in Bangladesh and other countries [14, 2123]. However, no research has been conducted in Bangladesh, especially on the impacts of industrialization on plant communities of the forest region. However, there is still a research void, notably regarding the effects of industrial activities on forest vegetation in the Bhawal Sal forest. Therefore, it is crucial to determine how industrialization has affected forest plant communities and conduct a comparative analysis of species diversity and population structure changes. This study’s findings will assist forest policymakers in assessing the possible causes of industrialization, understanding the overall impacts on the vegetation, and providing critical information for further research and forest ecosystem conservation. We hypothesized that industrialization had negative impacts on forest plant communities and population structure for this study. Thus, this study was undertaken to address the following research questions:(1)What are the compositions of plant species in three distant sites of the Bhawal Sal forest?(2)What are the impacts of industrialization on plant communities and tree population structure along the distance gradient of the Bhawal Sal forest?

2. Materials and Methods

2.1. Study Area

The present study was conducted at three sub-divisions: Park and Baupara under the National Park Range and Bishawakuribari under the Bhawal range in Bhawal National Park (BNP) at Gazipur district in Bangladesh. Geographically, it stands at 24°02 to 24°11°N latitude and 90°21 to 90°28°E longitude (Figure 1) in the Gazipur district of the Dhaka Forest Division. It is located 40 kilometers north of Dhaka on the Dhaka–Mymensingh highway [24]. It is one of Bangladesh’s most enormous belts of tropical moist deciduous forests. Sal (Shorea robusta) is a common species in this deciduous forest [25]. All the physiochemical characteristics of Bhawal forest soils are low compared to the forest covered areas [26]. It is mixed with manganese ferrous iron ore, commonly known as “Bhawal Kankar.” The yearly rainfall is 1500 mm and temperatures range from 11.5°C to 38.5°C. The annual relative humidity (RH) is 85.2 percent, with total evaporation of 1023.5 mm [19].

2.2. Data Collection Methods

The frequency of industrial activities in this forest was evaluated subjectively by observing the location and utilizing a Google map of the area. Distance gradient analysis was done to analyze the influence of industrialization on vegetation cover and tree species composition. This study dealt with 30 different categorized industries considered as the sampling point which covered 8.47% of the total number of industries inside the forest. Each sampling point’s spatial bounds were specified and geographic coordinates were obtained through Geographical Positioning System (GPS). A total of 90 quadrate plots were taken randomly of which 30 plots from near to industries (0 m distance) which considered as Site-1. Another 30 quadrate plots were selected within a distance of 150 m from Site-1 which was classified as Site-2 and the rest 30 plots were taken from another 150 m distance from Site-2, which is considered as Site-3 (Figure 2). Disturbance status of three sites were presented at Table 1. Each quadrate plot was 10 × 10 m2 (area = 100 m2) for the inventory of mature trees ≥10 cm·dbh and saplings (5–10 cm·dbh) [8]. Within the 100 m2 plot, a subplot of 5 × 5 m2 area was used for shrubs and climbers and another subplot 2 × 2 m2 for seedling [8]. Grasses and herbs were not inventoried in this study. The girth at breast height (gbh) of all the mature trees were measured at 1.37 m from the ground. Suunto Clinometer was used to measure the total height of all mature trees. However, the known tree species were recognized in the field, while forest officials identified the unfamiliar tree species. The identification of unknown species was done with the help of https://www.worldfloraonline.org/ [27].

2.3. Analysis of Data

The collected data were analyzed to determine phytosociological characters of the vegetation and biodiversity indices following equations that were enlisted in Tables 2 and 3. The abundance-to-frequency ratio was used to explain the species’ distribution pattern [41]. If the ratio is <0.025, it indicates regular distribution, random distribution between 0.025 and 0.05, and contagious distribution if it is >0.05 [42].

The total number of plant species found at various sample sizes was used to generate species area curves [8]. PAST software program version (4.08) was used to perform rank abundance curve and a hierarchical cluster dendrogram [43]. Also, the location map was constructed using ArcGIS software (Version 10.8). One-way ANOVA test was done to compare the difference in taxonomic diversity and forest structure by Minitab 17 software.

3. Results

3.1. Impacts of Industrialization on Species Composition, Taxonomic Diversity at the Three Distant Sites in the Study Area
3.1.1. Species Composition, Stand Density, and Basal Area

The field inventory of the 0.9 ha tropical moist deciduous forest yielded a total of 46 plant species. Of the 46 plant species, 24 were trees, 8 were shrubs, 5 were climbers, and 9 were saplings (Table 4). The number of genera and families was found to be lower in Site-1 (adjacent to industries) than in Site-3 (320 m distance) due to industrial activities (Table 5). The density (257 stems ha−1) and basal area (8.06 ± 0.602 m2·ha−1) of mature tree species at Site-1 was significantly lower compared to density (713 stems·ha−1) and basal area (24.52 ± 4.06 m2·ha−1) of mature tree species at Site-3. The density of shrubs, sapling, seedling, and climbers was statistically lowest at proximity to industries but increased with distance from the industrial zone (Table 5).

3.1.2. Cluster Analysis

A hierarchical cluster dendrogram was constructed which depicted the distribution of identified number of plant individuals across all study’s sample plots (Figure 3). The data matrix contained 90 plots and 46 plant species. Three sites were differentiated from each other based on their number of individuals (Figure 3).

3.1.3. Species Area Curve

Species-area curves showed that the number of sampling plots abundantly captured the array of plant species available at three distance sites (Figure 4). In Site-3 (320 m distance to industries), the number of plant species continued to increase in up to 30 plots, but in Site-1, it increased to 15 plots slowly then it increased rapidly. The maximum number of species was captured at all the plots of Site-3 that is less affected natural forest site (Figure 4).

3.1.4. Rank Abundance Curve

A rank abundance curve, also known as a Whittaker plot, is a graph that depicts a community’s pattern of species distribution based on species richness and evenness (Figure 5). The number of distinct species on the graph represents species richness, whereas the slope of the line that fits the graph represents species evenness. The steep gradient in the Site-2 and Site-3 rank abundance curves on species abundance data indicated low evenness, as the high-ranking species have substantially larger abundances than the low-ranking species (Figure 5). The analysis revealed that Site-3 had the highest abundance compared to the other sites.

3.1.5. Species Richness, Overlapping, and Similarity Index

Industrial activities along with variety of factors contribute to species extinction. Within the whole study area, 24 mature tree species were found (Table 6). Site-1 has the fewest mature tree species (6), followed by Site-2 (13) and Site-3 (16). The lowest plant species richness (12) was found in Site-1 compared to Site-3 (37). According to the result, the greatest 15 plant species (7 mature trees, 4 shrubs, 3saplings and 1 climbers) overlapped in between Site-2 and Site-3 that followed by in between Site-1 and Site-2 where 3 is mature tree, 1 is shrubs, 3 is saplings. The highest similarity index (0.32) of all plant species was found within Site-1 and Site-2 followed by other combined sites of the study area (Table 6).

3.1.6. Diversity Indices

Industrial activities had a huge impact on the diversity indices in the study sites. In almost all cases, statistically diversity indices were the lowest at close distance to industries, but it increased gradually with increasing distance. After statistical analysis it was found that the Shannon–Wiener index, Simpson’s index, Margalef’s index, Pielou’s evenness index, and Menhinick index of mature trees of Site-1 had the significantly lowest diversity compared to other sites (Table 7). The Simpson index, Margalef’s index, Pielou’s evenness index, and Menhinick index of sapling, shrubs, climbers, and seedlings at Site-1 were comparatively lower than other two sites (Table 7).

3.1.7. Distribution Pattern of Plant Species at Three Sites

The analysis of the distribution pattern of all the recorded species is shown in Table 8. In cases of all the sites, the majority of species in the industrial area of forest were contagiously distributed (Table 8), with only a few species having regular or random distribution. However, a comparison shows that the species in Site-2 were more regularly distributed, whereas those in Site-3 were more randomly distributed. In Site-1, mostly recorded species had a contagious distribution pattern (Table 8).

3.2. Impacts of Industrialization on Structural Composition of Tree Species in Three Sites
3.2.1. Girth Class of Tree Species

The girth classes demonstrated an increasing trend with increasing tree density in all three sites at certain level (Figure 6). At adjacent to industries (Site-1), the lowest stand density (4.0 stems ha−1) was found in the girth class (40.1–45 cm) whereas the highest stand density (77 stems ha−1) was found in the girth class (60.1–65 cm) (Figure 6). The results of Site-2 revealed that tree density (stems ha−1) of tree species gradually grew up to a girth class range of 61.1–70 cm (Figure 6). The lowest stand density (3 stems ha−1) was belonged to (>80 cm) the girth class while the girth class (65.1–70 cm) had the highest stand density (314 stems ha−1). In Site-3 (320 m distance to industries), the girth class (60.1–65 cm) had the highest stand density (523 stems ha−1) with (Figure 6).

3.2.2. Diameter and Height Class of Tree Species

The diameter distributions of the sampled trees followed a fluctuating distribution in all individual plots in three sites. With an increase of the diameters at breast height (DBH), the tree density increased at certain point (Figure 7(a)). The highest percentage of trees of Site-1 (adjacent to industries) was belonged to 15.1–20 cm diameter class range. However, the highest percentage of tree species of Site-3 (320 m distance to industries) was belonged to 20.1–25 cm diameter class followed by Site-2 (Figure 7(a)). The highest percentage of comparatively larger trees in terms of diameter was found at Site-3.

The height of mature trees of near to industrial sites (Site-1) did not exceed 20 m, whereas the tree height in the Site-3 that was low industrial affected site exceeded 20 m (Figure 7(b)).

3.2.3. Relative Density, Relative Frequency, Relative Dominance, and Importance Value Index (IVI)

The importance value index of the tree species indicates the overall dominancy of different species in forest area. Swietenia macrophylla King. f (IVI = 91.06) was the dominant species in the Site-1 (adjacent to industries) that was severely industrial affected site with relative density (29.87%), relative frequency (33.33%) and relative dominance (27.85%). In Site-2 (160 m distance to industries), Shorea robusta Gaetrn (IVI = 124.45) was the dominant species with relative density (51.14%), relative frequency (21.11%) and relative dominance (52.20%), followed by Artocarpus heterophyllus Lamk, Anthocephalus cadambo (Roxb) (Table 8). Shorea robusta (IVI = 173.08) was the most dominant species of Site-3 (320 m distance to industries) with relative density (71.96%), relative frequency (32.86%), and relative dominance (68.21%), followed by Microcos paniculata L, Dillenia pentagyna Roxb, Albizia lebbek (L) Benth, and Careya arborea Roxb (Table 9).

4. Discussion

It is well accepted that industrialization in the forest zone and human disturbances reduce species richness and contribute to species extinction [44]. The total number of species in the investigated area (46 plant species under 38 genera and 27 families: 24 trees, 8 shrubs, 5 climbers, and 9 saplings) is similar to the findings of [19]. 49 plant species (24 mature trees, 4 shrubs, 7 climbers, and 14 herb species); 43 plant species (22 trees, 10 shrubs, 6 climbers, and 5 herbs) were enlisted at Bhawal National Park [25]. However, it looks to be poor in plant species diversity compared to other South Asian tropical deciduous forests. Compared to 139 tree species (100 genera and 40 families) in the Madhupur National Park, Bangladesh [26]; 135 species (105 genera of 45 families) in the North Central Eastern Ghats of India [45]; and 89 plant species (34 trees, 15 shrubs, 25 herbs, and 15 climbers) in the tropical moist deciduous forest of Assam, northeast India [46], 58 plant species were recorded the Subba et al. [47] in the Sumbuk reserve forest in South Sikkim, India, composition of the current study is relatively low. However, the species composition of Site-1 (12 plant species; 6 mature tree species, 2 shrubs, and 4 saplings) is quite low, indicating that industrial activities hurt species composition. The present result is similar to the findings of other studies, such as Rabha [48] found 18 species (14 families) in disturbed Sal forests. Borah et al. [5] reported that the lowest number of species (42 species from 26 families) in the disturbed forest; 14 tree species in the severely disturbed sites in the southern part of India [49]; only 10 species were recorded in the significantly disrupted forest in Nepalese Sal forests [1]. This may be due to the rapid industrialization, urbanization, logging, forest resource exploitation, and other disturbances that have resulted in species extinctions in diverse plant ecosystems.

Among the three sites, Site-1 (adjacent to industries) appears to have the lowest tree stand density and basal area compared to Site-3, which is 320 m away from industries (Table 4). This study’s finding raises concerns about the viability of the Sal forest ecosystem. The tree density obtained in this study’s industrially damaged site is comparable to earlier results from extremely disturbed stands (338 stems ha−1) by Bhuyan et al. [50]; (246 trees ha−1) in the disturbed tropical forests of Borah et al. (2012); and (515 stems ha−1) in the mined regions by Sarma et al. [14]. Tenzin and Hasenauer [51] found that in Bhutan’s broad-leaved forests, tree density declined from 373 N·ha−1 in the natural forest zone to 114 N·ha−1. The study found that crushing and trampling of seedlings during leaf litter collection reduced species richness and density near industries. The tree basal area (8.06 ± 0.602 m2·ha−1) of Site-1 measured in this study is comparable to values reported in previous studies, such as (14.7 ± 4.1 m2·ha−1) in the frequently used Bhawal National Park by Rahman and Vacik [25]; (18.44 m2·ha−1) in the disturbed forest in India [5]; and (12.22 m2·ha−1) in the highly disturbed forest of the Western Himalaya, India [52]. This could be due to industrial activity and waste disposal channels that harm forest regeneration. A continuous pattern type species-area curve was found in this study that is similar to Behera et al. [53]; Rahman et al. [8]; Rahman and Vacik [25]; and Yeom and Kim [54]. The rank abundance curve of the study exhibited a steep gradient similar to Ismail et al. [55], indicating a pattern of species distribution in which one species is highly abundant and the rest are moderately to relatively rare [56]. The Shannon–Wiener diversity index values for three sites’ tree species forests are all within a narrow range (1.72–2.23), which is similar to the findings of Mishra et al. [57]; Tripathi et al. [58]; Borah and Garkoti [59]; Borah et al. [60]; Borah et al. [5]; Bajpai et al. [61]; Tahir et al. [62] that fall within the range (0.56–3.57). Industrial development and human disturbances are probably responsible for the smallest value of the Shannon–Wiener diversity index in degraded forests. The proximity to an industrial site had a lower Simpson’s index (0.82) for mature trees than other sites, which was supportive of the Simpson’s index value (0.94) in severely disturbed oak forests in Iran by Eshaghi et al. [17]. The Shannon–Wiener index, as well as Simpson’s diversity and its components (richness and evenness), had a relationship with the distance gradient. In similar findings, Simpson's concentration of dominance (0.04) and similarity index (0.07) were found in the severely disturbed site of a tropical wet evergreen forest reported by Bhuyan et al. [50]. A progressive reduction of all biodiversity indices from 320 m distance to the proximity of industries revealed how industrial activities influenced natural communities. Majumdar and Datta [63] found Simpson’s dominance index (0.30 to 0.09), and Pielou’s evenness index (0.06 to 0.84) in disturbed forest in Northeast India. Eshaghi et al. [64] revealed that extensively disturbed stands had the lowest Margalef’s index (0.794), which is lower than in the current study. Awan et al. [65] found lower values of species diversity (2.88) in the Kashmir Himalayas, Pakistan due to severe grazing pressure. As a result of increased industrialization within the forest, human settlement and road networks are fast expanding, resulting in resource exploitation and encroachment, diminishing biodiversity. The majority of the plant species in this study are clumped together or spread out contagiously. The contagious distribution of the species indicates natural vegetation is fragmented due to mining industries, as reported by Sarma [23]. Nandy and Das [66] discovered 51.06% and 62.05% regular distribution patterns of plant species in Badshahitila and Loharbond, respectively, but 64.29% contagious distribution patterns in the Sultani area of natural forests in the Barak valley, Northeast India, which supports the current findings. The similarity index of these three sites for plant species (0.23, 0.14, and 0.32) is similar to the Jaccard similarity index of 0.35 (35%) in the sacred groves of Manipur, North East India, reported by Khumbongmayum et al. [67].

Anthropogenic activities such as development of industries, human settlement, construction of roads network, labor gathering, etc., triggers the changes the phytodiversity and structure of forest. These factors have impacts on biodiversity, species richness, regeneration percentage, and forest community composition. Midha and Mathur [68] mentioned several activities as the cause of fragmentation of habitats in the Terai ecosystem. The shift in dominance (based on IVI) in plant communities was the most notable change in species composition due to the rise of illegal industries, waste water discharge channels and other factors. The main species at Site-1 was Swietenia macrophylla King. f (IVI = 91.06), whereas Shorea robusta (IVI = 173.08) was the most common species at Site-3. Similar results were found in some earlier works by Rahman et al. [26]; Borah and Garkoti [59]; and Sharma et al. [69]. Mondol et al. [19] reported that based on IVI value Shorea robusta (121.584) was the dominant species in this Sal forest which supports the findings of this study.

In the present study, the largest number of tree individuals were found at 40.1–45 cm girth at a distance of 0 m from industries. However, at a distance of 320 m from industries, the highest tree density was recorded in the 60.1–65 cm girth class. However, the tree density in all the girth classes of Site-1 was low, which could be due to the rampant and random clearing of forest for industrial development purposes that has led to the drastic change in the tree population structure of the forest. The present findings are similar to earlier reports from different tropical forests of India [5962]. The highest percentage of trees was found in the 15.1–20 cm diameter range of Site-1. With an increase in the diameter at breast height (DBH), the relative abundance of trees decreased. Rahman et al. [8] reported that the diameter distributions followed a reverse J-shaped distribution in Madhupur and Bhawal National Parks, which supports this study. The height of the mature tree at Site-1 did not exceed 20 m; whereas, at Site-3, the height of the tree exceeded 20 m. Rahman and Vacik [25] found a similar finding, which is supportive of the present study. Due to industrialization, the vegetation of this park is now in an endangered condition. Within the study area, the species richness of all plant groups tended to decrease close to the industries. The construction of infrastructure for industries may degrade the forest land through erosion, sedimentation, and pollutant runoff. As a result, root oxygen uptake is reduced, resulting in decreased water and nutrient absorption and slowed water transport through soils. Camping and strolling have the most visible effects on the environment, including the crushing, shearing off, and uprooting of vegetation [70]. These effects may cause changes in the vegetation, which indicates possible reasons for lower diversity and richness in the vicinity of the industries of the present study. Massive human pressure from industry settlement resulted in a decrease in species richness and diversity in this forest. Illicit industrial set-ups trigger other disturbances such as picnic sites, logging of native species, seedling damage, leaf litter sweeping, and livestock grazing.

5. Conclusions

The significantly lower biodiversity and stand features close to industries corroborate the growing industrial activities and their harmful influence on biodiversity. The disturbance for the establishment of industries has limited the chances of regeneration of species, thereby reducing the number of species. Industrial activities are “a double-edged sword for the socioeconomic movement,” in that they are reviled and revered at the same time. It has been criticized for its repercussions and praised for fostering long-term economic growth. Implementing biodiversity conservation activities in Bhawal National Park is difficult due to its proximity to the mega-city. However, industrialization’s consequences can be mitigated in numerous ways. The forest authorities should prosecute industrialists who set up businesses illegally by grabbing land. Scientific industries have to be appropriately established to minimize further damage to the vegetation. However, additional comprehensive research activities are required to assess the effects of industrial activities on soil properties and their function in ecological processes. This study’s findings will aid in monitoring and protecting tropical moist deciduous forests in Bangladesh.

Data Availability

The authors confirm that the data supporting the findings of this study are available within the article and its supplementary material. Raw data that support the findings of this study are available from the corresponding author, upon reasonable request.

Conflicts of Interest

The authors declare that thy have no potential conflicts of interest.

Acknowledgments

The authors would like to thank the authorities of Bhawal National Park in Bangladesh for their invaluable assistance throughout data gathering and would like to extend their gratitude to the Ministry of Science and Technology, Government of the People’s Republic of Bangladesh, for awarding the National Science and Technology (NST) fellowship for patronizing financial support, and the Ministry of Science and Technology, Government of the People’s Republic of Bangladesh gave a fellowship award named National Science and Technology (NST) fellowship for (54000 BDT).