*Article* **Spatial Pattern and Key Environmental Determinants of Vegetation in Sand Mining and Non-Mining Sites along the Panjkora River Basin**

**Kishwar Ali 1,\*, Nasrullah Khan 2, Rafi Ullah 2, Muzammil Shah 2, Muhammad Ezaz Hasan Khan 1, David Aaron Jones <sup>3</sup> and Maha Dewidar <sup>1</sup>**

	- Doha P.O. Box 24449, Qatar

**Abstract:** A specific set of environmental conditions characterizes plant species patterns and distribution on Earth. Similarly, riparian vegetation can be impacted by anthropogenic activities like mining practices involving the removal of vegetation cover, which destroys the structure and diversity of the habitat, adversely affecting the ecosystem services. In this study, we explored the role of environmental variables and biotic intervention in deriving spatial patterns and distribution of riparian vegetation at mining and non-mining sites along the most depleted Panjkora River basin in NW Pakistan. Vegetation data and its determining factors at 28 mining and non-mining sites (14 each) were sampled using 10 m <sup>×</sup> 10 m (100 m2) systematic plots at 50 m intervals along transects in a downstream direction from the upper catchments to the bottom junction with the Swat River. We recorded 186 species in both mining and non-mining sites, belonging to 70 families comprising 174 angiosperms, 3 gymnosperms, and 9 Pteridophytes. Results show that annual or perennial therophytic life forms predominated in the Panjkora River system, indicating anthropogenic disturbances. At the same time, the aggressively invasive species, such as *Xanthium strumarium* and *Cannabis sativa,* further heightened plant community disturbances. Generally, the species diversity was higher in non-mining sites and may be attributed to habitat fragmentation. Likewise, the Canonical Correspondence Analysis (CCA-ordination) revealed that geographic coordinate (i.e., latitude *r* = 0.80; longitude *r* = 0.75) and elevation (*r* = 0.95) were more meaningful predictors than soil texture (i.e., silt%, *r* = −0.30), nutrients (i.e., potassium, *r* = −0.35; phosphorus, *r* = 0.38) and soil pH (*r* = −0.50) in shaping the spatial pattern and vegetation structure. Our result implies that the present vegetation composition and spatial assemblages are due to heavy anthropogenic interventions, especially mining activities. Therefore, the heavily degraded fragile riparian system of the Panjkora River and its tributaries needed to be conserved and restored by predicting the composition of communities in response to changing climatic conditions.

**Keywords:** Panjkora River basin; Canonical Correspondence Analysis; mining activities; socioeconomic aspects

#### **1. Introduction**

Plant species are not evenly or randomly distributed throughout the world; instead, they are distributed along distinct geographical units due to specific climatic and environmental factors [1,2]. Plant species distribution is influenced by biotic and abiotic variables, including environmental changes, soil types, species migration, habitat characteristics, terrain, hydrology, geology, tectonic plate movement, etc. Due to all of these dynamics, ancient habitat topography and climate are drastically altered, mountains are raised, and

**Citation:** Ali, K.; Khan, N.; Ullah, R.; Shah, M.; Khan, M.E.H.; Jones, D.A.; Dewidar, M. Spatial Pattern and Key Environmental Determinants of Vegetation in Sand Mining and Non-Mining Sites along the Panjkora River Basin. *Land* **2022**, *11*, 1801. https://doi.org/10.3390/land 11101801

Academic Editor: Eva Papastergiadou

Received: 4 September 2022 Accepted: 12 October 2022 Published: 14 October 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

as a result, ecological factors (competition, hybridization, and pollination) significantly impact the patterns of phytogeographical distribution [3]. Its phytogeographic analysis traces specific plant species allocation, diversification, evolution, and origin [4].

