Spatial distribution of the four invasive plants and their impact on natural communities’ dynamics across the arid and semi-arid environments in northwest Pakistan

Introduction Non-native species are globally successful invaders with negative impacts on vegetation communities’ social, economic, and ecological values. Hence, the current research was carried out to assess the spatial distribution patterns and vegetative diversity of the four non-native species in severely invaded areas of the semi-arid parts of northern Pakistan. Methods The research was conducted using data from 1065 plots spread across 165 sites. These sites represented habitats throughout Northern Province, such as farm countryside, highlands, and abandoned places in rural and urban areas. Results and discussion The communities were floristically diverse, represented by 107 plant species, and dominated mainly by annual and perennial life forms with herbaceous habits. Similarly, the floristic structure shows significant variation tested by the χ2 test (P< 0.05) for plant status, life forms, life cycle, and habitat base distribution. In addition, the diversity indices show significant variation having the highest diversity in C-III (P. hysterophorus-dominated sites) and lowest in C-IV (S. marianum-dominated sites, i.e., primarily pure communities), indicating non-native species may increase or decrease site diversity. The diversity communities were further supported by higher quantities of soil nutrients, i.e., organic percentage (2.22 ± 0.04). Altitude, soil nutrients, and texture were shown to be the environmental factors most associated with communities that non-native species had invaded. Recommendation It is recommended that relevant, additional soil and climatic parameters be integrated into species distribution models to improve our understanding of the ecological niches of different species and to make a collective approach for preserving and conserving native plant communities.


Introduction
The relationship between plant communities and their environment is a crucial subject in plant ecology (Soliveres et al., 2015) that changes community structure and biodiversity, resulting in characteristic vegetation (Baldeck et al., 2013;Paker et al., 2014).In this regard, the effects of invasive species invasion and environmental factors at various scales combine to produce differences in community compositions (Baldeck et al., 2013;Havel et al., 2015;Soliveres et al., 2015).The species are introduced intentionally to transport the species having high yields of crops or accidentally as a byproduct of another process (Fofonoff et al., 2003;Gaertner et al., 2016).The species invasion has rapidly increased due to rapid trade, transportation, tourism, industry growth, technological advancements, and higher urbanization rates (Ullah et al., 2022).These factors make it easier for alien species to invade different areas and endorse pressure on the natural ecosystem that provides ecosystem services to the world's rapidly expanding human population (Pauchard et al., 2009).Therefore, invasive species have been a significant area of study for scientists worldwide (Gaertner et al., 2016;Fan et al., 2017).
An invasive species may be able to invade and spread in a new habitat despite biotic and abiotic filters (Wilson et al., 2007).Environmental parameters of any area, notably temperature, relative humidity, day length, and geographical and topographic variables, operate as an abiotic filter for newly invading invaders (Begum et al., 2021).The filter's efficiency is determined by how closely environmental conditions correspond to the species' native and invaded habitats (Thuiller et al., 2005).Invasive species establishment influences environmental conditions and natural plant communities directly or indirectly, including soil nutrients and physical and chemical characteristics (Weidenhamer and Callaway, 2010).Invasive plants affect evapotranspiration, causing severe threats to sustainable water supply (Levine et al., 2003) and ecosystem services such as resources and energy (Eviner et al., 2012).Invasive species may have deleterious impacts on native populations via several mechanisms, including genetic effects, pathogen introduction, and habitat degradation (Lee, 2002;Alpert, 2006).
Invasive species alter soil nutrients, particularly phosphorus (P), nitrogen (N), and potassium (K), which have a significant impact on how the soil and plants interact (Sun et al., 2019;Yinglan et al., 2019).Many studies have shown a relationship between soil N, P, and K levels and the diversity of plant species, and unimodal patterns between N, P, and K have emerged as the conventional wisdom (Pausas and Austin, 2001;John et al., 2007;Merunkováand Chytry, 2012).Species richness did not significantly alter after exogenous applications of N and P to Rengen Grassland, according to (Hejcman et al., 2010).However, it is essential to note that P and K fertilizer applications significantly changed the species composition (Hejcman et al., 2010).The total NPK contents of the soil considerably encourage the establishment of non-native shrubs and plants in the Yilong Lake Basin of the Yungui Plateau in China (Cui et al., 2009).In the hilly Ili River Valley area of China, Xu et al. (2011) revealed that the levels of TNTN, TKTK, and available phosphorus (AP) in the soil had a noticeably favorable impact on the local flora variety (Xu et al., 2011).In some instances, the soil's pH may indicate how many nutrients are present, which can change the variety or composition of the flora (Ravolainen et al., 2013).Although there is no discernible relationship between the richness of shrubs and trees and soil pH, annual herbs, perennial herbs, and semi-shrub plants rise linearly with soil pH (Heikkinen and Neuvonen, 1997).In addition, invasive species alter the soil organic carbon (SOC), N, and pH, changing community structure and composition (Zuo et al., 2012;Cardoso et al., 2013).
However, how invasive plant communities react to floristic and vegetative characteristics and environmental factors is still being determined.Several studies (McNeish & McEwan, 2016;Subedi et al., 2019;Ullah et al., 2021;Ullah et al., 2022) strongly emphasize how invasive species traits related to environmental factors and some research focuses on how invasive species affect native plant ecosystems (Maron & Marler, 2007;Beaury et al., 2020).However, the interaction of non-native communities with one another has received very little attention in the study area, where mainly the research focuses on the invasion of single species (Khan et al., 2020;Ullah et al., 2021;Khan et al., 2022;Ullah et al., 2022).Therefore, the present study focuses on the interactions between four non-native species communities with the hypothesis that hypothesized that the local environmental conditions associated with each community would provide a unique filter evident in the community floristic structure and composition of the resident species along with the disturbance and increases or decreases of diversity.In addition, the variation in communities across the spatial and soil variables was also questioned.The four non-native communities, i.e., Datura innoxia (Mill.),Xanthium strumarium L., Parthenium hysterophorus L., and Silybum marianum L. were focused in the study because these plant species rapidly spread in the region severely affecting the native communities' diversity.The study's aims were 1) to assess four non-native species communities and their characteristic floristic variations, 2) to assess the environmental variables affecting these communities and evaluate diversity and similarity indices, and 3) to evaluate essential factors in maintaining four non-native species communities.

