Do Habitats Show a Di ﬀ erent Invasibility Pattern by Alien Plant Species? A Test on a Wetland Protected Area

: Biological invasions are deemed to be the second most important global driver of biodiversity loss, right behind habitat destruction and fragmentation. In this study, we aimed at testing if community invasibility, deﬁned as the vulnerability to invasion of a community, could be associated with the characteristics of a given habitat, as described by the composition and structure of its native species. Based on a probabilistic sampling of the alien ﬂora occurring in the temperate wetland Lake Doberd ò (Friuli Venezia Giulia region, NE Italy) and using a null-model-based approach, the observed occurrence of Invasive Alien Species (IAS) within sampling units was randomized within habitats. While testing the degree of invasibility for each habitat within the wetland, our null hypothesis postulated that habitats are equally invaded by IAS, as IAS can spread homogeneously in the environment thanks to their plasticity in functional traits that makes them able to cope with di ﬀ erent ecological conditions. The obtained results comparing observed IAS frequencies, abundance and richness to those obtained by the null model randomizations show that, for all habitats, invasion was selective. Speciﬁcally, a marked preference for habitats with an intermediate disturbance level, a high nutrients level and a medium-high light availability was observed, while an avoidance was detected for habitats characterized by lower levels of nutrients and light availability or extreme conditions caused by prolonged submersion. This method allows us to provide useful information using a simple-to-run simulation for the management of the IAS threat within protected areas. Moreover, the method allows us to infer important ecological characteristics leading to habitat invasion without sampling the environmental characteristic of the habitats, which is an expensive operation in terms of time and money.


Introduction
Biological invasions are deemed to be the second most important global driver of biodiversity loss, right behind habitat loss and fragmentation [1][2][3]. Indeed, the introduction and spread of Invasive Alien Species (hereafter IAS) [4,5] have heavy repercussions on ecosystems, leading to negative impacts on native species, habitats, landscape, ecosystem production, naturalistic value, human health and the economy [6][7][8][9]. For these reasons, research interest on biological invasions has exponentially increased in recent years [10,11].
including rare and endemic ones, in relatively small areas; moreover, wetlands act as ecological sinks, accumulating nutrients transported by water, and, for this reason, they can also be more susceptible to IAS invasion [61,62].
In this study, we analyzed the possible effect of habitat features on community invasibility on a temperate wetland in NE Italy using null models. Specifically, we aimed at estimating the degree of invasibility for each habitat within the wetland PA. Our null hypothesis postulated that IAS presented a plasticity able to cope with differences in biological and environmental habitat characteristics, allowing them to spread without exclusive preferences; in other words, we expect them to be randomly spread across the wetland PA.

Study Area
The study was carried out at Lake Doberdò (Figure 1, Lat. 45.831574, Long. 13.562023), a karstic lake in north eastern Italy. It is part of the regional natural reserve "Laghi di Doberdò e Pietrarossa" and is included in the Natura 2000 network (SAC-IT3340006/SPA-IT3341002).
Metabolites 2020, 10, x FOR PEER REVIEW 3 of 17 accumulating nutrients transported by water, and, for this reason, they can also be more susceptible to IAS invasion [61,62].
In this study, we analyzed the possible effect of habitat features on community invasibility on a temperate wetland in NE Italy using null models. Specifically, we aimed at estimating the degree of invasibility for each habitat within the wetland PA. Our null hypothesis postulated that IAS presented a plasticity able to cope with differences in biological and environmental habitat characteristics, allowing them to spread without exclusive preferences; in other words, we expect them to be randomly spread across the wetland PA.