Because riparian vegetation is regarded as one of the indicators to be used in management practices, ecologists and land managers are aware of the variety of ecosystem services it provides [5]. Variations in the riparian environment substantially influence the composition and diversity of plant species [6]. Plant species variety has increased significantly due to variations in the riparian ecosystem's soil physicochemical characteristics and microclimate [7]. Natural and artificial calamities drastically destroy an area's vegetation [8,9]. Over-exploitation is the primary factor causing the observed reduction in an area's vegetation, which often progresses more rapidly than the pace of species regeneration and natural restoration. Numerous reasons, such as the removal of trees for use as lumber and wood, the expansion of agricultural land, the rapid spread of invasive species, and grazing pressure, contribute to the ongoing loss of species diversity in these environments [7].

The health of the local ecosystems and the protection of the streams and accompanying fauna depend on the riparian vegetation [10]. Despite their significance, anthropoid activities have fragmented these forest corridors lowering biodiversity and several ecological processes, making riparian forests one of the most degraded ecosystems in the world [10,11]. Floristic and structural studies are necessary to understand a given community better and track the profile of habitats with unique environmental features in response to human stresses [11]. When addressing abiotic elements such as edaphic attributes, these studies allow the comprehension of how these abiotic factors, in addition to floristic composition, environmental quality, vegetation structure, and history, might impact the vegetation of these sites [12]. Moreover, the formation of subgroups in vegetative forms with specific affinities and edaphic factors influences vegetation significantly at the local scale [13].

Similarly, local elements such as soil, topography, and vegetation in the region significantly influence riparian vegetation, resulting in a varied floristic composition with species characteristic of river banks and distinct phytophysiognomies [14]. Additionally, due to regular changes, mostly brought on by variations in the inundating regime and natural and human instabilities, riparian vegetation often displays substantial variety in its species structure and geographic distribution [15]. Because of the abundance of forest species, these forests are often better off in terms of vegetation than non-riparian habitats [16].

Sand mining has become a significant source of employment for a substantial number of people worldwide and has brought about harmful anthropogenic alterations in the river basins where this activity takes place [17]. River sand and gravel have long been utilized as aggregate for building and road construction, and demand for these resources is still growing. Sand is primarily obtained in Indonesia via river or stream mining [18]. The natural balance of a stream channel may be drastically altered by excessive sand removal [19]. By disrupting the subterranean bed, the procedure impacts the land, causing water imbalances and harmful environmental issues in the affected regions [20,21]. A drop in the groundwater level reduces the amount of water accessible for home and agricultural use. People relying on agricultural land for living also lose their jobs due to mining operations destroying their farmland [22]. Sand mining negatively impacts the vegetation's structure and content and causes societal unrest [23].

Management and successful restoration efforts in susceptible areas affected by mining depend on an understanding of species composition, vegetation structure, and interactions between vegetation and edaphic components [24]. Past sand mining studies have focused on the advantages and disadvantages of sand mining methods on the environment and local people's means of subsistence in various places. For example, the studies conducted by [20–24] have focused on the mining activities and their socioeconomic importance. However, these studies lack vegetation characteristics and their relationship with the corresponding environmental factors. Therefore, a detailed study was conducted to compare

the riparian vegetation diversity between mining and non-mining communities along the Panjkora River basin, Khyber Pakhtunkhwa, Pakistan.

The significant environmental factors that operate in riparian vegetation were also evaluated using multivariate analysis. We hypothesized that mining had resulted in a deliberate change in habitat structure that is well evident from the floristic composition and vegetation structure of communities and its diversity and can be assessed by comparing it with control sites (non-mining sites). This research will add to our understanding of the region's mining communities and will provide insight into how we can better manage these communities for socioeconomic benefit and conservation.

#### **2. Materials and Methods**

#### *2.1. Study Area*