Study area
The study area was Pakistan's north-west province, Khyber Pakhtunkhwa (KP), as depicted in Figure 1.The province boundaries are bordered by the Himalayan, Hindukush, and Karakorum mountain ranges, having a lowlands district (Peshawar) at 327 meters above sea level (asl) and uphills at an elevation of 7708 m asl (Tirch Mir Mountain peak) (Dawood, 2017).The province's plains and mountains generate a climatic gradient varying from south to north and north to west, such as lowlands, warm climates, and cold in the northern highlands (Ali et al., 2018).Based on temperatures, the warmest month is June, having temperatures ranging from 35.01 ± 0.96°C to 18.20 ± 0.58°C, while January is the coldest month, having temperatures ranging from 14.45 ± 2.10°C to 0.83 ± 0.82°C (Shah et al., 2017), having harsher winters with freezing-point temperatures.Moreover, the relative humidity varied from 55.32 to 78.41 percent, while the yearly precipitation ranged from 379 to 743 mm (Ali et al., 2018).The area climate is vital in determining vegetation structure (Deo and Sahin, 2015).The sampled locations were between 34.59 and 34.85 N°latitudes and 360-1200 meters above sea level (Figure 1).The vegetation structure and communities are complicated and much dissected due to the essentially rough and changeable terrain (Ahmad et al., 2020;Ullah et al., 2021).Hence, the region sustained varied biodiverse flora and fauna (Xu et al., 2019).

Studied species
The study considers four different non-native species communities, which are naturalized, invasive, or exotic (Plat 1).These communities were considered invaded because the non-native species were dominant in the population, i.e., having an importance value index (IVI) of 50% or > 50% for non-native species; however, in exotic species communities, the IVI may exceed more than 70% for a single species, e.g., Silybum marianum communities, while in some cases entire homogenized population was reported.

Datura innoxia (Mill.)
Angel's trumpet (D. innoxia) is an herbaceous perennial plant of the family Solanaceae that prefers warmer climates across the globe (tropical and sub-tropical regions).The plant grows up to 1.5 m in height and has spreading branches.The leaves are simple, glabrous, roughly oval, dark green, and shallowly lobed.Datura has trumpet-shaped flowers with a pleasant scent and a range of hues from white to yellow and light to dark purple; large, solitary, hermaphrodite; insect-pollinated flowers (Srivastava and Shukla, 2016).The specie is native to South and Central America (Waoo et al., 2014)and is recently naturalized in the region, spreading along with roadside areas, walls, and arable fields, disrupting the native plant ecosystems (Abeed et al., 2022).

