Landscape Heterogeneity Drive Plant Assemblage Dynamics and Invasibility of Semi‐Natural Grasslands Under the Long‐Term Invasion of Ageratina adenophora

1. Introduction

The grasslands, cover approximately 26% of the global terrestrial area, playing a crucial role in ecosystem service of maintaining biodiversity, carbon sequestration, climate regulation, and biogeochemical cycles [1,2,3,4]. In China, grassland occupy approximately 40 % of the national land area, equivalent to more than 263 million hectares [5]. Semi-natural grasslands, serve as habitats for indigenous herbaceous plants, are characterized as either unimproved or grazed herd pastures that account for 80.59 % of the national grassland area and are predominantly distributed across Tibet, Inner Mongolia, Xinjiang, Qinghai, Gansu, Sichuan, and Yunnan Province [6]. It has been widely recognized that the functional services of ecological barriers offered by grassland play a positive role in ecological conservation and defense against alien plant invasion [7,8,9]. However, both theoretical understanding and practical implementation remain insufficient, constraining the effectiveness of related projects and policies [10].

Yunnan Province, one of the global biodiversity hot spots, located in southwestern China and bordering South and Southeast Asian countries, contains semi-natural grasslands that not only provide pasture and habitat for livestock and diverse wildlife but also serve as ecological barriers against alien invasive organisms migrating from Vietnam, Laos, and Myanmar [11,12]. However, with the rapid growth of China-ASEAN international travel and merchandise importation, the potential risk of invasive alien pests unintentionally introduced into Yunnan and inland increased.

Ageratina adenophora, one of the most aggressive alien Asteraceae species, has raised serious concern among both the public and the academic community worldwide [13,14,15]. In China, A. adenophora was first introduced from Myanmar to Lincang, a prefecture of Yunnan Province, in the 1940s, and rapidly spread inland along the southwestern monsoon [16]. The impact of A. adenophora invasion on plant communities exhibits a complex process and uncertain consequences. Generally, following the alien invasion exacerbate, it leads to a decline in the diversity of native plant species [17,18]. However, in certain cases, the invasion of A. adenophora can also exert a positive effect by increasing the species richness within the native range [19,20]. Consequently, the question arose: besides of geographical proximity and climatic similarity, how the landscape and invasion stage drive local flora composition dynamics? Thus, considering the contradictory perspectives and the insufficient researches focusing on semi-natural grasslands, a multidimensional study addressing the long-term invasion impact of A. adenophora on plant species diversity of grassland ecosystems is required.

It is widely recognized that maintaining the level of species diversity is crucial for sustaining ecosystem stability and integrity [21,22,23]. In addition, species diversity represents a fundamental characteristic of grasslands and an important indicator for assessing the functions of grassland ecosystems [24,25,26]. Therefore, assessing plant species diversity in grasslands facilitates a comprehensive understanding of floral community dynamics and structural attributes [27,28]. Although numerous government-initiated programs and related research have recently achieved considerable progress in managing grasslands alien species, studies focusing on diverse tropic/subtropic semi-natural grasslands addressing diversity loss in long-term invasion scenario, remain comparatively limited. In addition, the invasibiliy and resilience of diverse grassland types to A. adenophora invasion is still not well described yet. Hence, to address this research gap concerning the long-term impact of A. Adenophora on grasslands, the present study, conducted in the semi-natural grasslands of Chengjiang County in central Yunnan, integrates field surveys and laboratory bioassays to: 1) clarify the invasibility of diverse grassland landscapes under long-term A. adenophora invasion; 2) determine the phytotoxic effects of A. adenophora on coexisting indigenous plant species; and 3) assess the feasibility and propose biological solutions using indigenous plant replacement for A. adenophora management.

2. Material and Methods

2.1. Study Area

Chengjiang County, located in central Yunnan Province, China, at geographic coordinates 24° 29′−24° 55′ N and 102° 47′−103° 04′ E, has a mountainous administrative area of 773 km². It is renowned for Fuxian Lake, the largest national deep freshwater lake, which plays an important role in local ecological preservation and regional climate regulation.