The Panjkora River basin in Khyber Pakhtunkhwa, Pakistan, originates from tributaries of the eastern Hindu Kush Mountains and was the site of the current research (Figure 1) [25]. The Panjkora River begins as a torrent from the Hindukush Mountain's enormous icecaps. It regenerates at the intersection of Barawal and Jandool streams to the east and Ushirai and Niag streams to the west. It originates in District Dir Lower, its flow travels 113 kilometres, has a catchment area of 5905 km2, and enters the Chakdara Valley before joining the Swat River near Qalangi [26]. Alternatively, it is separated from Kabul, and the Swat River basins by western and eastern mountain ranges and the extension of these mountain ranges then forms the Panjkora River basin. Geographically, the region stretches from 71◦13 8 to 72◦22 13 E longitude and 34◦39 30 to 35◦47 17 N latitude. These mountains are covered with snow and have numerous glaciers in valleys that reach above 4000 m above sea level (a.s.l.) climatically; the area under the present work has a cool to warm summer with temperatures ranging from 16 ◦C to 32 ◦C. There is a significant drop in the temperature between December and February, and is characterized by snowfall. Rainfall is reported annually and varies from 823 to 2149 mm [27].

**Figure 1.** Study area map showing sampling points across the River Panjkora in Khyber Pakhtunkhwa, Pakistan.

#### *2.2. Vegetation Sampling and Diversity Indices*

The phytosociological vegetation attributes were assessed in mining (Experimental) and non-mining (Control) sites during the field survey for two years starting in 2020 to 2021, along the Panjkora river basin by using the quadrat method following [24]. Vegetation was recorded using a 10 × 10 m quadrat for trees, 5 × 5 m for shrubs, and 1 × 1 m for herbs in the selected stands. A total of 60 stands and 600 quadrates, which include 30 mining and 30 non-mining sites, were selected for data collection and assessment. Various quantitative variables were used including density, frequency, height, cover, importance value (IV), and importance value index (IVI), employing techniques as described by [28], were recorded in each quadrate. Herb and shrub height and cover were recorded using a forestry tap, while tree diameter was recorded at breast height (dbh). The IV was calculated as the sum of the relative frequency, relative density, and relative cover [28]. A 10-m buffer zone was removed from the stand borders, as recommended by [29], to lessen the impact of edge effects. The plants recorded were identified using *Flora of Pakistan* [30,31] plant taxonomic nomenclature.

Floristic characteristics were evaluated following the *Flora of Pakistan*, and an inventory was compiled to record species details [30]. Dried plant specimens were stored on Herbarium vouchers at the Department of Botany at Malakand Khyber Pakhtunkhwa University in Pakistan. Following [32], species density was used to compute three main diversity indices: species richness (S), evenness index (E), and Shannon-Wiener diversity index (H ). The number of species present determines the species richness of a stand. The H and E indices were calculated using the following formulas:

$$E = \frac{H\nu}{\ln S} \tag{1}$$

$$H\ell = -\sum\_{i=1}^{S} pi \text{ } \ln pi \text{ } \tag{2}$$

where *pi* = proportion of the species (*i*) to a total number of species, ln = natural logarithm of pi, and lnS = natural logarithm of species richness.

#### *2.3. Environmental Variables and Soil Analysis*