Xanthium strumarium L.
Rough cocklebur (X.strumarium), an annual herb of the Asteraceae family, which grows to a height of approximately 1.5 m, has a tap root and is typically reproduced by seeds (Chikuruwo et al., 2017).The visible morphological feature is oval, triangular alternating leaves, a hairy stem, and a racemose inflorescence with pistillate heads below the staminate inflorescence (Löve and Dansereau, 1959).It is a physiological short-day plant that needs 10 hours of photoperiods in the south and 7.5 hours in the north to blossom (Abdulrahman and Winstead, 1977).The specie originated from Central or South America, progressively spreading toward North America, occupying shorelines and railway tracks and dispersing widely via water and animals (Ullah et al., 2022).In KPKP, Pakistan, the specie invaded via Afghan refugee moment and pronounced impacted maize, soybean, ground nut, and cotton crops (Hashim et al., 2013).Growing as a weed, especially in soybeans and cotton, the fruits (spiny bur) are up to 2 cm long with two stouts at the end and two achenes (Hicks, 1971).As a result, the output of both soybeans and cotton is drastically reduced (up to 70%) (McWhorter and Anderson, 1976).

Parthenium hysterophorus L.
Bitterweed (P.hysterophorus), also known as white top and bitter weed, is regarded as noxious due to its large number of seeds per plant (roughly 25,000), ability to spread quickly (Haseler, 1976), fierce competition with crops (Khosla and Sobti, 1981;McFadyen, 1992;Asif et al., 2017), and potential health risks for both humans and animals (Chippendale and Panetta, 1994).The species is considered native to Central South America or the Gulf of Mexico (Jabeen et al., 2015), is considered a noxious weed of national importance, and has rapidly spread in the arid and semiarid environment of KPKP, Pakistan (Ali et al., 2018).The primary cause of the widespread and gregarious expansion of invasive exotic P. hysterophorus is human habitat disruption, which creates empty spaces that can be quickly occupied.P. hysterophorus has been identified as a possible significant weed of Pakistan's natural and agricultural ecosystems because of its ability to thrive in semi-arid climates (Adkins and Navie, 2006).

Silybum marianum L.
Milk thistle (S. marianum), of the Asteraceae family, is an annual or biennial plant widely introduced outside its native range (southeast coast of England), becoming an invasive weed in North America, Iran, Australia, and New Zealand (Abenavoli et al., 2018).It is gorgeous in flower and vegetable seed packs (Mishra et al., 2013).Once established, the milk thistle turns into a competitor for resources and water, growing in large, thick patches that shade out other plants (Berner et al., 2002).The species is considered ruderal or weedy due to its tendency to spread in patches along roads and waste areas (Gabay et al., 1994).Milk thistle became a severe weed for winter crops, including Berseem clover (Trifolium alexandrinum L.), Wheat (Triticum aestivum L.), Sugar cane (Saccharum officinarum L.), Oat (Avena sativa L.), and Barley (Hordeum vulgare L.) in KPKP, Pakistan (Marwat et al., 2006).