The county includes 3,600 hm² of grassland, encompassing four grassland types: tropical grasslands (TG), tropical shrub–grasslands (TS), warm-temperate grasslands (WG), and warm-temperate shrub–grasslands (WS) (Figure 1; Table S1). The region experiences a monsoon climate characteristic of the central subtropical plateau, with an annual precipitation ranging from 872 to 1028 mm (Figure 2a) and annual temperature range of 11.14–17.35 ℃ (Figure 2b). It is a typical hilly area in the central Yunnan Plateau, with a maximum slope of 65° (Figure 2c) and predominantly east–northwest slope orientation (Figure 2d). The highest elevation, 2802 m, is found at Mount Liangwangshan, while the lowest elevation, 1282 m, occurs at the confluence of the Nanpanjiang and Haikou Rivers.

2.2. Investigation Methods and Data Collection

2.2.1. Plot Design and Investigation

The survey was conducted in the Chengjiang grasslands between June and September 2023. Based on the grassland types (TG, TS, WG, and WS) and corresponding acreage，pairwise plots- invasion plot and corresponding non–invasion plot belonging to identical grassland types with similar ecological characteristics (elevation, slope, and slope aspect) were sequentially investigated. Specifically, in each 0.5 hm² plot, the simple random sampling method was applied to establish five small quadrats (1 × 1 m) [29]. The invasion plot area covered approximately 0.5% of each grassland type to ensure maximum representativeness. During the survey, the abundance, height, and coverage of plant species were documented, along with ecological indicators, longitude, latitude, and elevation. Plant height was measured using a scale with an accuracy of 1 mm, whereas abundance was visually assessed by counting the above–ground parts of individual plants.

2.2.2. Meteorologic and Geographic Data Sources