A 3 kg sample was taken from two opposing corners and the middle of each site using auger borings to define the soil parameters in each location. In general, the top layer of soil is nutrient-rich. As a result, samples were collected between 0 and 30 cm deep [33], bagged, and carefully mixed to minimize variability [28,34]. A digital pH meter and an EC meter were used in the field to measure electrical conductivity and soil pH in a soil-water solution (1:5). By air-drying the samples and running them through a 2 mm filter following USDA guidelines, the physiochemical makeup and textural characteristics of the soil were ascertained [35]. The Walkley–Black technique was employed to estimate organic matter, and wet combustion with chromic acid digestion followed by dry combustion was utilized to quantify Total and Organic Carbon [36]. Micro-Kjeldahl was used to measure total nitrogen, while Yadav et al. [37] methods were used to assess accessible phosphorus (P2+) and exchangeable potassium (K+). The geometric approach was used to determine lime (calcium carbonate; percent), and geometric measurement of CO2 evolution was done following [38]. Using an online calculator (https://www.nrcs.usda.gov (accessed on 20 November 2021)), we evaluated other soil properties, including field capacity (FC), available water (AW), conductivity (s/cm), wilting and saturation point (WSP), bulk density (BD), following [39]. The soil analysis was performed at the Swat Agriculture Research Institute (SARI).

#### *2.4. Data Analysis*

From the 28 stands, floral and phytosociological data and relevant environmental factors, i.e., geographic attributes, soil texture, and nutrients, were data-banked for analysis and interpretation. Following the flowchart in Figure 2, various floristic attributes, i.e., life form, leaf form, status in Pakistan, and Rainkiaer life form, were assigned to species following already published literature conducted in areas that resemble the study area [40,41]. The relative phytosociological qualities were converted into an importance value index (IVI). After numerically classifying mining and non-mining sites, each stand's species was allocated to a phytosociological category, either mining or non-mining sites. Recent literature was consulted to determine whether a particular species may indicate a group [31,42]. The means of the environmental and topographical factors were determined using a paired *t*-test to identify variances.

**Figure 2.** Flow chart of the methodology.

Using canonical correspondence analysis (CCA) of species IVI, the link between mining and non-mining vegetation and environmental factors was investigated. In the ordination graph, the vectors of the significant variables were also shown. All quantitative vegetation and environmental data analyses were conducted using Graph Pad Version 7.0 and Microsoft Excel [43].

#### **3. Results**

#### *3.1. Floristic and Phytosociological Diversity*

A total of 186 plants, representing 70 families of angiosperms, gymnosperms, and pteridophytes, were identified in mining and non-mining sites. The families with the highest species were Asteraceae (20), Poaceae (19), Lamiaceae (9), Amaranthaceae (8), Rosaceae (6), Brassicaceae (5 each of Cyperaceae, Cucurbitaceae, Polygonaceae, and Solanaceae), and Moraceae (4 species) (Table S1). One hundred seventy-four of the total plant species reported were angiosperms, compared to 9 pteridophytes and three gymnosperms. Based on life form, the flora was dominated by herbs 134 (72%), followed by shrubs 32 (17%), while the rest were tree species, and one was a climber (Figure 3). Based on life span (LS), perennials (P) made up the majority of the flora, accounting for 116 (62%), followed by annuals (A), 63 (34%), biennials (B), 1 (>1%), annual/perennials (A/P), 4 (2%), and annual/biennial/perennial (A/B/P), 1 (>1%).

**Figure 3.** Life form (**A**) and Habit based spectrum (**B**) of the flora recorded in the mining and non-mining sites along the River Panjkora basin.

Results revealed that in the study area, dicots dominated with 140 (75.26%), followed by monocots 34 (18.27%), while the rest were either gymnosperms or Pteridophytes. The spectrum of flora was diverse based on Status in Pakistan (SP) ranging from native to exotic: Casual (CA) 1 (>1%), Cosmopolitan (CO) 1 (>1%), Cultivated (CU) 18 (10%), Exotic (EX) 1 (>1%), Invasive (I) are 14 (8%), Indigenous (ID) 2 (1%), Introduced (IN) 2 (1%), Naturalized (N) are 21 (11%), Native (NV) 109 (59%), and Weed (W) 17 (8%) as indicated in Figure 4.

**Figure 4.** Spectrum based (**A**) on seed and status in Pakistan (**B**) of the flora recorded in the mining and non-mining sites along the River Panjkora basin.

Based on Raunkiaer's life form (Figure 5), 87 (47%) of the species were therophytes (Th), 31 (17%) were hemicryptophytes (H), 18 (10%) were megaphanerophytes (Megp), 17 (9%) were nanophanerophytes (Np), 14 (7%) were geophytes (G), 8 (4%) were mesophanerophytes (Mesp), 6 (3%) were chamaephytes (Ch), 4 (2%) were microphanerophyte (Micp) and 1 was parasite (P). The leaf size spectrum was similarly diverse, ranging from microphyllous to aphyllous, including microphyll 67 (36.02%), nanophyll 50 (26.88%), mesophyll 47 (25.26%), leptophyll 10 (5.37%), macrophyll 3 (4.30%), megaphyll 3 (1.61%), and aphyllous 1 (0.53%).

**Figure 5.** Raunkiaer life form (**A**) and leaf size class spectrum (**B**) of the flora recorded in the mining and non-mining sites along the River Panjkora basin.

In the mining sites (Group I), the IVI means value ranges from 4.52 ± 0.70–1.02 ± 0.48 representing *Cynodon dactylon* and *Maytenus royleanus*, respectively. In mining sites, the dominant plants were *Cynodon dactylon (*4.52 ± 0.70), having co-dominance of invasive *Xanthium strumarium* (4.23 ± 1.12) and *Cannabis sativa* (2.7 ± 0.87), exotic *Populus nigra* with IVI of 3.41 ± 0.98), and native *Mentha longifolia* having IVI *of* 3.12 ± 0.91. In addition, some of the species have IVI between 3–1, like *Calotropis procera*, *Centaurea cyanus*, *Chrozophora tinctoria*, *Citrus limon*, *Cirsium arvense*, *Corynephorus canescens*, *Punica granatum*, *Ricinus communis*, *Solanum surattense*, *Marsilea quadrifolia*, *Ranunculus sceleratus*, *Maytenus royleanus*, *Tagetes minuta*, *Withania somnifera*. The IVI mean values of most plant species were less than 1, as presented in Table S2.

