Distribution of Five Aquatic Plants Native to South America and Invasive Elsewhere under Current Climate

: Biological invasions and climate pose two of the most important challenges facing global biodiversity. Certainly, climate change may intensify the impacts of invasion by allowing invasive plants to increase in abundance and further expand their ranges. For example, most aquatic alien plants in temperate climate are of tropical and subtropical origins and the northern limits of their ranges are generally determined by minimum winter temperatures, and they will probably expand their distributions northwards if climate warms. The distribution of ﬁve invasive aquatic plants in freshwater systems across continents were investigated. Their global distributions in the current climate were modeled using a recently developed ensemble species distribution model approach, speciﬁcally designed to account for dispersal constraints on the distributions of range-expanding species. It was found that the species appear capable of substantial range expansion, and that low winter temperature is the strongest factor limiting their invasion. These ﬁndings can be used to identify areas at risk of recently introduction of neophytes, and develop future monitoring programs for aquatic ecosystems, prioritizing control efforts, which enables the effective use of ecological niche models to forecast aquatic invasion in other geographic regions.


Introduction
Freshwater ecosystems are vulnerable to both alterations in climate and biological invasion [1]. In future years, water temperature is expected to increase, and the modification of precipitation events may alter flow regimes [2,3]. As modifications in climate alter aquatic ecosystems by increasing nutrient loading from the catchment areas to the lakes [4,5], influencing streamflow, amount and duration of ice cover, water quality and chemistry, then aquatic invasive species are predicted to invade previously uninhabitable areas [6,7], and even aquarium plants that are not currently invasive would likely become invasive due to the fact that most plants traded in aquaria are of tropical and subtropical origins [8][9][10]. Consequently, processes, such as flowering, nutrient uptake, and competition between species and the metabolic events, can be altered [11,12]. In addition, compared to terrestrial plants, aquatic plants are shown to have a higher probability of becoming invasive in new environments [13]. For example, Egeria densa Planch. (Hydrocharitaceae); Eichhornia crassipes (Mart.) Solms (Pontederiaceae) (Pontederia crassipes Mart.); Myriophyllum aquaticum (Vell.) Verdc. (Haloragaceae); Pistia stratiotes L. (Araceae), and Salvinia x molesta D.S. Mitch. (Salviniaceae) are five aquatic plant species native to South America and highly invasive elsewhere [14][15][16][17][18]. They have been introduced into Europe, North America, Australia, and New Zealand, where they have been increasingly expanding their ranges [19][20][21][22], and have been reported as the most damaging alien plants in South Africa's freshwater systems [10].
In both terrestrial and freshwater ecosystems, species distribution models (SDMs) can be effectively used to generate risk maps that predict the relative suitability of habitat, if the introduction occurs outside the native range (e.g., [23][24][25][26][27][28][29][30]). Introduction events across continents (primary pathway) or spread after a single introduction (secondary pathway) [31] may promote the capacity of the species to respond to selection in new environments [32,33]. SDMs have been used to predict potential ranges under current and future climate conditions [34,35], for native and invasive species [36][37][38][39]. SDMs have generally been proven to reflect the correct response to climate change and other environmental gradients [40]. Species distributions are affected by many abiotic and biotic interactions [30,41,42], including anthropogenic factors [38]. However, climate is often regarded as a critical determinant of plant species ranges and ranges are correlated more strongly with climate than with any other factor especially on broad spatial scales [41,43]. This approach may only be appropriate if assumption of equilibrium is not violated (i.e., potential dispersal constraints and environmental tolerances are often ignored), and assuming there are no missing values in the climate space (i.e., not truncated niche) [44]. Distribution models identify sets of variables associated with the presence of invasive species to project their requirements onto the geographic space [45]. Matching the climate of the native range to non-native areas has provided very valuable information as to potential risks of establishment and spread. However, such an approach assumes that the same interactions between biotic and climatic factors that limit range size in the native range operate in the introduced range [46].
One assumption underlying SDMs is the principle of niche conservatism. This concept assumes that species tend to preserve their climatic and ecological requirements over time and space [47] or only change slowly over hundreds to millions of years [32,48]. However, during biological invasions, alien species could expand beyond their predicted climatic envelope into the invaded range (i.e., climatic niche shift), occupying new niches that are not predictable from knowledge of the native range alone. The species then change environmental tolerances to deal with the climate change [37]. Assuming niche conservatism, then invasive ranges can be predicted with models fitted with data from the native range [49]. However, the reliability of predictions provided by SDMs is questionable if niche shifts have truly taken place in different areas [47,48,50]. Many studies have investigated invasive plant species distribution through projection by SDMs [51][52][53] but only a few have considered non-native aquatic plant species [10,[54][55][56]. Indeed, few aquatic plants have been studied with SDMs at the continental scale, for example Hydrilla verticillata [57], Elodea canadensis [54], E. densa, Ludwigia sp.pl., and M. aquaticum [58]. Nevertheless, some studies considered the distribution of some of these species and other invasive aquatic species (e.g., Ludwigia grandiflora Greuter and Burdet, and Myriophyllum heterophyllum Michaux) at a country scale or based on a limited set of algorithms or bioclimatic variables useful mostly for terrestrial plants [59,60].
Invasive plants have the potential for major further global range expansion. Therefore, this study aims to disentangle the role of hydroclimatic predictors on the occurrence of five South American invasive aquatic plants across continents, and determine which predictors influenced projected habitat suitability and potentially limit range expansion into unoccupied regions. Thus, the following questions were asked: (a) what is the potential for further invasive range expansion of these five species in the current climate? (b) What are the key climatic constraints shaping their occurrences? The predictions of high-risk areas for invasion provide valuable information for prevention, better surveillance, and rapid responses [61], as well as informing about the potential effects of climate on invasion risk.

Study Species and Occurrence Data
Five well-known hydrophytes native to South America and invasive elsewhere were selected, namely Egeria densa, Eichhornia crassipes, Myriophyllum aquaticum, Pistia stratiotes, and Salvinia molesta [62][63][64] for distribution modeling. Among them, two are categorized as submerged (E. densa and M. aquaticum) and three as free-floating (E. crassipes, P. stratiotes, and S. molesta). These five-species were selected based on the high level of invasiveness reported by [65] in the risk assessment for aquatic species in South America. In addition to the native distribution, we distinguish locations where the species has become invasive within South America, but outside the historic native range (classified as alien in) were distinguished. Importantly, three of them (E. crassipes, M. aquaticum, and S. molesta) have been included in the list on invasive alien species of (European) Union concern (Regulation (EU) No. 1143/2014 of the European Parliament and of the Council of 22 October 2014, on the prevention and management of the introduction and spread of invasive alien species, hereafter, IAS Regulation). In addition, four of them are regulated by current Italian legislation (E. densa, E. crassipes, M. aquaticum, and S. molesta), identifying the most relevant regulated species. The Regulation establishes aspects as early warning and surveillance procedures, rapid eradications to avoid establishment and control systems [66].
Dubious records, such as those where the species has been introduced outside its climatic niche for attempts of establishment (e.g., those where winter frosts are known to kill the species or populations known to occupy climatically anomalous microhabitats) or cultivated were removed. For example, repeated introductions, of aquatic plant species, failed to establish in Germany up until 2005. However, since 2008, an established population of P. stratiotes has been permanently present in thermal sections of the River Erft [72]. In particular, S. molesta has been found in Austria, Belgium, France (Corsica), Germany, Italy, the Netherlands and Portugal, but it is not clear if reports represent established populations [73]. Botanic garden records were also excluded as well as occasional field observations (e.g., ephemeral presence of E. crassipes and P. stratiotes in parts of Germany, Russia, the Netherlands, and the Czech Republic [72,74]). Literature documenting failed establishments (e.g., non-permanent populations occurring in cold parts of USA that proliferate for summer and then die and survive or establish during the warmer months of the year, and may have an extremely limited (or no capacity) to overwinter under natural conditions [14,18]) was also considered [15,[17][18][19]63,75].
After data cleaning, a total number of geographic records equal to 652 for E. densa, 3210 for E. crassipes, 1485 for M. aquaticum, 1688 for P. stratiotes, and 704 for S. molesta was obtained ( Figure 1). Then, the records were gridded at a 0.25 × 0.25 degree resolution for global modeling (Table S1).

Climate and Freshwater System Data
Freshwater-specific climate data were collected from near-global environmental information for freshwater ecosystems (https://www.earthenv.org/streams; [76]). The dataset consists of near-global, spatially continuous, and freshwater-specific environmental variables in a standardized 1 km grid. The same author delineated the sub-catchment for each grid cell along the HydroSHEDS (Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales, https://www.hydrosheds.org [77]), river network and summarized the upstream environment (climate, topography, land cover, surface geology, and soil) to each grid cell using various metrics (average, minimum, maximum, range, sum, inverse distance-weighted average, and sum). The hydroclimate layers were found to correlate more strongly with water temperature and flow rates than the equivalent bioclimatic variables (WorldClim database [78]), typically used in distribution modeling.
Then, it was considered useful and novel for this specific aquatic research the use of the 19 long-term variables, there after called "hydroclimatic variables". All variables were subsequently averaged across single lakes and reservoirs of the Global Lakes and Wetlands Database [79] that are connected to the river network. Figure 1. Distribution records of the five aquatic species that were modeled, plotted as white dots located inside the native range (red polygon) and in the invasive range for Egeria densa, Eichhornia crassipes, Myriophyllum aquaticum, Pistia stratiotes, and Salvinia molesta. Regions highlighted in orange are considered highly unsuitable for each of the five-investigated species and, therefore, included in the model background (pseudo-absences). Global maps exclude parts of the world that may be environmentally suitable, but are out of dispersal range (grey).
Based on natural history data of each aquatic species, the following hydroclimatic and environmental variables were selected for the modeling: hydro 06 = minimum upstream temperature of coldest month (reflecting exposure to frost); hydro 10 = mean upstream temperature of warmest quarter (reflecting the growing season thermal regime) and hydro 18 (log10 (1 + hydro_18 transformed mm) = upstream precipitation of warmest quarter (seasonal drying out of waterbodies may reduce suitability). Then, the variance inflation factor (VIF, "usdm" R package [80,81]) was applied in order to exclude multi-collinearity between variables. Since the model have aquatic habitat requirements, projected suitability maps were filtered to only include cells containing aquatic ecosystems using the Global Lakes and Wetlands Database (GLWD) [79].

Background Domain
To generate potential distributions of the species in their native and invaded ranges, presence-background modeling (presence-only) were used, drawing the background points (pseudo-absences), using a biological-informed pseudo-absence selection approach [82]. This allows algorithms to contrast the environment at the species occurrence locations (e.g., data selected from suitable areas) against a sample of background environmental conditions (i.e., pseudo-absences). The background points are sampled from both an accessible region [83] and an unsuitable area for the species based on its environmental tolerances. As such, both available environmental data and prior biological knowledge were used to fit the model (Table 1).
Specifically, the background region comprised: (1) the accessible native range (i.e., minimum convex polygon created around native occurrences [84]); (2) accessible parts of the invaded range (i.e., the area that the species has potentially visited, called "M"). In other words, the area within which presences may exist and within which absences are meaningful, in that they represent sites with the broader background likely to have been tested by the species for suitability, but not occupied yet. In the Biotic, Abiotic, and Movement (BAM) diagram, A represents the geographic region with the appropriate set of abiotic factors for the species and may be regarded as the geographic expression of the fundamental niche; B is the region where the right combination of interacting species occurs and M are those parts of the world "accessible" to the species in some ecological sense, without barriers to movement and colonization [85]; and (3) highly unsuitable regions (i.e., areas outside "M", where the species cannot occur). Table 1 shows the hydroclimatic variables used to define the unsuitable background region, based on the climatic conditions and environmental limits for the five species documented in the literature. According to the literature, it was assumed that the key environmental constraints on plant distributions are mostly determined by: mean minimum temperature of the coldest month ranges between <−3 and 0 • C (tolerance freezing temperatures and prolonged exposure to lethal frosts); mean temperature of the warmest quarter ranges between <10 and 11 • C (below the minimum growth temperature) and precipitation of the warmest quarter < 4-5 mm (only small and seasonally dry habitat is available, which is expected to be of low suitability) [62,63,72,86,87].
Areas outside "M" will be predicted at lower suitability levels and inclusion of these areas in the model makes the model look better than it actually is. Therefore, accessible area is used to attempt to account for dispersal constraints and the background sampling protocol is designed to exclude "suitable-but-not-reached" areas [83]. According to [82], dispersal limitations in unsuitable areas could underestimate occupancy in those regions. Nevertheless, this should not impede the prediction of the region suitable for invasion.
To control for spatial biases in recording effort, the target group approach was used [88]. For this, the background sampling was weighted in proportion to the tracheophyte record density on the Global Biodiversity Information Facility -GBIF-(i.e., weighted selection of background data, with higher probability sites being selected most often), which was assumed as a good proxy for recording the effort of the focal species. To sample as much of the background environment as possible, twenty background samples of 10,000 randomly generated points (i.e., pseudo-absences) were selected, which is sufficiently large to ensure model convergence ( Figure S1).

Species Distribution Models SDMs
Species distribution models to predict potential suitable areas for the invasion of the five aquatic species were used, employing BIOMOD (biomod2 R package [89]). Each dataset (i.e., background samples from only the accessible area and background samples from both the accessible and unsuitable area) was randomly split into 80% for model training and Variable importance was obtained, and variable response functions were produced using BIOMOD. Model predictive performance was assessed by calculating the Area Under the Receiver-Operator Curve (AUC) and True Skill Statistic (TSS) for model predictions on the evaluation data, which were reserved from model fitting. AUC can be interpreted as the probability that a randomly selected presence has a higher model-predicted suitability than a randomly selected background location [90]. This information was used to combine the predictions of the different algorithms to produce ensemble projections of the model. For this, the two algorithms with the lowest AUC were first rejected and then predictions of the remaining algorithms were averaged, weighted by their AUC (Table S2).  [60].
<4 mm; germination occurs when substrates are exposed as water recedes and as dry substrates are moistened when water levels rise. Seeds survive in wet mud and are long-lived [63].
<0 • C; below the minimum growth temperature. Consistent with the temperature inhibiting seed germination [75].
Salvinia molesta <0 • C; prolonged exposure to lethal frosts [86]. Buds are killed when exposed for more than two hours to temperatures below −3 • C.
<10 • C; below the minimum growth temperature.

Current Potential Distribution
Acceptable numbers of grid cells with presences were obtained for modeling the five invasive aquatic species (645 grid cells for E. densa, 1077 for E. crassipes, 637 for M. aquaticum, 817 for P. stratiotes, and 237 for S. molesta). Figure 1 shows their current distribution and the regions considered highly unsuitable for these species based on the rules in Table 1. The areas estimated to be climatically suitable for the invasive aquatic species under current climatic conditions were plotted at global scale ( Figure 2). The projection of E. densa, E. crassipes, and S. molesta models indicated high suitability throughout the tropical and subtropical regions of the world (areas in yellow, Figure 2). In addition, areas suitable for the distribution of M. aquaticum and P. stratiotes were found in subtropical areas. Non-native occurrences of the species were largely consistent with these projections. Among the five species, the potential geographical range for E. crassipes and S. molesta was predicted to be extremely broad. In general terms, their region is restricted by cold stress (Figure 3). No climatically suitable areas for the species were identified in the Andean cordillera in Chile (some high-altitude areas are too cold for persistence), but a high suitability in southeast North America, possible causing an invasion threat where there are sources of standing water. In Europe, models predict that E. densa and M. aquaticum may be able to invade widely in northern regions and relatively frost-free parts of Central Europe. In Australia, the suitable areas were located in the north (E. crassipes, P. stratiotes, and S. molesta) and east (M. aquaticum).

Model Evaluation
Across the models, the discrimination ability (AUC and TSS) between aquatic species records and grid cells without records ("pseudo-absences") is shown in Table S2. The minimum and maximum AUC values from model outputs generated by ensemble model for the species ranged from 0.867 to 0.970. These results indicate a relatively high performance of the species distribution models. However, the inclusion of the unsuitable region inflates AUC values, as those places have no presences, and it is easy for the models to predict low suitability there. The variable importance of the fitted model algorithms and the ensemble suggested that suitability was most strongly determined by the minimum upstream temperature of coldest month ranged between 53.3% and 85% among the five species ( Figure 3, Table S2). The response curve plots show that the ensemble model estimated biologically reasonable curves, with suitability limited by harsh frosts, low growing season temperatures particularly in the case of E. crassipes, low cover of large wetlands, and low precipitation in the growing season ( Figure 3). The function also indicated that suitability was reduced if minimum temperatures were too high.

Distribution in the Native Range
In this study, suitable areas for invasion under the current climate for five invasive aquatic plants were predicted, suggesting that extensive regions are predicted to be suitable for these species. They are all classified as invasive, which means they can have massive environmental and economic negative impacts over large areas.
In South America, where the species are native, but also introduced outside their native ranges, the highest suitability was obtained across Amazonia (e.g., E. densa, E. crassipes, P. stratiotes, and S. molesta). The climatic potential range appears to extend into cooler areas of Argentina than they have been reported (e.g., E. densa, E. crassipes, M. aquaticum, and S. molesta), in the lower slopes of the Andean cordillera in Colombia, and in evergreen forests in southern Venezuela-except for M. aquaticum-all characterized by dry and warm climates.
The proportion of suitable environments for aquatic species is predicted to decrease in the future in their native continent (e.g., E. densa and M. aquaticum) [58], and the potential bioclimatic range is conditional upon multiple factors that can impede the establishment in new environments [45].
Importantly, these species have their native range in South America and are recorded as invasive aliens in many regions inside their native continent, i.e., they act as alien in South America [65,68]. Distribution in the native continent may not be at equilibrium with the environment because the capacity of alien species to spread and occupy all climatically suitable habitats is limited by their dispersal ability constraints and environmental tolerances [44,91,92], which is often a feature of invasive aquatic species and is strongly associated with dispersal barriers at the boundaries of watersheds [93,94]. A strength of the modeling approach is that unoccupied but potentially climatically suitable parts of the native and non-native continents were excluded from the pseudo-absence sampling; thus, these regions were not erroneously modeled as unsuitable.

Distribution in the Non-Native Range
All of the major clusters of non-native records (e.g., Australia, Europe, southern USA, and South Africa) were projected to be suitable. Southern USA and the Mediterranean Basin were predicted as suitable regions for all the five species to invade. Northern Australia and India were projected to be suitable for E. crassipes, P. stratiotes, and S. molesta. The potential impacts of climate change in South Africa reveal, an according to [10], which may favor the expansion of E. crassipes and S. molesta. In North America, it is assumed that the presence of E. crassipes is due to repeated re-invasion from anthropogenic release by humans, rather than propagation by seed [63]. According to [95], bioclimatic models suggested an increase in emergent aquatic plant species followed by an expansion in their non-native distribution. The results therefore highlight that further range expansion appears likely for all species and there is a risk of colonization of new areas, but it is also an opportunity to manage the invasion process.
The northern limit of the potential distribution for E. crassipes in Europe is shaped by cold stress, according to the model, with climatically suitable habitat encircling the Pyrenees in southern France. All of the countries of the Mediterranean Basin and North Africa are at significant risk, in terms of degree of climatic suitability. The possible reason for the quiterestricted invasion of E. crassipes in Mediterranean regions (Italy, Portugal, Spain, Turkey, and Corsica, e.g., [16,[96][97][98]) is very likely due to removal and control interventions taking place, and the fact that regrowth can occur mostly in spring and summer, with warmer climate conditions. Although prolonged cold weather may kill many plants, portions of the mat or the seeds may remain viable and allow regeneration when favorable conditions return. However, according to future climate change projections, a great potential for future range expansion lies in Europe [63]. Nonetheless, the standing water resources throughout the Mediterranean are critical to human survival and agriculture production due to the summer drought period that typifies the Mediterranean climate [15]. On the contrary, the climatic potential non-native range for M. aquaticum and E. densa appears to extend into cooler areas of north Europe. Importantly three of them (E. crassipes, M. aquaticum, and S. molesta) are included in the list of invasive alien species of (European) Union concern (Reg. no. 1143/2014).
In general, for all five species, the potential non-native range extends mostly towards west Europe, where summers are cool, but winters are mild (i.e., low probability to find frost and more likelihood to survive). The predicted range of M. aquaticum shows that this species has suitable climatic areas at latitudes, which are further north than those suitable to the others aquatic plants. This result is consistent with models assessed by [60] in Ireland and [58] at global scale, where their distribution could expand. In accordance with [58], M. aquaticum generally has the lowest habitat suitability throughout the world. The same authors reported that for E. densa and M. aquaticum, the proportion of suitable environments is predicted to decrease in the future in their native range, as well as in other areas of the world, where they have not yet been reported. However, the same authors reported-and this is in accordance with the results obtained in this study -that E. densa and M. aquaticum have a higher cold tolerance that may allow to thrive in Nordic countries.

Climatic Limits to Occurrence
The predictions on the current climate ranges, of the five invasive aquatic species investigated in this study, reflect their current distribution and climatic associations. This suggests there are suitable (but not yet occupied) regions where the establishment is possible under the climatic and environmental conditions, but where invasion has not been reached through dispersal. According to [58], one reason could be the time, since introduction into their invasive range has not always allowed them to be in equilibrium with the environment, and likely has not realized all of their suitable environments yet, which is also seemingly the case in the native range, where they are restricted to particular watersheds. Since in many cases the species are not at equilibrium, model performance is not accurate during the invasion process, and false-absences may be misleading if they represent areas that have yet to be colonized [99] and are treated as pseudo-absences in model building. If suitable, unoccupied regions out of dispersal range are included in the background area, then the model will confuse lack of presence of the species due to dispersal limitations with a lack of presence due to environmental unsuitability, which may bias the fitting of species-environment and the prediction of potential distributions in novel regions. Nevertheless, if the accessible area is restricted by polygons around concentrated regions of occurrences (e.g., invasive range)-as in the biologically-informed approach, to reduce this bias within dispersal range of the occurrence [82,83]-the model can result in better model performance, more than if the models were built using a geographic background much larger than the core area of the species distribution [84].
The models developed here did not consider local abiotic factors, such as water body size, water depth, or water quality (pH, nutrients, salinity, and turbidity), or biotic interactions, which can influence the suitability of aquatic environments for the aquatic plants [74,[100][101][102][103][104], or the rate of spread. In fact, these factors can be analyzed and modeled only at a finer scale. However, in general, climate is considered the driver most strongly associated with species distributions at global scale, but may be the least important factor associated with the regional or local distribution of aquatic invasive species [60]. In particular, water availability depends on precipitation regimes [1,55,105,106], while temperatures affect growth, survival, and reproduction [95]. According to [60], it is possible that a stronger association with climatic variables would emerge if climate variables more specific to freshwater habitats were available and used (e.g., maximum, minimum, and mean water surface temperatures). Unlike other distribution modeling studies for aquatic species, the models applied in this research have used the available hydroclimatic variables, which better reflect temperature and flow regimes in aquatic environments [76]. Nevertheless, wherever the water temperature is artificially warm enough (e.g., warming associated with industrial outflows, thermal waters), these species might establish. In addition, nutrients (eutrophication) and water drawdown can have positive effects on the overwintering of alien aquatic plant, as in the case of E. crassipes, at the northern margin of its distribution in China [107], or can promote their establishment and spread [108].

Implications for Climate Change
Invasion ecologists have attempted to predict possible effects of climate change on plant invasions and range expansions by projecting current species temperature limits onto maps of future temperature conditions [109]. In theory, the hydroclimatic variables used in this study can be produced for climate change scenarios, allowing a mapping of potential future distributions, but these are not currently available [76]. Nevertheless, examination of the species response curves suggests potential responses of the species to climate change. For example, the Mediterranean Basin is predicted to shift towards drier and warmer conditions [110]. The model response curves suggest that this warming is likely to increase suitability for all species. Possible mechanisms for this include increased survival through warmer winters, and longer growing seasons, permitting increased reproductive output. However, on the other side, lower precipitation may reduce suitability for the species in this region as a result of drying out of waterbodies or competition for water uses, such as crop irrigation [111]. In northern Europe, cold winters act as a survival barrier that prevents warm-water-adapted species from establishing self-sustaining populations [7]. Indeed, the distribution of E. crassipes and P. stratiotes is currently thought to be limited by cold, hard freezes, or ice cover [86]. Given the response curves of the five modeled species, aquatic invasive species may be better able to overwinter if temperatures increase, snowfall is reduced or seasonal ice cover melts earlier in the year [112]. In these regions, better overwinter survival may reduce the effectiveness of current management by mechanical means or hand pulling, as has been seen for Eichhornia species in North America [44]. Shorter winters at high latitudes will allow longer growing seasons, increasing the depths of macrophyte establishment [113]. Indeed, the occurrence of several warm winters in the late 1940s was responsible for the invasion of some non-native species (i.e., the spread of some introduced species appears to be temperature dependent) [114].
The niche of aquatic plants is not always conserved among different regions [115]. Therefore, habitat suitability in non-native range could be extended because of the current climate, resulting in the expansion of the distributions of all species, when they could disperse, but it is possible that the niche expansion could allow colonization of regions projected to be unsuitable. Studying niche shifts during current invasion could inform about the potential to undergo further niche shifts. The potential lack of niche conservatism in aquatic plants could be due to: (1) an unsuitable or incomplete set of explanatory variables; (2) genetic or epigenetic variation; or (3) phenotypic plasticity within species. According to [116], genetic variation may explain niche shift between continents, whereas phenotypic plasticity may be a more important factor in within-continent comparisons.

Conclusions
This study employs global-scale modeling approaches and predictions specifically developed for invasive non-native aquatic species distributions, using previous knowledge of the species tolerances and dispersal constraints. It is also one of the first studies to develop species distribution models using climatic predictors specifically produced for freshwater ecosystems. The models predicted that five highly invasive aquatic plants, natives to South America, have the potential for major further global range expansion in the current climate. Although the areas predicted as suitable for a species may not necessarily mean that the populations will establish there, it does offer useful information for detecting areas of potential invasion and spread [100] and, therefore, for prevention and early detection. A strong barrier to further range expansion is represented by low winter temperature, low summer temperature, and low precipitation, suggesting that warming of temperate regions may allow further poleward range expansion, while altered precipitation patterns in more arid regions could also affect suitability. As such, the distribution models provided in this study are useful to support risk assessment and management of these aquatic invasive plant species as they help in identifying regions that are currently at risk of invasion. It is suggested that a global-scale approach is useful for characterizing the species climate associations, but that further fine-tuning assessment of niche differences across the global range and assessment of the role of habitat conditions in determining local distribution patterns would also be extremely beneficial.
Supplementary Materials: The following are available online at https://www.mdpi.com/2673-4 133/2/1/3/s1, Figure S1: Selected background grid cells used in the modeling of Egeria densa, Eichhornia crassipes, Myriophyllum aquaticum, Pistia stratiotes, and Salvinia molesta. Table S1: Number of occurrences before (raw data) and after cleaning procedure and numbers of grid cells with presences used for modeling. Table S2: Summary of the cross-validation predictive performance (AUC and TSS) and variable importance of the fitted model algorithms and the ensemble (AUC-weighted average of the best performing seven algorithms) for five South American occurring elsewhere as invasive non-native species. Variable importances are given as percentages for Hydro6 (mean minimum temperature of coldest month), Hydro10 (mean temperature of warmest quarter) and Hydro18 (precipitation of warmest quarter) and Global Lakes and Wetlands.

Funding:
The study was funded by the COST Action TD1209 "European Information System for Alien Species" (Alien Challenge).

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.