Climate variables were obtained from the World Climate Database (https://www.worldclim.org/), extracting two variables: annual mean temperature and annual precipitation. The Spatial Analyst tool in ArcGIS 10.8 was used for sampling and extracting climate data. Topographic variables were derived from the Geospatial Data Cloud (http://www.gscloud.cn/), where DEM data, slope, and aspect analyses were performed using the 3D Analyst tool in ArcGIS 10.8.

2.3. Effect of A. Adenophora Aqueous Extract on Indigenous Plants’ Seedling Growth

Herbaceous species with high occurrence frequency found in both invasion and non–invasion plots were selected as recipient plants against A. Adenophora extracts stress, including Imperata cylindrica, Saccharum arundinaceum, Eragrostis ferruginea, Calamagrostis epigeios, and Rumex hastatus. Plants of A. adenophora and seeds of the recipient plants were collected from the invasion plots. Fresh leaves of A. adenophora were shade–dried, ground, and sieved through a 40–mesh sieve. An aqueous extract was prepared by soaking 5 g of leaf tissue in 100 mL of distilled water for 48 hr, followed by filtration through two layers of gauze to obtain a 50 g/L extract. Extract concentrations of 25, 14, and 7 g/L were prepared by dilution with distilled water, while distilled water served as the control. Each concentration was dissolved in agar at 10 g/L with heating. Then, 20 mL of the agar solution was poured into 9-cm-diameter petri dishes containing five layers of filter paper to create agar–based germination substrates. The seeds of herbaceous plants were surface–sterilized in 0.1% (w/v) NaClO solution for 10 min, rinsed thoroughly with distilled water three times, and air–dried. Fifty seeds of each herbaceous species were placed on filter paper for the germination assay, with five replicates per extract concentration, including the control. The seeds were then cultivated in an intelligent artificial climate chamber (PRX–1200B, Ningbo Saifu Experimental Instrument) under 16 hr of illumination, 8 hr of darkness, at a temperature of 25 ± 1 ℃, and relative humidity of 90% ± 0.8 (%RH). Germination counts were recorded daily for each treatment. On the 20th day of cultivation, 30 seedlings per treatment were randomly selected for root length and shoot height measurements. An additional 20 seedlings per treatment were harvested to determine fresh weight. Germination and seedling vigor were calculated based on the following formulas [30]:

$$Germination energy\left(GE,\mathrm{\%}\right)={\displaystyle \frac{N_3}{N_{total}}}\times 100\mathrm{\%}$$

(1)

$$Germination rate\left(GR,\mathrm{\%}\right)={\displaystyle \frac{N_{germ}}{N_{total}}}\times 100\mathrm{\%}$$

(2)

$$Germination index\left(GI\right)=\sum \left(G_t\times D_t-1\right)$$

(3)

$$Vigor index\left(VI\right)=W\times GI$$

(4)

Where N₃ is the number of seeds germinated on the 3^rd day; N_total is the total seeds tested; N_germ is the cumulative number of germinated seeds; G_t is the number of seeds germinated on t day; D_t is the number of days to germination; W is the average fresh weight (g) of 20 seedlings.

2.4. Data Analysis

The Kernel Density Estimation (KDE) analysis was conducted using ArcGIS 10.8 to clarify and visualize the spatial distribution and relative abundance of A. adenophora in grasslands. The variation in density was measured using an established distance decay function, which was employed to explore the distribution and changing characteristics of density hot spots within spatial regions [31]. The invasion intensity index and community invasibility index were calculated to assess the invasion intensity of A. adenophora and the degree of community invasibility [32].

$$Relative abundance\left(P_i\right)={\displaystyle \frac{n_i}{N}}\times 100\mathrm{\%}$$

(5)

$$Community invasibility index\left(CII\right)=1-(Max{P}_i-P_i)$$

(6)

$$Invasion intensity index\left(III\right)=P_i/Max{P}_i$$

(7)

where n_i is the number of individuals of species i in the plot; N is the total number of all species in the individual plot; Pi is the observed relative abundance of A. adenophora in one invaded plot; MaxPi is the maximum of the relative abundance of A. adenophora among all invaded plots.

The following formula is utilized for the importance value and alpha diversity calculation [33].

$$Relative height\left({RH}_i\right)={\displaystyle \frac{Mean height of species i}{Total height of all spcies}}\times 100\mathrm{\%}$$

(8)

$$Relative density\left({RD}_i\right)={\displaystyle \frac{number of individuals of speciesi}{Total number of individuals of all species}}\times 100\mathrm{\%}$$

(9)

$$Relative coverage\left({RC}_i\right)={\displaystyle \frac{Coverage of speciesi}{Total coverage of all species}}\times 100\mathrm{\%}$$

(10)

$$Importance value\left(IV\right)={\displaystyle \frac{{RH}_i+{RD}_i+{RC}_i}{3}}$$

(11)

$$Richness index\left(S\right)=Species number in plot$$

(12)

$$Shannon Wiener index(H\prime )=-{\sum }_1^{s}Pi\mathit{ln}Pi$$

(13)

$$Pielou index\left(E\right)={\displaystyle \frac{-{\sum }_1^{s}Pi\mathit{ln}Pi}{\mathit{ln}S}}$$

(14)

$$Simpson index\left(D\right)=1-{\sum }_{i=1}^s{Pi}^2$$

(15)

The similarity index was utilized to compare differences in species abundance and community composition across various plots along latitude and longitude gradients to validate whether biodiversity distribution adheres to gradient variation patterns. The Jaccard similarity index was applied to examine differences between all plots, and species were quantified based on presence–absence data. β-diversity was also assessed using the Sorensen similarity index [34,35].

$$Jaccard similarity index={\displaystyle \frac{c}{a+b+c}}$$

(16)

$$Sorenson similarity index={\displaystyle \frac{2c}{a+b+2c}}$$

(17)

where a denote the number of species recorded in invasion plots; b is the number of species in non-invasion plots; c is the number of species common to both plots.

Microsoft Excel 2019 was used for data organization and for calculating Alpha and Beta diversity. IBM SPSS Statistics 26 was utilized to perform a one-way ANOVA test, with results expressed as mean values and standard errors, as well as Pearson and Spearman correlation analyses. Origin 2021 was employed for plotting figures, while ArcGIS 10.8 was employed for extracting data and mapping the climate and topography of the survey plots, as well as for kernel density estimation (KDE) analysis.

3. Results

3.1. KDE Analysis of the Local-Scale Geographic Spatial Pattern of A. adenophora

The KDE analysis revealed that A. adenophora distributed in a mosaic and unevenly scattered pattern across six townships in the grasslands. Higher density values, ranging from 5152 to 7728, converged in the region spanning three townships, Yousuo, Jiucun, and Fenglu, which are located either in central municipal townships or rural–urban fringe areas. In contrast, the areas of Haikou, Longjie, and Yangzong, which are mainly rural and distant from downtown or major communication lines, exhibited comparatively lower density, with an approximate average kernel density value of 3560 (Figure 3a). This result confirmed that transportation convenience and anthropogenic disturbances exacerbated the propagation of A. adenophora through human-mediated seed dispersal [36].

3.2. Comparative Analysis of A. adenophora Invasion Degree in Different Grassland Types

A total of 28 invasive plots of A. adenophora were investigated. Statistically, the invasion plots were categorized as 27.14% mildly invasive, 18.57% moderately invasive, and 54.29% severely invasive plots. Among the different grassland types, TG (10 plots) and TS (14 plots) exhibited the highest proportions of severe invasion (66.00% and 60.00%, respectively). In contrast, WG and WS, with 2 plots each, were predominantly mildly invaded, at rates of 60.00% and 50.00%, respectively (Figure 3b). The results indicated that the invasion degree in both TG and TS was higher than that in WG and WS, demonstrating that TG and TS experienced greater damage caused by A. adenophora invasion.

3.3. The Community Invasibility Variation of Individual Grassland Types

The grassland community invasibility analysis revealed that the community invasibility index followed a descending order-TS > TG > WS > WG (p < 0.05) (Figure 3c). In contrast, although TS exhibited the highest average invasion intensity index, the differences among grasslands were not significant (Figure 3d). This finding indicated that, compared to WG and WS, both TS and TG were more vulnerable to A. adenophora invasion, possessing stronger community adaptability and invasion competitiveness.

3.4. Plant Community Comparison Between Invasion and Non-Invasion Grasslands

Although both invasion (I) and non-invasion (N) plots share several plant species (common species), the statistics consistently demonstrated a reduction in plant species numbers in invasion plots across all grassland types (Figure 4a). Cluster analysis showed that plant species can be divided into four complementary clades (Figure 4b, c). Specifically, for WT and WS, A. adenophora invasion significantly altered community composition, with both grassland types comprising more species in non-invasion plots than in the corresponding invasion plots. For example, in WG, 33 species were identified, 23 and 29 species were documented in the invasion plots and non-invasion plots, respectively (Figure 4a). Remarkably, three species were present in invasion plots but absent in non-invasion plots, e.g A. adenophora, Phragmites australis and Ipomoea purpurea. In contrast, nine species were present in non-invasion plots but absent in invasion plots, including Rubus parvifolius, Pennisetum alopecuroides, Malvastrum coromandelianum, Dodonaea viscosa, Lespedeza bicolor, Cynoglossum furcatum, Origanum vulgare, Lactuca tatarica, and Capillipedium parviflorum (Figure 4b). Similarly, in WS, 29 species were recorded, with 21 and 25 species in invasion and non-invasion plots, respectively (Figure 4a). Invasion plots exclusively contained four species (A. adenophora, D. viscosa, Potentilla discolor, and Clinopodium megalanthum), while eight species (Erigeron canadensis, Eragrostis pilosa, Gentiana cruciata, Pistacia weinmanniifolia, Dipsacus asper, Bothriochloa ischaemum, Argentina lineata, and Aegopodium podagraria) occurred only in non-invasion plots (Figure 4b). These results indicated that WG and WS consisted of heterogeneous plant species among inter-grassland types and that A. adenophora invasion squeezed native plants from habitats and distinctly altered community composition.

Differing from WG and WS, the TG and TS grasslands comprised more abundant plant species. The greater number of plant species in these two grassland types indicated that they possessed higher plant species richness. Figure 4c indicates that in TG, 85 species were found, with 38 and 80 species in invasion plots and non-invasion plots, respectively (Figure 4a). Exceptionally, A. adenophora, Phytolacca acinosa, Clematis yunnanensis, Cynoglossum amabile, and Artemisia selengensis existed only in invasion plots, whereas Cymbopogon goeringii, Rubus rosifolius, Oenothera biennis, Trifolium repens, Urena lobata, Lotus corniculatus and Fagopyrum dibotrys were identified only in non-invasion plots. In addition, E. ferruginea, R. hastatus, Bidens pilosa, S. arundinaceum, S. scandens, I. cylindrica, A. argyi, and 31 other species coexisted in both plot types (Figure 4c). In TS, 98 species were identified. Among them, 42 and 89 species were discovered in invasion and non-invasion plots, respectively (Figure 4a). Specifically, nine species were found only in invasion plots, such as A. adenophora, Dodonaea viscosa, Buddleja officinalis, Datura stramonium, Phytolacca acinosa, Clematis yunnanensis, Equisetum arvense, Urtica fissa, and Ipomoea purpurea. In contrast, 56 species existed only in non-invasion plots. Similar to TG, 33 species were identified in both invasion and non-invasion plots (Figure 4c).

The investigation results demonstrated distinct plant species compositions across different grassland types, reflecting the drastically diversified ecogeography in the survey area. Generally, TG and TS fostered a more diverse flora than WG and WS (Figure 4a). In all grassland types, over 20 plant species coexisted in both invasion and non-invasion plots.

3.5. Alpha Diversity

Alpha diversity analysis showed that the Richness index (S) was higher in non-invasion than in invasion plots across all grassland types. Specifically, for TG and TS, there were highly significant differences (p < 0.01) between the pairwise plots (Figure 5a). Likewise, regarding the Shannon-Wiener index (H′), all grassland types exhibited a trend similar to the Richness index between non-invasion and invasion plots, for WS, there was a significant difference (p < 0.05) between them (Figure 5b). The Simpson index (D) was higher in non-invasion plots than in invasion ones across all grassland types, indicating that A. adenophora invasion decreased species diversity within individual grassland communities. Except for WG and WS, there were highly significant difference (p < 0.01) between invasion and non-invasion plots in TG and TS (Figure 5c). The Pielou index (E) of non-invasion plots was higher than that of invasion plots, except for WG, meanwhile it exhibited a significant difference (p < 0.05) between invasion and non-invasion plots in TG, and for TS, there was a highly significant difference (p < 0.01) between them. Among these four grassland types, the difference between the Pielou indices of WG plots showed the highest variation, resulting in an opposite trend (Figure 5d). Based on these results, it was implied that both TG and TS were more vulnerable to A. adenophora invasion than WG and WS.

3.6. Beta Diversity

The invasion and non-invasion plots of TG and TS exhibited moderate dissimilarity with Jaccard similarity indices of 0.2185 and 0.2012, respectively (Table 1). WG and WS also exhibited moderate dissimilarity with Jaccard similarity indices of 0.2698 each. This indicates that the invasion of A. adenophora remarkably altered grassland community assembly, as reflected by the greater differences in species composition and the reduction in the number of commonly shared species.

Sorenson similarity analysis showed moderate dissimilarity between the invasion and no-invasion plots of TG and TS, with similarity indices of 0.3587 and 0.3350, respectively (Table 2). There was moderate dissimilarity between WG and WS, with Sorenson’s similarity indices both being 0.4250. These values were consistent with Jaccard’s similarity indices, indicating that the invasion of A. adenophora resulted in a decrease in plant species diversity within the A. adenophora-invaded grassland community.

3.7. The Regression Analysis Between A. adenophora Importance Values (IV) and Grassland Flora Alpha Diversity

The regression analysis results showed that, within the invaded plots, there was a highly significant linear relationship between A. adenophora IVs and the Simpson index (R² = 0.8922) (Figure 6c), followed by the Shannon-Wiener index (Figure 6a) (R² = 0.8226), the Richness index (Figure 6b) (R² = 0.6336), and the Pielou index (Figure 6d) (R² = 0.5958). Specifically, it demonstrated that both of Shannon-Wiener index and Pielou index inverted U-shaped correlated with IVs of A. adenophora, and once IVs increase and reach the shreshold 0.36 would decrease these two indices (Figure 6a,d). These findings revealed that A. adenophora invasion exerted a sophisticated effect on grassland community dynamics, the IVs of A. adenophora, the extent of invasion significantly shaping flora composition assemblage.

Meanwhile, A. adenophora’s importance values were significantly positively correlated with longitude (p < 0.01). In contrast, they were significantly negatively correlated with annual precipitation, with a coefficient of 0.17 (p < 0.01) (Figure 6e). These results indicated that the colonization and invasion degree of A. adenophora were determined by multiple ecogeographical factors, including altitude, average annual temperature and annual precipitation.

3.8. The Inhabitory Effect of A. adenophora Aqueous Extracts on Seedling Growth of Native Plants

For screening potential native plant species in grassland, five herbaceous species with economic value and ecological importance, I. cylindrica, S. arundinaceum, E. ferruginea, C. epigeios, and R. hastatus, identified in both invasion and no-invasion plots (Figure 4b,c), were selected for laboratory bioassay. The results showed that germination rates (Figure 7a), germination energy (Figure 7b), germination index (Figure 7c), and vigor index (Figure 7d) of the five plants decreased with increasing concentrations of A. adenophora leaf aqueous extract. Specifically, the aqueous extract exerted an inhibitory effect on seedling height and root length of all tested plants (Figure 7e,f). However, S. arundinaceum exhibited the best integrative performance across five indices: germination rate, germination energy, germination index, vigor index, and seedling height inhibition rate. These results indicated that S. arundinaceum demonstrated the highest tolerance to A. adenophora leaf aqueous allelochemicals.

The seedling growth index of indigenous plant species under A. adenophora leaf aqueous extract stress likely reflects these species’ richness and competitiveness within A. adenophora-invaded plots. Spearman correlation analysis was conducted between the seedling growth indices of the screened native plant species and corresponding field data to verify this hypothesis. Table 3 indicates that, except for Rumex hastatus and Imperata cylindrica, germination energy, vigor index, seedling height, and root length were correlated with relative density and height of S. arundinaceum, C. epigeios, and E. ferruginea (Table 3).

Combined with seedling growth indices (Figure 7; Table S2), it exhibits that A. adenophora leaf aqueous extracts imposed allelopathic detrimental effects on native plant seedling growth.

4. Discussions

4.1. The Impact of A. adenophora Invasion on the Grassland Plant Diversity

A considerable body of literature has reported that A. adenophora is capable of adapting to a diverse range of habitats, such as uncultivated land, cropland, woodland, and roadsides [17,18,37]. Generally, within three to five years of invasion, native herbaceous plants are directly or indirectly affected by competitive exclusion, leading to the exclusive community formation by A. adenophora. This process results in a decrease in the number and abundance of native species within the invaded habitats. In the present survey, all A. adenophora-invaded grassland plots exhibited a trend of decreasing plant species richness, indicating that the invasion has significantly altered the species composition and structure of the grassland ecosystem. Therefore, the declining or absent native plant species in the invaded plots were squeezed from the original grassland patches (Figure 4b,c). Alpha diversity analysis revealed that in A. adenophora-invaded plots, the Pielou index, Richness index, Shannon–Wiener index, and Simpson index were significantly lower than those of the non-invaded plots across different grassland types. Similar studies have also demonstrated that plant species index decreased significantly or very significantly in habitats severely invaded by A. adenophora [38,39]. Under conditions of abundant ecological niches in specific ecosystems, A. adenophora invasion can facilitate the colonization of indigenous plant species [40]. Similarly, in the present investigation, some indigenous plant species such as D. viscosa, P. discolor, and C. megalanthum were found in invasion plots but were absent in non-invaded ones (Figure 4b,c).

Beta diversity indices are defined to assess the degree of change in species composition within biocoenosis across environmental gradients, reflecting not only the distribution of species diversity within a region but also the relationship between species and their environment [41]. In the present research, Beta diversity analysis of plants showed that the Jaccard similarity index and Sorenson similarity index exhibited moderate dissimilarity between the invasion and non-invasion plots of TG and TS, whereas they exhibited high dissimilarity in WG and WS. In conclusion, the dissimilarity of Beta diversity across different grassland types is probably associated with the sophisticated invasion process of A. adenophora and the specific species composition of each grassland type. In addition, low species similarity within an ecosystem indicates high variation in species composition, a low number of shared species, and high variability in species turnover [42]. The present findings provide complementary evidence that A. adenophora possesses remarkable adaptive, competitive, and invasive abilities compared to co-occurring plant species in grasslands of differing ecosystem types and species richness, potentially resulting in varied consequences depending on the invaded grassland type.

4.2. The Relationship Between A. adenophora Importance Values and Flora Community Diversity

The investigation of invasion plots showed that the importance values of A. adenophora range from 0.09 to 1.00 across all grassland types. Nonlinear curve-fitting regression analysis indicated that as the importance value of A. adenophora increased and reached the threshold of 0.36, the Shannon-Wiener and Pielou indices peaked, after which the alpha diversity index gradually declined (Figure 6a and Figure 6d). This finding implied that, within a specific threshold, an A. adenophora community with a lower importance value facilitates the maintenance of indigenous flora diversity. This is consistent with previous research showing that biological invasions usually decrease biodiversity on a global scale but increase species diversity on a regional scale [43,44]. Although substantial observations have been documented in numerous contexts, the mechanisms underlying these phenomena are worth investigating through interspecific allelopathic assays between endangered grassland plant species and A. adenophora.

4.3. The Allelopathic Effect of A. adenophora on the Native Plants

The aggressive competitiveness of A. adenophora is attributed to its secretion of allelochemicals that suppress adjacent plant species in the ecosystem. It revealed that volatile compounds in its root secretions negatively affect seed germination in other plants [45]. Previous research reported that the aqueous extract of A. adenophora leaves inhibited Triticum aestivum, Lens culinaris, Pinus roxburghii, Brassica campestris, and Quercus leucotrichophora in terms of seed germination and seedling development efficiency [46,47]. Despite the inhibitory potential of A. adenophora, in the present study, both invasive and non-invasive plots commonly shared a certain number of plant species (Figure 4b,c). Therefore, it was rationally speculated that these coexisting indigenous plants can be more tolerant to the hydrosoluble allelochemicals in A. adenophora leaf tissue. One of the recipient native plant species, S. arundinaceum showed strong tolerance against A. adenophora stress. A correlation analysis between the seedling growth indices obtained from the bioassay and their growth indices (relative density, relative height, and relative coverage) indicates that the indicators of seed germination and seedling growth of native plants indirectly reflect the in situ population conditions (relative height, relative cover, and others) in grasslands. Therefore, it is rationally speculated that the use of A. adenophora-indigenous plant allelopathic bioassay results can estimate individual native plant species communities in situ within an A. adenophora-invaded ecosystem. However, more intensive research is required to provide supporting evidence.

Using indigenous plants with economic or ecological value in replacement control approaches have been validated a promising alternative for IAP management. For example, Sonneratia apetala and S. caseolaris have been utilized to control Spartina alterniflora through alternative planting while promoting the restoration of native mangrove forests [48]. In the present study, the recipient plant S. arundinaceum, characterized by conspicuous tillering potential, multipurpose use, strong tolerance to diverse eco-geographical conditions, and genetic relation to sugarcane, combined with its superior tolerance to A. adenophora leaf aqueous extracts, makes it the optimal candidate for replacement control of A. adenophora. However, long-term localized field experiments are still required for further validation.