The non-mining sites (Group II) were diverse, with more plant species than the mining sites in Group I. The mean IVI of species ranges from 4.65 ± 1.28–0.09 ± 0.09 representing *Cannabis sativa* and *Zea mays*. At the same time, co-dominated by *Pinus roxburghii* with IVI of 2.86 ± 1.23, *Cynodon dactylon* having 2.8 ± 0.94 IVI, *Dodonaea viscosa* with IVI of 2.1 ± 1.43 and *Olea ferruginea* with IVI of 2.21 ± 1.12, belonging to native species (Table 1). Similarly, the mean IVI of most plant species was less than 1.

**Table 1.** Mean importance value index of dominant plant species in Group I (mining sites) and (Group II) non-mining sites.


Note: Only twelve major dominant species having an importance value (IVI) higher than 1 are selected, and their IVI is presented as Mean ± SE.

#### *3.2. Environmental Variables and Its Influence*

Table 2 shows the various environmental factors impacting non-mining and mining areas. Non-mining sites have flat to moderately sloping topography, reasonably deep soil (up to 90 cm in depth), and typically well-drained soil. The usual soil type was sandyloam, with minor variations. All vegetation had somewhat acidic soil, likely since their parental material is comparable (igneous rocks). From species-poor (mining) to species-rich (non-mining) plant types, the value of nitrogen and phosphorus rises while that of lime, potassium, and organic matter declines. According to [1], for *Dodonaea viscosa* in the same location, the available phosphorus was quite high. It was never over the 2.5 mg/kg detection level in other plant natural communities. Numerous soil characteristics, such as wilting point, electric conductivity, bulk density, field capacity, and saturation point, decreased with species-rich to species-poor vegetation types and showed statistically insignificant differences (Table 2). Still, they did not affect these or other vegetation physiognomies significantly.

**Table 2.** Descriptive value (Mean ± SD) of environmental variables affecting the vegetation sites of mining and non-mining sites.


Elevation and latitude have a strong, significant association (*p* < 0.01; *r* = 0.93), elevation and longitude (*p* < 0.01; *r* = 0.84), whereas potassium has a weak, statistically significant relationship (*p* < 0.05; *r* = −0.41) with elevation as depicted from Figure 6. However, a weak negative correlation was shown between latitude and potassium (*p* < 0.05; *r* = −0.43), sand and silt (*p* < 0.01; *r* = −0.68), and latitude and longitude (*p* < 0.01; *r* = 0.87). In contrast, a positive association between saturation point and clay particles (*p* < 0.01; *r* = 0.94), while a significant negative correlation between clay particles and bulk density (*p* < 0.01; *r* = −0.94). The soil nutrients like potassium and nitrogen significantly correlated with electrical conductivity having *p* < 0.05; *r* = −0.36 and *p* < 0.05; *r* = −0.42, respectively. Likewise, potassium was found correlated with pH and lime, having *p* < 0.05; *r* = 0.41, and *p* < 0.05; *r* = 0.36, respectively. In addition, organic matter and electrical conductivity negatively correlated with *r* = −0.47; *p* > 0.01. Moreover, an important soil hydraulic property, i.e., wilting point, was significantly correlated with electrical conductivity (*p* < 0.05; *r* = 0.36). A negative correlation between available water and bulk density, electrical conductivity, and wilting point (*p* < 0.01, *r* = −0.69; *p* < 0.01, *r* = −0.90).