Phytosociological data assessment
The phytosociological assessment (Communities floristics composition and classification) was entailed in 15 districts (Figure 1), where 21 sites of D. innoxia, 45 of X. strumarium, 24 of P. hysterophorus, and 75 of S. marianum were randomly sampled to gather vegetation data.The data set consists of 21 stands and 210 plots (21 × 10 = 210 plots), 45 stands and 450 plots (45 × 10 = 450 plots), 24 stands and 240 plots (24 × 10 = 240 plots), and 75 stands and 750 plots (75 × 10 = 750 plots) phytosociological stands and plots for D. innoxia, X. strumarium, P. hysterophorus, and S. marianum, respectively that were taken for a reasonable survey of the whole research region.We designated the D. innoxia community as a community I (C-I), the X. strumarium community as community II (C-II), the P. hysterophorus community as community III (C-III), and the S. marianum community as community IV (C-IV) for simplicity and to prevent repeating the binomial.Plot sizes were kept between 3-5 m 2 , and changes up to a factor of five were considered acceptable (Morsdorf et al., 2010).Following (Curtis and Mcintosh, 1950), the Importance value index (IVI) was calculated using relative frequency, density, and cover, in each stand.A 10-m border area from the stand's edges was eliminated to minimize edge effects (Martıńez-Falcoń et al., 2018).The plant species were identified following Flora of Pakistan (Ali and Qaiser, 1986) and the naming of plant taxa according to Aeschimann et al. (2004).The number of species, cover, and density were assessed to evaluate species diversity and invasion impacts.Based on species density, the following three crucial ecological diversity indices were calculated, i.e., species richness (S), Shannon-Wiener diversity index (H′), and evenness index (E).The numeral of species in the stand was calculated to define the species richness.The following equations were used to determine the H′ and E (Equations I and II) as described by Strong ( 2016): Where pi (species proportion), i (total species), In (pi) (natural logarithm of species proportion), and In(S) (species richness natural logarithm).
Moreover, community dissimilarities and similarities indices were calculated to evaluate the floristic and vegetation differences better.The Sørensen similarity coefficient and dissimilarity coefficient were calculated using equations III and IV following Erenso et al. ( 2014

Soil analysis
The soil characteristics of each invaded stand were assessed by auger-bored 3 kg sample from the center and two opposing corners of each stand.The topsoil layer often contains a high concentration of nutrients.Therefore, samples were collected between 0 and 30 cm deep (Castilloń et al., 2015), bulked, and mixed well to obtain a homogenous soil mixture (Kamrani et al., 2011;Martıńez-Falcoń et al., 2018).In the field, the pH was measured using a digital pH meter using a soil-water solution (1:5), and electrical conductivity (ECEC) was assessed using a conductivity meter.The soil's physiochemical properties and textural characteristics (Sand, silt, and clay percentage) were determined using air-dried samples by running them through a 2-mm filter (Reis et al., 2015).The Walkley-Black method was used to determine the organic matter (Nelson and Sommers, 1996), while the micro-Kjeldahl method was used to assess nitrogen (%).Then, exchangeable potassium (K + ) and phosphorus (P 2+ ) in mg/Kg were calculated (Yadav et al., 2002).Geometric measurements of CO 2 evolution and calcium carbonate (%) were measured following (Reis et al., 2015).The online calculator (https://www.nrcs.usda.gov)was used to evaluate soil hydraulic characteristics, including field capacity (FCFC), available water (AWAW), bulk density (BD), and wilting and saturation point (WSP) by following (Yadav et al., 2002).The hydrometer technique at Agriculture Research Institute (ARI) Swat assessed textural characteristics, i.e., sand, silt, and clay percentage (Gangwar and Baskar, 2019).

Data analysis
The relative phytosociological attributes were converted into an importance value index (IVI) (Hailu, 2017) to assign each species a specific position in the four non-native communities.In addition, PC-ORD ver.6.0, were used for further categorization of vegetation communities and was assessed using Ward's agglomerative approach by selecting Euclidean distance (McCune, 1997).After the numerical categorization based on IVI, each species in each community was given a phytosociological category if a particular species might be regarded as dominant in a community and of a group in each community (Curtis and Mcintosh, 1950;Nabi et al., 2012;Khan et al., 2013;Biondi et al., 2015;Nasrullah et al., 2015).Each community's environmental and soil variables and diversity indices (Table 1) were subjected to descriptive statistics and analysis of variance (ANOVA) to assess the variation among the non-native invaded communities.The ANOVA results were further tested by Post hoc (Tukey Honestly significant difference) test to evaluate the inter-communities variations.Canonical correspondence analysis (CCA) assessed the relationship between floristic variation within invaded species-dominated vegetation and environmental factors following Ali et al. (2022).Fifteen soil parameters, four environmental variables, and three types of diversity indices were utilized in the CCA, and a Monte Carlo permutation test was utilized to examine the post hoc interpretation of the ordination axes.The ordination bi-plots showed vector representations of the essential variables (Legendre & Gallagher, 2001).The figure was drawn using Graph Pad Ver.7.0 (Barua et al., 2020).
The floristic comparisons of the four non-native communities are presented in Table 2.The highest variations in communities' comparison were presented by the status of the communities having P< 0.001 for c 2 statistics.The associated plant invasive species were dominated in C-IV, i.e., 59% (13 species), followed by C-II, having 44% weightage.In contrast, the C-I and C-III were dominated by native species having 59% and 60% weightage of the community species.In life forms, the four non-native communities were dominated by herbaceous species having significant variation with c 2 value of 24.49 and P< 0.05.Similarly, the communities show significant differences based on the life cycle duration, having a c 2 value of 7.02 and P< 0.05.The annual species dominated C-I, II, and IV, while perennial species dominated C-III.Based on the C-I, II, and III habitats, the species was primarily found in roadside areas,   followed by wasteland and arable land.However, in C-III, most species were inhabited by arable land and wasteland.The non-native species communities were categorized into four ecologically valuable groupings (C-I, C-II, C-III, and C-IV) based on dominant species in terms of IVI (Table 3).Notably, the secondmost dominant species in each of the four communities differ uniquely, except for C-I and II, which are comparable but somewhat different in IVI, i.e., 5.72 ± 0.34 and 7.91 ± 0.88, respectively, for P. hysterophorus.Moreover, D. innoxia, with a codominance of C. sativa (17.28%) and X. strumarium (5.47%), characterizes the floristic assemblage in C-I.However, invasive P. hysterophorus (5.27%) and Chenopodium album (3.45%) also contributed to this community.The C-II was dominated by invasive X. strumarium (45.33%), having co-dominant species, i.e., D. innoxia (5.14%) and P. hysterophorus (7.91%), as well as numerous additional new species from their previous and following communities (C-I, III, and IV).P. hysterophorus (46.77%) dominates the C-III based on IVI dominance (Table 3).Apart from C. dactylon (16.72%), none of the species had IVI higher than 5% in this community.Interestingly, the group IV vegetation consisted of S. marianum (76.30%) with S. nigrum (1.75%) and A. sativa (1.81%).
The four non-native communities were further categorized using Ward's agglomerative cluster into groups (Figure 3).C-I consists of three groups having P. hysterophorus co-dominant species of group I.In contrast, invasive species might be categorized as part of group II, i.e., the Datura-Parthenium-Cannabis community.The co-dominant species of group III was Medicago denticulate, also known as the Datura-Medicago group.This native species shows that group III is more stable than groups I and II.The C-II, X. strumarium invaded regions were divided into four ecologically significant groups based on different ecosystems and species composition.I was distinguished by the presence of Xanthium-Cannabis, with C. album and C. dactylon predominating.The overwhelming prevalence of invasive P. hysterophorus and C. album in group II set it apart.In contrast to group IV, this includes X. strumarium-C.sativa, group III has X. strumarium-D.innoxia.The co-dominant specie between groups I and IV was similar, but since additional distinctive species are present, they are considered separate groups.Three separate groups comprise the C-III of P. hysterophorus, with group-I being dominated by P. hysterophorus and two significant related species, namely C. sativa and C. dactylon.Similar to group I, groups II and III likewise have C. dactylon as the co-dominant species.In C-IV, most stands had a pure population of S. marianum, which was grouped as group-I in the cluster dendrogram.Nevertheless, three co-dominant species, namely Phalaris minor, Stylosanthes humilis, and Fumaria indica, were in the heavily invaded group (group II).Melilotus indicus, Solanum nigrum, and Euphorbia helioscopia were the three codominant species in partially invaded communities, i.e., group III.

Communities' diversity and environmental variables
C-I has the lowest species richness, i.e., nine species/stand, while C-III has diverse stands having higher species richness, i.e., nineteen species/stand.The Shannon-Wiener diversity index changes considerably (p< 0.05), from 2.06 ± 0.71 (C-III) to 0.74 ± 0.91 (C-IV).The species evenness index, in contrast, varies in marginal ranges (0.77 ± 0.08 (C-I)-0.37 ± 0.45 (C-IV), showing considerable variance (P< 0.05) as depicted in Figure 4. Similarly, Sorenson similarity and dissimilarity indices were calculated as presented in Table 4.The highest similarity was shown between C-I and C-II with a Sorenson similarity index of 0.55, while the lowest was between C-I and C-IV (0.13).

FIGURE 3
Ward's agglomerative cluster classifying the communities C-I to C-IV.Khan et al. 10.3389/fpls.2023.1207222Frontiers in Plant Science frontiersin.org In environmental factors, latitude, longitude, and aspect angle were the most prominent factor that differed significantly across communities (P< 0.001).The C-III lies at a higher elevation of 803.42 ± 7.05 m, whereas C-IV stands at 675.32 ± 7.12 m, showing non-significant variation.Similarly, in soil texture, C-IV had a higher percentage of silt (32.43 ± 0.11) than C-I, II, and IV, with a significant difference indicated by ANOVA (F-value 6.24 and Pvalue 0.0005).The organic matter in C-I was lower (1.29 ± 0.04%), whereas higher in C-III (2.22 ± 0.04%) than in II and IV; however, it differed non-significantly (F-value 0.24 and P-value 0.86).Likewise, C-IV had the highest CaCO 3 percentage mean value (11.51 ± 0.05) compared to the other communities.In contrast, the C-III had the highest CaCO 3 percentage (0.56 ± 0.01), indicating a significant difference demonstrated by ANOVA (F-value 58.3 and P-value 1.1×10 -).The soil of C-IV has the highest mean value of nitrogen (2.90 ± 0.13%) compared to C-I, II, and III, with a significant difference demonstrated by ANOVA (F-value 2.4 and P-value 0.05).In contrast, the soil hydraulic parameters (Wilting point, Field Capacity, Bulk density, Saturation point, and Available water) of C-I-IV vary in very narrow ranges and differ non-significantly, i.e., P > 0.05 (Table 5).

Communities-environment relationships
The CCA Axis 1 accounted for 16.3% of the variance in C-I, followed by axes 2 and 3 and 42.2% of the variation on the three axes (Table A1).The first ordination axis showed strong positive correlations with organic matter (r =0.33), the Evenness index (r =0.32), and the Shannon-Wiener diversity index (r =0.41).On the other hand, significant negative correlations exist on axis 1 for essential variables such as lime % (r = -0.28),silt percentage (r = -0.26),and species richness (r = -0.44).These factors are illustrated in the CCA biplot, represented by arrows showing their Biplot scores (Figure 4).The CCA analysis of C-II revealed that axis 1 has 5.9% of the total variations, followed by axes 2 and 3 with 5.8% and 5.6% variations (Table A2).The correlation coefficient of nutrient contents like organic matter (r = 0.37), and phosphorus (r = 0.45),

Sorenson (Czekanowski) Similarity (Ss)
Sorenson dissimilarity index Species diversity indexes of C-I to C-IV communities.Different letters in the columns (A and B) depict significant differences (P< 0.05).
were significant on axis 3.In contrast, the diversity indices, i.e., the Evenness index (r = 0.23) and the Shannon-Wiener diversity index (0.20), revealed significant positive relationships on axis 1.However, elevation gradient with a negative influence (r = -0.42) on X. strumarium-dominated vegetation on the first axis.In addition, clay contents showed a comparatively strong positive correlation on axes 2 and 3 (r = 0.41 and 0.31, respectively) and may have a substantial role in maintaining C-II (Figure 5).In C-III, the CCA axes of the vegetation data explained a cumulative variance of 34.5% (13.5%, 12%, and 9%, for axes 1, 2, and 3, respectively).The first axis revealed a strong positive association (r=-0.41)with latitude and a negative correlation (r = 0.23) with longitude, while axis 2 had the strongest positive correlation (r = 0.41) with % clay.The findings imply that the vegetation composition is governed by the topographic gradient from lowland to highland mountains with low sand particles.In contrast, the second axis had modest connections with field capacity (r = 0.33), bulk density (r = 0.33), and available water (r = -0.32),as given in Table A3.The CCA in C-IV shows heavily loaded data on axis 1, with an Eigenvalue of 0.43 and a variance of 27.3, corresponding to 98% of Pearson's correlation.The three CCA axes explained 37% variation, with axis two and axis 3 contributing 6.5 and 3.3%, respectively (Table A4).The ordination biplot showed a definite deflection and separation of fully invaded areas from severely and partly invaded sites.The biplot showed that environmental variables (aspect degree and elevation), nutrients (CaCO 3 %, potassium, and phosphorus), and soil texture (sand, silt, and clay) factors all had an impact on the communities (Figure 4).

Discussion
In the floristics spectrum, 107 plant species were recorded in the studied region, belonging to 25 families.The Asteraceae, Amaranthaceae, Solanaceae, Brassicaceae, and Fabaceae had the most species diversity, while the other families had the least.This relationship between plant families and invading regions having low species richness has been reported globally.For instance, (Qureshi et al., 2014;Seifu et al., 2017;Iqbal et al., 2020) found that the Asteraceae, Poaceae, Amranthaceae, and Fabaceae were the dominant families of non-native species-dominated communities.Even though our data concurred with these findings, differences in taxonomic ranks at the family level may be related to geography and climate (Kosanic et al., 2018).
In the floristics composition, some of the species reported were of conservation concern, e.g., Calotropis procera, Artemisia scoparia, and Indigofera gerardiana were recently reported to be rare species in the region (Ilyas, 2015;Faisal and Khan, 2022).Native species become more susceptible to extinction due to nonnative species disrupting the variety and composition of native plant ecosystems (Zimmermann, 1991).The ability of non-native species to alter the physical and chemical composition of local soils, in turn altering community diversity, is an essential factor that assesses its potential to become naturalized (Iqbal et al., 2021).The communities were dominated by annual and perennial plant species, which is related to the opposing theories put forth in the literature that 1) annual plants, in particular, have significant effects on the homogenization of communities (Qian and Guo, 2010) and 2) perennial plants encourage community disturbance and invasion (Milchunas et al., 1988;Tappeiner et al., 1991).However, given the propensity for bigger plants to displace smaller ones over time, little herbs may be disadvantaged (Timsina et al., 2011).In areas where herbs or shrubs predominate, invasion of certain plant species was dominant on the sides of highways and waste fields, as well as in agricultural fields.Similar predictions have been made by (Burleigh, 2005;Brisson et al., 2010;Richardson, 2011;Davies et al., 2013;Salgado et al., 2019) that agricultural fields, railroad lines, and road pavements would be the most vulnerable to plant invasion.Low diversity indices were recorded in the study region, possibly related to the high invasion rate at the assessed sites.Several researchers, e.g., (Seifu et al., 2017;Qureshi et al., 2019;Iqbal et al., 2020;Iqbal et al., 2021) However, the number of species in a community is little impacted by non-native naturalized species (Awasthi et al., 2012), such as in the communities of D. innoxia.Similarly, D. innoxia and Cannabis sativa are native invasive herbaceous plants with high growth rates, which may encourage a rise in native, fast-growing species due to compositional changes.
In environmental variables, elevation was a prominent factor in every community, i.e., C-I-IV, across which the vegetation characteristics vary (S Table 2-5).The species were most abundant at lower elevations, as high elevations have more diurnal variation in the climate (Zia et al., 2018).Along the elevation gradient in non-native communities, our results showed distinct and substantial variations in community composition across environmental variables.For instance, as the altitude rises, the percentages of OM, N, and available water decrease while the percentages of lime, silt, and lime increase.Elevation impacts soil characteristics (Blackburn et al., 2011), while (Kosanic et al., 2018) reported that a plant species' capacity to invade the soil indirectly or directly impacts the chemical composition and availability of nutrients.According to research, invasive and naturalized species may use the nutrients in the soil to control their establishment and reproduction (Qian and Guo, 2010;Kosanic et al., 2018).The consistency of soil nutrients may also result from alien naturalized/invasive plants, which might encourage future invasions into populated regions (Qian and Guo, 2010).Carbon, nitrogen, phosphorous, and potassium increased in invaded areas (Bekele, 2000).In nutrient-rich situations, plant growth rate and biomass production increased concurrently with changes in plant biomass (Tappeiner et al., 1991), consistent with our findings.Environmental variable alterations have been shown to assist invaders like X. strumarium (Gallardo et al., 2017) and P. hysterophorus disruption (Milchunas et al., 1988), and this may also be true for D. innoxia.A combination of soil instability and fertilizer inputs is crucial for promoting and establishing non-native plants (Salgado et al., 2019).The relationship between plant invasion and soil nutrition is well established (Davies et al., 2013) and is shown by the increased degradation and microbiological activity in P. hysterophorus-infected soils (Brisson et al., 2010).D. innoxia has experienced the same enhancement of microbiological activities as evident from nutrient-rich zones that were primarily invaded.Based on our findings, invasion and naturalization intensities may be restricted at higher elevations, as evident from decreased non-native species IVI in higher elevation sites.
Although the precise mechanism and year of the community invasion in the area are unknown, non-native species have unintentionally been introduced into several Khyber Pakhtunkhwa regions over the past two decades during the movements of Afghan refugees, flood runoff, and tourist visits for recreational purposes (Gaertner et al., 2016;Gallardo et al., 2017).These non-native species have been widely established, naturalized, and multiplied, negatively affecting biodiversity and agriculture (Marwat et al., 2010).Depending on how well they adapt to the new environment, non-native plant species in new habitats may either grow or perish (Ullah et al., 2022).It is challenging to foresee whether these would flourish and naturalize in the new setting (Sherley, 2000;Marwat et al., 2010;Gallien and Carboni, 2017).The successful, invasive plants often spread, altering vegetation and environmental factors (Zimmermann, 1991;Le Maitre et al., 2011;Ratnayake, 2014;Ullah et al., 2022) and endangering native plants and animals (Zimmermann, 1991;Qureshi et al., 2019).Our research confirmed prior assertions that these invasive non-native plants have effectively supplanted local plant species, colonized the semi-arid parts of Khyber Pakhtunkhwa, Pakistan, and established themselves as significant invasive species there.