Study Area
The study was carried out at Lake Doberdò (Figure 1, Lat. 45.831574, Long. 13.562023), a karstic lake in north eastern Italy. It is part of the regional natural reserve "Laghi di Doberdò e Pietrarossa" and is included in the Natura 2000 network (SAC-IT3340006/SPA-IT3341002).
The study area (ca. 65 ha) includes the lake and all the surrounding vegetation directly dependent on the presence of lake water, including the meso-hygrophilous one. Bioclimate is temperate with lower meso-temperate thermotype and lower humid ombrotype, with relatively warm summers and dry and cool winters. The verage annual rainfall is around 1200 mm. The mean annual temperature is 15 °C, with an average temperature of 5 °C in January (coldest month) and 23 °C in July (warmest month) (source: https://www.meteo.fvg.it; reference period: 1999-2018).
The study area is characterized by a complex geomorphology. It consists of a thick succession of limestone and subordinately dolomitic rocks, dated between the Lower Cretaceous (Aptiano) and the lower Eocene; the platform deposits are surmounted by turbidite deposits (flysch) [63][64][65]. The soils of the study area are classified as "red soils", e.g., clayey or silty-clayey soils that cover the calcareous or dolomitic substrates with variable thicknesses, filling particular fractures and areas with depressed morphology [66]. The study area (ca. 65 ha) includes the lake and all the surrounding vegetation directly dependent on the presence of lake water, including the meso-hygrophilous one.
Bioclimate is temperate with lower meso-temperate thermotype and lower humid ombrotype, with relatively warm summers and dry and cool winters. The verage annual rainfall is around 1200 mm. The mean annual temperature is 15 • C, with an average temperature of 5 • C in January (coldest month) and 23 • C in July (warmest month) (source: https://www.meteo.fvg.it; reference period: 1999-2018).
The study area is characterized by a complex geomorphology. It consists of a thick succession of limestone and subordinately dolomitic rocks, dated between the Lower Cretaceous (Aptiano) and the lower Eocene; the platform deposits are surmounted by turbidite deposits (flysch) [63][64][65]. The soils of the study area are classified as "red soils", e.g., clayey or silty-clayey soils that cover the calcareous or dolomitic substrates with variable thicknesses, filling particular fractures and areas with depressed morphology [66]. Lake Doberdò is a karst lake with a typical peculiar hydrology: it has no superficial inflow or outflow; the superficial runoff is near zero and the water feed mainly comes from the underground waters [67]. The water fills the lake emerging through karst springs located on the north-western side and disappears through the swallow holes on the eastern side of the lake. The lake feeding is mainly due to the underground flows of the Isonzo river and the precipitations, recharging the karst groundwater [65,68,69]. The water level of the lake is extremely variable throughout the year, leading to a significant change in the extension of the lake surface area, from 200 m 2 during dry periods to 400,000 m 2 during wet periods [70].
The vegetation of the lake is characterized, starting from the center of the lake and going outwards, by aquatic communities (dominated by submerged and floating-leaved rooted plants, such as Potamogeton spp., Myriophyllum spp. and Nuphar lutea), by helophytic marsh vegetation dominated by tall sedges (Carex elata) and common reeds (Phragmites australis); behind these, a thin strip of mud soil along the lake shore hosts hygro-nitrophilous herbaceous communities. The zonation is completed toward the lake shores and landwards by hygrophilous willow shrubs and woods (with Salix cinerea, Salix alba), and a meso-hygrophilous woodland with the field elm (Ulmus minor), black poplar (Populus nigra) and narrowed-leaved ash (Fraxinus angustifolia subsp. oxycarpa), with its mantle with Ulmus minor and Paliurus spina-christi.

Sampling Design and Data Collection
The identification of the habitats occurring in the study area was based on the phytosociological analysis of plant communities, following the approach adopted in modern European habitat classifications [71][72][73][74]. According to a preliminary vegetation survey of the area, carried out in 2015, 33 plant communities were classified according to the phytosociological approach [75][76][77]. Due to the scarce extension of some communities (<1000 m 2 ) or to the small-scale heterogeneity in communities' spatial arrangements in some areas (where complex vegetation mosaics were present), only 20 vegetation types were effectively mapped ( Figure S1) and considered in this study as different habitat units (Table 1). IAS occurrences within the 20 vegetation types were sampled by means of a probabilistic sampling approach. Specifically, a stratified random approach was adopted using vegetation types as homogeneous strata. For each vegetation type identified, squared sampling units (hereafter plots) of 1 m 2 were randomly selected, with their number proportional to the total area occupied by the vegetation type. The number of plots was selected proportionally to the area of the communities as follows: 4 plots for all the communities with an extension up to 1.5 ha and 1 plot was then added every 0.7 ha up to a total extension of 5 ha; for communities bigger than 5 ha, the proportion 1 plot/ha was used. This sampling design was chosen to ensure the best compromise for representing small communities without oversampling the larger ones, thus maintaining a good compromise between statistical robustness and sampling effort. A total of 123 plots were sampled within the whole study area. We used such a small plot size (1 m 2 ) in order to reduce the margin effect due to the scarce extension and jagged shape of some vegetation types.
The plots were materialized in the field and only IAS abundance was recorded within each plot, measured as percent cover based on visual estimation. The data were collected during the late summer-autumn 2015. The nomenclature, taxonomy and classification of the invasion status of naturalized and invasive alien plant species follow Galasso et al. [78], and the syntaxonomic nomenclature of the plant communities up to the level of alliance follows Biondi et al. [79]. All alien species have been included in the null model analysis irrespective of their invasion status. Alien plants occurring in only 1 plot were excluded from further analyses.
To determine if the distribution, abundance and richness of IAS in the main habitats were significantly different from random expectation, we used a simulation approach based on Monte Carlo randomization [80]. The randomization is designed to produce a pattern that would be expected in the absence of a particular ecological mechanism [81]. To the best of our knowledge, this is the first attempt to use null models to test IAS distribution in relation to habitat invasibility, thus drastically reducing the amount of data required to test the null hypothesis, compared, for example, to correlative studies.
A similar approach was used by Bajocco and Ricotta [82] to identify land-cover types where fire incidence was higher (preferred) or lower (avoided) than expected from a random null model. This method was originally designed to study resource selection by animals [83,84]. Bajocco and Ricotta [82] considered fire as an "herbivore" with variable preferences for different resources (i.e., land-cover types); we considered IAS and habitats in the same way. The effect of the method was that the random frequencies were only dependent on the habitat extension, i.e., larger areas have higher probability to have a IAS.
Observed occurrences for each IAS in the sampled plots were randomly reassigned to the 7 habitats by the simulation, so that the probability of the species to be assigned to a given habitat was proportional to the area of that habitat itself. The null hypothesis was that IAS were randomly distributed in the study area while the alternative hypothesis was that IAS were clustered (or absent) according to the habitat type. Comparing the result of 9999 randomizations to the observed occurrence data, we evaluated the probability (pseudo p-value) that observed occurrences were significantly different from those of the simulations. For each habitat, p-values (two-tailed test, α = 0.05) were calculated as the proportion of Monte-Carlo-derived values that were as low or lower (or as high or higher) than the real ones. The same method was also applied for the randomization of IAS abundances.
On the basis of the IAS frequencies resulting from the null model simulations, we also calculated the IAS random richness for every habitat, reshuffling 9999 times the matrix assembled with the random frequencies of each IAS per habitat (keeping the IAS probability resulting from the null model fixed). In this case, the null hypothesis was that the observed IAS richness per habitat was not different from the one resulting from the simulated one.
Concerning life forms, there is a prevalence of perennial species (seven species), while annual species are 5. Species mainly originated from North America (10).
The most abundant IAS was B. frondosa, collected in 48.78% of the total sampled plots (60 out of 123), followed by B. vulgata and X. italicum, both collected in 16.26% of total sampled plots (20 out of 123). Overall, the average number of IAS per plot was 1.20 ± 1.29 (mean ± SD), while the average number of IAS per plot within each habitat was: 0.69 ± 0.48 for aquatic vegetation, 1.13 ± 1.06 for marsh vegetation, 4.29 ± 0.95 for lake shore herbaceous vegetation, 1.4 ± 1.07 for willow shrublands and woodlands, 0.78 ± 1.01 for meso-hygrophilous shrublands and forests, 0.67 ± 0.87 for nitrophilous forests and 2.2 ± 0.45 for plantations.

Null Model Simulations
The null models' results are shown in Tables 2 and 3. Table 2. For each habitat type, the observed frequencies (expressed as %) of the Invasive Alien Species (IAS) are shown with the 95% confidence intervals (C.I.) obtained from 9999 Monte Carlo simulations. In bold and with asterisks, the observed IAS frequencies are significantly different than expected (*** p < 0.001; ** p < 0.01; * p < 0.05).  Table 3. For each habitat type, the observed mean abundance (expressed as % cover of the plot) of the IAS are shown with the 95% confidence intervals (C.I.) obtained from 9999 Monte Carlo simulations. In bold and with asterisks, the IAS mean abundance is significantly different than expected (*** p < 0.001; ** p < 0.01; * p < 0.05). The number of occurrences (Table 2) was higher than expected from a random null model for E. nuttallii in aquatic vegetation, X. italicum in marsh vegetation, A. artemisiifolia, B. vulgata, C. campestris and X. italicum in lake shore herbaceous vegetation, B. vulgata in willow shrublands and woodlands, O. stricta in meso-hygrophilous shrublands and forests, R. pseudoacacia in nitrophilous forests and A. fruticosa in plantations, while the frequencies were lower than expected from a random null model for B. frondosa in aquatic vegetation, R. pseudoacacia in marsh vegetation and B. vulgata and X. italicum in meso-hygrophilous shrublands and forests.

Habitat Type
Where frequencies were significantly higher (or lower) than null models, the cover of IAS within the plots followed the same trend (Table 3), or, rather, IAS that had higher frequencies also had higher abundances and vice versa, except for B. frondosa in lake shore herbaceous vegetation, where its frequency was not significant but its abundance was higher than random, B. vulgata in willow shrublands and woodlands where its frequency was significantly higher but its abundance was not, B. frondosa and O. stricta in meso-hygrophilous shrublands and forests habitats, where the former had a significantly lower abundance and the latter had a significantly higher frequency. The species that did not differ by the random expectation, both in terms of frequencies and abundances, were P. quinquefolia and V. ×ruggerii. Overall, considering the results provided by Tables 2 and 3, there were 11 preferences and five avoidances for habitats by IAS.
The results of IAS richness per habitat analysis (Table 4) show that the observed IAS richness in aquatic vegetation, marsh vegetation and meso-hygrophilous shrublands and forests was significantly lower than random richness resulting from the null models.

Discussion
In this study, we used a null model simulation to test the association between alien species and different habitats in a wetland PA, taking into account only the habitat as a possible determinant for the observed resulting distribution. A first analysis of the alien floristic data shows that a third of the sampled IAS belongs to the Asteraceae family, confirming it as a successful invader taxon [85][86][87][88], among which, in particular, Bidens frondosa invaded indistinctly almost every habitat in the study area, except for the aquatic vegetation. Moreover, our results pointed out that IAS frequency, abundance and richness were not randomly distributed within habitats and that some types of habitat were more prone to invasion by IAS (see Tables 2-4).
The differences we observed in IAS distribution across habitats might be indirectly inferred by considering both habitat features and IAS functional traits and referring to IAS autoecology described in other studies [13,19,41,89]. Analyzing one habitat at a time and both the IAS frequencies and abundances that were significantly higher (or lower) than expected by chance, starting from the aquatic vegetation, we observed that the only successful invader here was Elodea nuttallii, a perennial submerged rooted hydrophyte native to North America. This species was detected for the first time in Lake Doberdò in 2011 [90] and seems to have completely replaced another IAS, Elodea canadensis, which was previously reported in the area [91]. E. nuttallii asexual propagation is reported to be improved by nutrient rich sediments [92] and its regeneration and colonization abilities are slightly higher than E. canadensis, but these characteristics alone seem insufficient to explain the displacement of E. canadensis by E. nuttallii in the area [93]. However, the replacement of E. canadensis by E. nuttallii has been observed in many European areas [94][95][96]. Conversely, B. frondosa occurrence and abundance in this habitat was significantly lower than random expectation, a feature easily explainable taking in account that the species cannot live submerged for long periods. In this case, it was pretty obvious that the only successful IAS would have been a plant adapted to live in the water, but, although trivial, this result served to verify the null model efficiency.
In the marsh vegetation, the observed occurrence and abundance of Xanthium italicum was higher than expected, while the occurrence and abundance of Robinia pseudoacacia was lower than random expectation. X. italicum was mostly found in those communities occupying the dryer parts of the marsh vegetation, particularly in the areas where nutrients accumulate [97]. For Lake Doberdò, these situations occur in the marginal areas of the lake, which dry out during the summer and are subject to the deposition of nutrients due to the direction of the flowing water. X. italicum fruits are easily dispersed by clinging to the fur of animals and its seeds can tolerate prolonged submersion without showing significant effects on final germination during the dry season [98], and, for these reasons, it is particularly competitive in this environment. Concerning the distribution of R. pseudoacacia in marsh vegetation, its lower observed occurrence and abundance is related to the species avoidance strategy to protracted submersion [99], a characteristic in common with most of trees with the exception for those adapted to wet conditions (e.g., Salix spp.).
The lake shore herbaceous vegetation is distinctly the most invaded habitat as shown both in terms of IAS significant occurrences and abundances (Tables 2 and 3). Indeed, this habitat type showed highly significant frequencies for four species (Ambrosia artemisiifolia, Bidens vulgata, Cuscuta campestris and X. italicum) and highly significant abundances for all the previous ones plus B. frondosa. Its high indiscriminate invasibility is probably related to three characteristics of the habitat, namely recurrent disturbs linked to the peculiar hydrological regime (water level oscillations, duration and seasonality of floods), the high availability of nutrients transported by water flow especially in spring period [100] and high solar energy availability. These characteristics are also confirmed by the observed poorly structured herbaceous native communities, which showed high dynamism in their structure during the year, shifting from the predominance of associations of the Agrostietea stoloniferae class in spring to that of associations of the Bidentetea tripartitae class in the late summer, which leads to an environment more susceptible to invasion [29]. We could detect the invasion success of A. artemisiifolia, an annual plant native to North America, as related to both high solar radiation (due to the available spaces in these plant communities) and the nutrients level in the habitat [101]. B. vulgata and B. frondosa were facilitated by the high dynamism in the resident plant communities, in which the native B. tripartita was replaced by both B. vulgata and B. frondosa. The latter is known to be more competitive than its native congener for many characteristics, such as a higher competitiveness where the nutrients level is high, a high plasticity of the vegetative phase length that allows it to grow taller than B. tripartita and to foster seed production and, finally, cause the greater resistance to drought during summer [102][103][104]. Similar characteristics are probably also present in B. vulgata, which is often reported in co-occurrence with B. frondosa [105,106]. Furthermore, Bidens seeds are easily dispersed by animals and humans for tens of km, explaining their wide distribution in the study area. For the last two species, X. italicum and C. campestris, which are strongly related because the latter grows often as parasite on the former, they have been found as strongly associated to the lake shore herbaceous vegetation, both IAS finding a suitable habitat because the first is a nitrophilous species, as stated before, and the second normally grows as a parasite by twisting around the stems of species of Polygonum, Xanthium and several other psammophilous plants [107].
As regards the willow shrublands and woodlands habitat, we observed an occurrence higher than expected for B. vulgata but its abundance was not significantly different than random. Its frequency suggests a certain degree of affinity for this habitat, due to nutrients accumulation, although lower than in the lake shore herbaceous vegetation. However, the not significant abundance could indicate the presence of limiting factors, for example, a lower light availability due to higher shading.
In the meso-hygrophilous shrublands and forests habitat, both the frequency and abundance of B. vulgata and X. italicum were lower than expected, while only the abundance of B. frondosa was significantly lower than random expectation. Conversely, Oxalis stricta occurrence was, here, significantly higher. O. stricta is a perennial herbaceous plant with overwintering buds at ground level, native to North America. It prefers moist soils, both in sunny and partially shaded areas [108]. Its occurrence here was favored for two reasons: the low or absent competition with the sparse herbaceous vegetation [30] and the high coverage of native shrubs that precludes potential invasions by other IAS. This is confirmed by the low frequencies and abundances of the other IAS. The habitat resistance to invasion is also favored by the slope that let water wash away nutrients from the soil after the lake floods.
R. pseudoacacia was highly associated with the nitrophilous forests habitat where it is known to add or replace the native Sambucus nigra [109], showing high cover values also due to asexual reproduction through root suckering. Here, other IAS were limited in their spreading (see Tables 2 and 3) due to shading conditions and to competition with highly abundant R. pseudoacacia.
For what concerned the poplar plantations, we observed the significant occurrence and abundance of Amorpha fruticosa, a woody shrub native to North America. A. fruticosa can tolerate dry soils, but it is most abundant along river or lake banks and at the edges of flooded forests. The plant grows well in medium to wet, well-drained soils in full sun to light shade and is tolerant of occasional flooding, lasting less than 30 days [110][111][112][113][114]. This explains its observed high frequency in the habitat, which is more favorable to the species because of the moist soil, the light shade and the occasional flooding events. Its high potential invasive behavior is due not only to its high tolerance of various environmental conditions but also to its reproduction strategies through self-seeding and suckers to form thickets [110]. A. fruticosa usually colonizes degraded wet habitats but also invades natural plant communities where it competes with native vegetation leading to a decrease in species diversity [115]. For this reason, it represents a serious threat for the diversity conservation, especially in wetlands.
Comparing observed and simulated IAS richness (Table 4) we can highlight that some habitats were less prone to a non-selective IAS invasion than others. While this may be quite obvious for aquatic vegetation, where only one of the sampled IAS was a hydrophyte, this is quite surprising for marsh vegetation and meso-hygrophilous shrublands and forests habitats. Although both habitats presented a relatively high number of sampled IAS, 7 and 5, respectively, they represented together the 77% of the study area, and, for the species richness-area relation [116], we could have expected to spot a higher IAS richness. These habitats were less prone to invasion, in case of the marsh vegetation because of extreme conditions caused by periodical and long-lasting floods that require special adaptions for plants to survive, and in that of meso-hygrophilous shrublands and forests because of a medium-low nutrients level and low light availability, which do not favor the spread of IAS, which are generally nutrient/light-demanding species.
In general, what we observed in the study area was that invasion was promoted in the habitats where one or more of the following conditions were satisfied: an intermediate disturbance level, caused by the water fluctuation or by the high dynamism of the communities within a habitat, a high level of nutrients and a medium to high light availability. These factors favored IAS invasion as they are known to be more competitive than natives both in a condition of intermediate disturbance [117] and when resources are easily available [30,41,89,118]. On the contrary, invasion by terrestrial IAS was avoided in the habitats where the communities presented a stable structure during the year, without seasonal changes, a medium nutrients level, moderate to low light availability or extreme conditions caused by long-lasting floods. Such different conditions in the small Lake Doberdò area were reflected as well by the native communities, represented by 33 different associations adapted to those peculiar environments. However, the high resource availability and the advantageous biotic interaction observed for some habitats may become a key driver along the subsequent stages of the invasion, leading to a loss in native diversity [19,41,89,119]. This is even more serious considering that the habitats are part of a wetland PA, in a global scenario of a dramatic decline in wetland areas due to extensive habitat degradation and fragmentation driven by human activities [56][57][58]. An early warning and eradication of IAS may impede new invasions and/or halt the ongoing ones: this would be particularly urgent in wetland PAs if we want to preserve their biodiversity, their conservation purpose and prevent the degeneration of native communities.
Regarding the possible weakness of the null model simulations, in some cases they were not capable of highlighting the actual distribution pattern of species; e.g., E. nuttallii is a hydrophyte and its distribution across the non-aquatic habitats was never significantly lower than random, as was expected. In particular, IAS avoidances may have been masked by confidence intervals in some cases (see Tables 2 and 3) as the confidence intervals have almost always the lower limit equal to 0. This probably happened because the random reassemble of the matrix was only weighted by the habitat area. Moreover, the null model suffered also for the low frequencies of certain species (e.g., P. quinquefolia and V. ×ruggerii) even though, overall, it provided useful results. In fact, we were not completely focused on highlighting the actual distribution of all IAS within the wetland, but rather on pointing out those IAS that were particularly threatening for the different habitats of the PA.
We also have to stress that the adopted sampling design was originally planned to evaluate the occurrence of alien species in different habitats in a rather complex area and to assess the conservation status of the native communities. A multi-scale approach should have been adopted considering that the pattern of variation in the alien/native species are strongly scale dependent [19]. The adopted plot size, then, somehow questionable for the characterization of plant community structure, can represent a limitation of this study and future analyses will take in account a multi-scale approach.
Nonetheless, the method tested in this study provides useful information to monitor and manage the IAS threat within PAs. By inferring important ecological property leading to habitat invasion without sampling the environmental characteristics of the habitats (an expensive and time-consuming operation in the field), we suggest using the applied method as a complementary tool in the analysis of biological invasions.

Conclusions
The results highlight an interesting linkage between IAS and invaded habitats, confirming that some habitats are more prone to biological invasion than others, and that environmental conditions can promote or avoid invasion depending also on the specific IAS ecology and biology. The study clearly shows the role of habitats in filtering the invasion depending on the species, besides the role of natural intermediate disturbance resulting on one side in an environment rich in biodiversity and on the other side advantaging the IAS colonization of the area [117].
Concerning the method, it enabled us to derive useful information, proving to be a low-cost analysis to routinely integrate PAs' management activities. In particular, it may help to detect highly endangered habitats and to invest time and money more efficiently to control IAS, and, when possible, also to prevent IAS colonization.

IAS
Invasive Alien Species PA Protected Area