**Figure 6.** Correlation heat map showing the relationship among the environmental variables of mining and non-mining sites. Note: Legends are same as that of Table 2.

#### *3.3. Species-Environment Correlation*

Canonical correspondence analysis explains the total variance of 20.9% on the ordination axes having the highest variance on axis 1 (8.9%). Similarly, the Pearson correlation on axis 1 (0.99) is higher than on axis 2 and 3. In contrast, the Kendall correlation is greater on axis three than on axis 2 and 1 (Table 3). In the axis's correlation, geographic variable i.e., elevation, latitude and longitude on axis 1 shows significant negative correlation (*r* = −0.95, *r* = −0.80, *r* = −0.75 respectively, at *p* < 0.01). Similarly, aspect angle, and pH shows significant negative correlation on axis 3 with *p* < 0.05; *r* = −0.33, and *p* < 0.01; *r* = −0.54, respectively. Similarly, soil nutrient i.e., Potassium (*p* < 0.05; *r* = −0.35) on axis 1 and Lime on axis 3 (*p* < 0.05; *r* = −0.25) shows significant correlation indicating that these factors plays important role in shaping the communities.

**Table 3.** Canonical correspondence analysis axes summary, correlation and biplot scores of the environmental variables.



**Table 3.** *Cont.*

Note: Legends are the same as that in Table 2.

Elevation, latitude, and longitude in axis 1 have negative biplot scores of −0.521, −0.464, and −0.433, respectively. The biplot shows that critical factors influencing the community structure and composition are potassium, electrical conductivity, aspect angle, pH, and silt (Figure 7).

**Figure 7.** CCA-ordination and biplot showing vegetation distribution in 30 stands along the Panjkora river basin. Note: St. (Stand number); Environmental variables legends are same as that in Table 2; Species codes are same as that in Table S2.

#### **4. Discussion**

The present research analyzed the vegetation of riparian locations throughout the Panjkora River basin and compared mining and non-mining vegetation to determine the impacts of mining on the vegetation. The qualitative analysis of flora and its composition provides an important source of information on the dynamics of ecological changes on a geographical and temporal scale [44]. The biological spectrums of flora depict the climatic conditions and possible environmental influences that might alter the structure and composition of the community [45]. The vegetation was diversified with 186 plant species and a broad floristic spectrum. Similarly, 215 species in the riparian vegetation throughout the Hindu Kush Mountains Range was recorded by [4].

In contrast, 105 species from mining sites in four districts of Khyber Pakhtunkhwa, Pakistan was reported by [46]. The overall sum of species included in the study was greater than the sum of species listed by other studies; including [47] recorded 92 species from a mid-elevation forest in the Central Himalayas of India, while [48] reported 123 species from a tropical forest in Manipur, north India. Similarly, 85 species from Dulahazara Garjan, Bangladesh by [49], and investigators in tropical forests in various regions of the globe, for example, 89 species were reported by [50], 52 species were recorded by [51], and 94 species was reported by [52]. Most plant species in this research demonstrate the transmissible distribution of vegetation, which is thought to be the fundamental characteristic of natural vegetation [53].