Conclusions
The findings of the research suggest that the area is invaded by different types of non-native species, which may be naturalized (D. innoxia), invasive (X.strumarium), and exotic (P.hysterophorus and S. marianum).In these species, S. marianum form pure communities having homogenized communities.In contrast, the P. hysterophorus communities were diverse, consisting of various native and non-native species suggesting that invasive species can decrease or increase biodiversity.The research indicated that the communities in the area are at risk, and proper strategies are needed to control and overcome the spread of these non-native species.By conducting qualitative and quantitative monitoring using species inventories (seasonally) and phytosociological methodologies, as well as mapping using ground-based techniques (through map overlays or GPS) and remotely sensed images (aerial photos, high-resolution multispectral digital data), future invasions of such invasive species may be prevented.Researchers need to engage with taxonomists, ecologists, and land managers to establish a plant detection network in vulnerable areas to prevent such species from becoming invasive.The CCA analysis of each community revealed the importance of elevation gradient in spreading these invasive species.Therefore, the environmental gradient must be considered while devising strategies for controlling these invasive species.
Moreover, the life cycle of these species in most cases was also non-overlapped, indicating the disturbance of communities in the whole growing period.In addition, there is also a need to highlight the use of these invasive species in ecosystem services as bioresources for different purposes.The biomass of invasive species may also be utilized for medicinal and phytochemical purposes, and some of these species, such as X. strumarium, have recently been reported for heavy metal pollution biomonitoring and phytoremediation (Ulalh and Khan, 2022).We also advise assessing these species' morphological, phenological, phytochemical, physiological, and reproductive biology for efficient management.In addition, as suggested by Dickie et al. (2014), we argue to assess the positive and negative traits of such species be investigated in a broad, systematic framework to prevent any inconsistencies between their risks and benefits.
is the Sørensen similarity coefficient, a = common species in both stands, b = unique/uncommon species in the first stand, and c = unique/uncommon species to the second strand The dissimilarity index can be calculated following Erenso et al. (2014) D S = 1:0 − S S (Equation IV) , reported that invaded sites have low diversity indices.In non-native populations, Shannon diversity and IVI showed a declining tendency along invasion gradients.For instance, colonies of S. marianum that have been entirely or severely invaded have low diversity indices.In invaded areas, Seifu et al. (2017); Qureshi et al. (2019), and Iqbal et al. (2020) demonstrated poor diversity of vegetation in sites invaded by non-native species.

FIGURE 5
FIGURE 5Canonical correspondence analysis (CCA) biplots for the plant communities (C-I to C-IV) and their environmental factors.

TABLE 1
Environmental and soil variables and diversity indices during the analysis of communities.

TABLE 2
Comparison of floristic characteristics (Status, life forms, Life cycle, and habitat) of the four non-native communities.

TABLE 3
Importance Values (Mean ± SE) of individual species in a community type were computed in the analysis.

TABLE 3 Continued
*** represents absence of a species in particular community.

TABLE 4
Similarity and dissimilarity indices of the four communities.

TABLE 5
Environmental and Soil characteristics of the C-I to C-IV communities.