Native species such *C. dactylon*, *Mentha longifolia* (L.) Huds, and *Dodonaea viscosa* (L.) Jacq predominates in the communities, whereas invasive species with greater IVI include *Xanthium strumarium* L. and *Cannabis sativa* L. The species that dominated were mostly therophytes, which make up 37.5 percent of the flora, and may damage the native plants, as revealed by [54]. The communities were also dominated by perennial and annual plant species, which is in contrast to opposing views expressed in the literature that (1) annuals, in particular, have strong effects on community homogenization [55] and (2) perennials promote community disturbance and favour alien [56,57] species, which are found in riparian communities. Smaller herb species like *X. strumarium*, *M. longifolia*, *D. viscosa*, and *C. sativa* may be at a disadvantage versus tree plant species like *Populus nigra* L. and *Pinus roxburghii* Sargent, which will eventually displace those smaller herbs [58]. This displacement eventually results in disturbed communities and will pose a risk to vegetation diversity in the future.

Understanding the ecological basis of vegetation is aided by the spectrum of diverse life forms, which depicts the physiognomy of the flora and vegetation, which is the final product of all living activities coupled with the environment [59]. According to Figure 1, therophytes, with 87 species (37.5%), predominated in the region, followed by hemicryptophytes, with 16 species (10%). The prevalence of therophytes suggests that this landform experiences frequent, minor flood disturbances and terrain modification, which has limited the occurrence of woody vegetation (phanerophytes) and forced them to relocate to riparian slopes and uplands with minimal flood disturbance. Most therophytes are annual plants that reproduce in the floodplain and marshes in the next season after surviving through the adverse season via seeds or spores. The relatively high proportion of therophytes in the riparian vegetation demonstrated that the landscape was often changed by human and natural disturbances, which encourages therophytic types of vegetation [60].

Similarly, a phytosociological study of the vegetation of Hayatabad, Peshawar conducted by [61], stated that therophytes were the leading class over others, supporting our findings. We reported therophytes followed by hemicryptophytes from the area. Because they are intimately tied to the succession pattern and life forms to the geomorphology of the river, riparian plant life forms and life history strategies are crucial to preserving these ecosystems [40]. Documenting functional vegetation's life cycle methods is also important for the ecosystem's non-native vegetation.

The area's leaf size spectrum was dominated by microphyll (36.02%), followed by manophylls, and progressed gradually with a declining percentage, with only one aphyllous species. A key indicator of climate change is leaf physiology, which is strongly affected by temperature and precipitation [62]. Nanophyll was listed as the area's major leaf size class by [63], which virtually supports our conclusion. Following our results, [64] reported that nanophylls were the dominant leaf size class, followed by Microphyll, leptophylls, and mesophyll at the University of Peshawar Botanical Garden in Azhakhel, Nowshehra. Leaf size spectra of different altitudinal zones indicated that microphylls steadily expand from lower to higher zones [65]. Similar to findings when microphylls dominate the zone at greater altitudes, nanophylls were somewhat higher in the lower zone. Microphylls were the predominant leaf size, followed by nanophylls in springs and monsoon, according to [40], who studied the leaf size spectra of Ganga Chotti and Bedori Hills. This information strengthens our current analysis since it considers the same altitudinal fluctuation. Our current results are supported by [41] study of the biological properties of plants in Tehsil Takht-e-Nasrati, which noted that microphylls were the predominant leaf size with 52.8 percent of species, followed by nanophylls (19.88 percent). The floral variety of the Frontier Region was recently explained by [66] and agreed with our results.

Species richness, which is identified as the total of various species in a particular region, is considered an important characteristic of a plant community [67]. In the mining site community in this study, the Shannon–Wiener index (H') ranged from 2.14 to 2.87, while in the non-mining community, its value ranged from 2.72 to 3.39. In mining site communities, the equitability or evenness index (J) varied from 0.86 to 0.95, but in nonmining site communities, it ranged from 0.83 to 0.96. If the value of (H') is high, the community should be more diverse. The H value of a community containing only one species should be zero (0) because the value of Pi should be equal to 1 and multiplied by Pi, which should equal zero (0). When the evenness is high, the H value will also be high. So, the H value not only indicates the number of species but also shows us the abundance of the species distributions in the total number of species among the whole community. Hence, the diversity index value of this study falls within the range reported for tropical forests [48,68]. Species richness along the river basin is influenced by edaphic, biotic and climatic factors [69]. Moreover, riparian forests connect different types of vegetation and adjacent flora, and soil and hydromorphic conditions create a spatial heterogeneity that favours species richness in these environments [11,70].

The link between the different vegetation indicators and the measured environmental variables was then ascertained using CCA. A permutation approach was used to assess the correlations between the canonical axes and the explanatory matrix and the importance of each species. The proposed associations between the response and explanatory factors were evaluated by normalizing the axis scores, focusing on the unit variance, and using axes scaled to best reflect each species. The findings indicate that the elements that support these communities include soil texture, potassium, aspect angle, and electrical conductivity. In both sites, the content and organization of the community were greatly influenced by the soils. The site soils were rich in organic matter, which is associated with high soil fertility. These soils also enable water and air to pass through them, allowing roots to penetrate more easily and supplying plants with nutrients and clay aggregation stability [71] revealing that the properties of the soil might explain the distribution of vegetation at local scales. The communities were found to be dominated by certain invasive generalist species, including *C. sativa*, *X. strumarium*, *P. nigra*, and *C. dactylon*, pioneer species of disturbed vegetation [72].

#### **5. Conclusions**

Riparian vegetation is considered the health indicator of aquatic vegetation along riparian corridors. The floristic diversity was dominated by perennials and therophytes that indicated the anthropogenic pressure, which mainly contributed to mining activities in the area. These floristic indicators provide baseline information for designing strategies for conserving and managing the riparian vegetation affected by the mining activities. Additionally, identifying plant communities' particular to the various mineral zones may serve as a foundation for phytoremediation strategies for mine waste rehabilitation. The

mechanisms underlying these vegetation variations would benefit from further investigation. These mechanisms may be related to the preferential uptake or tolerance of specific soil minerals, such as heavy metals, as well as variations in other soil properties such as pH, water-holding capacity, and the availability of macronutrient elements. The current research results also show that soil characteristics affect species richness, structure, dominance, and establishment patterns and that ecological character can also affect the composition. In these study regions, soil fertility and physical attributes may be limiting element that serves as an environmental filter for establishing species. To fully comprehend the relationship between riparian vegetation and its socioeconomic effect, we recommend comparing riparian vegetation's socioeconomic impact at mining and non-mining sites with livelihood assessments of residents.

#### **6. Limitations**

The current study focuses on floristic compositions, vegetation characteristics, and their relationships with the environmental factors that play a vital role in shaping plant communities and provide baseline information for conservation and preservation. Despite the importance of the study, there are still some limitations to the current research, which can be beneficial for the development of future research in the study area. The present study lacks information about the socioeconomic impacts of the mining operations on the livelihoods of the local people, which is important aspect and, most of the time, creates hurdles in the restoration of mining vegetation. In addition, the possible restoration process and mechanism need further comprehensive evaluation for better utilization of the information in rehabilitating vegetation on the mining sites. Moreover, riparian vegetation can promote the growth of stress-resistant plants. Identifying such plants and their possible stress-resistance potential can better provide information about the utilization of mining flora for the possible remediation and accumulation of toxic pollutants.

**Supplementary Materials:** The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/land11101801/s1, Table S1: Check List of plant species along with taxonomic attributes (Status in Pakistan, family, Raunkier life form and Leaf size classes); Table S2: Showing the IVI and Mean Standard error of 186 plant species found in both mining (Group I) and non-mining sites (Group II).

**Author Contributions:** Conceptualization, by K.A. and N.K.; methodology, by R.U.; software, by M.S.; validation and Data curation, by M.E.H.K. and M.D.; formal revision and validation, by D.A.J. All authors have read and agreed to the published version of the manuscript.

**Funding:** No funding was received for this research.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All the research data are presented in the manuscript and are available.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


MDPI St. Alban-Anlage 66 4052 Basel Switzerland www.mdpi.com

MDPI Books Editorial Office E-mail: books@mdpi.com www.mdpi.com/books

Disclaimer/Publisher's Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Academic Open Access Publishing

mdpi.com ISBN 978-3-0365-9261-9