AJBMAnales del Jardín Botánico de Madrid0211-13221988-3196Consejo Superior de Investigaciones CientificasAJBM201709-245210.3989/ajbm.2452ArticlesAn invasion risk map for non-native aquatic macrophytes of the Iberian PeninsulaMapa de riesgo de invasión de macrófitos acuáticos exóticos de la Península IbéricaInvasion risk map for non-native macrophytes of the Iberian PeninsulaRodríguez-MerinoArgantonio1*Fernández-ZamudioRocío2García-MurilloPablo1
Department of Plant Biology and Ecology, Faculty of Pharmacy, University of Seville, Profesor García González St., no. 2, 41012 Seville, Spain; argantonio.rodriguez@gmail.com
Doñana Biological Station, CSIC, Américo Vespucio Ave., s.n., 41092 Seville, Spain
Corresponding author
ORCID ID: A. Rodríguez-Merino (http://orcid.org/0000-0002-1568-5087); R. Fernández-Zamudio (http://orcid.org/0000-0001-5804-9518); P. García-Murillo (http://orcid.org/0000-0002-1761-9569).
Freshwater systems are particularly susceptible to non-native organisms, owing to their high sensitivity to the impacts that are caused by these organisms. Species distribution models, which are based on both environmental and socio-economic variables, facilitate the identification of the most vulnerable areas for the spread of non-native species. We used MaxEnt to predict the potential distribution of 20 non-native aquatic macrophytes in the Iberian Peninsula. Some selected variables, such as the temperature seasonality and the precipitation in the driest quarter, highlight the importance of the climate on their distribution. Notably, the human influence in the territory appears as a key variable in the distribution of studied species. The model discriminated between favorable and unfavorable areas with high accuracy. We used the model to build an invasion risk map of aquatic macrophytes for the Iberian Peninsula that included results from 20 individual models. It showed that the most vulnerable areas are located near to the sea, the major rivers basins, and the high population density areas. These facts suggest the importance of the human impact on the colonization and distribution of non-native aquatic macrophytes in the Iberian Peninsula, and more precisely agricultural development during the Green Revolution at the end of the 70’s. Our work also emphasizes the utility of species distribution models for the prevention and management of biological invasions.
Resumen
Los sistemas acuáticos son especialmente susceptibles a los organismos exóticos debido a su elevada fragilidad y a los impactos que provocan estas especies en este tipo de hábitats. Los modelos de distribución de especies, basados en variables ambientales y socioeconómicas, facilitan la identificación de las áreas más vulnerables ante la expansión de especies exóticas. Se utilizó MaxEnt para predecir la distribución potencial de 20 macrofitos exóticos en la Península Ibérica. Algunas de las variables estudiadas, como la estacionalidad de la temperatura y la precipitación del cuatrimestre más seco, ponen en evidencia la importancia de los factores climáticos en su distribución. Además, la influencia humana en el territorio se presenta como una variable clave en la distribución de las especies estudiadas. El modelo obtenido discrimina claramente entre áreas favorables y desfavorables con mucha precisión. Se utilizó el modelo para construir un mapa de riesgo de invasión de macrófitos acuáticos para la Península Ibérica que incluyó los resultados de 20 modelos individuales y que muestra que las áreas más vulnerables son las zonas cercanas al mar, las cuencas de los grandes ríos y las zonas con una alta densidad de población. Estos resultados vinculan la importancia del impacto humano en la colonización y la distribución de los macrófitos acuáticos exóticos en la Península Ibérica y, más concretamente, con la Revolución Verde de finales de la década de los setenta. Nuestro trabajo enfatiza la utilidad de los modelos de distribución de especies para la prevención y gestión de invasiones biológicas.
KeywordsAquatic plantsbioclimatic factorsbiological invasionsecological niche modelsfreshwater ecosystemsmap risk assessmentMaxEntnon-native speciessocio-economic factorsspecies distribution modelPalabras claveEcosistemas acuáticos continentalesespecies exóticasfactores bioclimáticosfactores socioeconómicosinvasiones biológicasmapa de evaluación de riesgosMaxEntmodelos de nicho ecológicomodelos de distribución de especiesplantas acuáticasINTRODUCTION
Invasive species are one of the main causes of biodiversity loss. At a global scale, they represent a major threat to the ecosystems functioning (Mack & al., 2000; Sala & al., 2000; Brooks & al., 2004). Non-native species may also cause negative effects (Ricciardi & Kipp, 2008; Pyšek & Richardson, 2010) on human health (Hulme, 2006; Chytrý & al., 2009), as well as important economic impacts (Pimentel & al., 2005). Some freshwater systems are considered biodiversity hotspots (Murphy, 2002; Strayer & Dudgeon, 2010; Brundu, 2015; Serrano & Díaz Paniagua, 2015) and are one of the most threatened ecosystems in the world (Collen & al., 2014; Brundu, 2015; Serrano & Díaz Paniagua, 2015). These systems are particularly susceptible to biological invasions, because of their propensity to shift away from natural conditions and feedbacks that alter colonized habitats (Willby, 2007; Aguiar & Ferreira, 2013; Brundu, 2015; Gallardo & al., 2015). Aquatic macrophytes play an important role in the structure and function of freshwater systems (Chambers & al., 2008; García-Murillo & Fernández-Zamudio, 2015) by providing a structurally complex environment (Rennie & Jackson, 2005; Dibble & al., 2006). They contribute to environmental heterogeneity (Harrel & Dibble, 2001) and to increase the diversity of ecological niches. Aquatic invaders features like high productivity, broad ecological tolerances, notable phenotypic plasticity, and a remarkable facility in producing propagules (Santamaría, 2002; Les & al., 2003), have led some invasive freshwater plants to belong to the group of the “100 of the World´s Worst Invasive Alien Species” (Lowe & al., 2004). In addition, the nutrient increase in many water bodies due to human activities and the frequent absence of natural enemies in this group of plants, have led in some cases to absolute dominance in the invaded habitats (García-Murillo & al., 2007; Ruiz & al., 2008).
An early detection of the arrival of non-native species can increase the success in their eradication before the establishment, preventing future invasions (Broennimann & Guisan, 2008; Williams & Grosholz, 2008; Crafton, 2015). For this reason, it should be necessary to identify the most exposed areas of invasion risk (Reshetnikov & Ficetola, 2011). But aquatic habitats in general and aquatic macrophytes in particular are difficult to be monitored (Brundu, 2015). So, the development and use of alternative methodologies for the prevention and control of exotic species are essential for the identification of areas with a high invasion risk. This kind of methodologies will allow us to manage potential non-native species while preserving native species (Gallardo & al., 2012).
Species distribution models have the potential to predict invasiveness and have become common in the study and management of biological invasions (Peterson, 2003; Thuiller & al., 2005). Significant recent advances have been achieved in the development of species distribution models (v.gr., Elith & Leathwick, 2009). Appropriate factors in modeling the potential distribution of species, as well as the use of suitable occurrence data, are essential to execute more accurate models. In our case, we have chosen the algorithm MaxEnt (Phillips & al., 2006), based on the maximum entropy principle, for modeling the potential distribution of non-native aquatic macrophytes. Several authors propose that MaxEnt model is better than other algorithms based on presence-only data (Elith & al., 2006; Elith & Leathwick, 2009; Mateo & al., 2010).
The Iberian Peninsula has been considered as a plant biodiversity hotspot (Molina & al., 2015), including aquatic plants (Chappuis & al., 2012). But over the last decades a significant transformation seems to have occurred in some important Iberian inland aquatic ecosystems. In essence, we have observed an expansion of some non-native aquatic plants and the decrease in some other native ones (Cirujano & al., 2014). The aim of this study is to predict the potential priority risk areas for invasion of aquatic plants in the Iberian Peninsula. To accomplish this objective we have employed a species distribution model. We firstly determined the influence of environmental and socio-economic factors over 20 non-native aquatic macrophytes at a global scale. Secondly, we overlapped the individual models to achieve a map that shows the higher vulnerable areas, due to the effect of multiple invasions. Finally, we compared the most vulnerable regions with the irrigated agricultural areas in order to find an explanation for the distribution of the studied species.
MATERIAL AND METHODSStudy area
The Iberian Peninsula is located in the southwestern Europe. It is restricted by the Atlantic Ocean and the Mediterranean Sea. The Pyrenees separate it from the rest of Europe, and the Strait of Gibraltar from Africa. The climate diversity of the study area and the rugged topography of the land along with the geographic isolation, are key elements to develop an outstanding biodiversity (López-López & al., 2011). Concerning aquatic plants, this territory shows a high diversity of aquatic ecosystems and water bodies. Thus, we can find several types of rivers, streams, creeks, lakes, ponds —temporary or permanent—, bogs, and marshlands.
Species selection
We have modeled the distribution of 20 non-native aquatic macrophytes which are currently established in the Iberian Peninsula (Table 1). The non-native species belong to 13 genus and 9 families and were selected from Cirujano & al. (2014) complemented with the European and Mediterranean Plant Protection Organization list —EPPO, see http://www.eppo.int — and the Delivering Alien Invasive Species Inventories for Europe list —DAISIE, see http://www.europe-alien.org.
List of non-native aquatic macrophytes selected for the Iberian Peninsula.
Family
Genus
Species
Azollaceae
Azolla
A. filiculoides Lam. (incl. A. caroliniana Willd.)
Araceae
Lemna
L. minuta Kunth
L. valdiviana Phil.
Pistia
P. stratiotes L.
Haloragaceae
Myriophyllum
M. aquaticum (Vell.) Verdc.
M. heterophyllum Michx.
Hydrocharitaceae
Egeria
E. densa Planch.
Elodea
E. canadensis Michx.
Lagarosiphon
L. major (Ridley) Moss ex Wager
Nymphaeaceae
Nymphaea
N. mexicana Zucc.
Onagraceae
Ludwigia
L. grandiflora (Michx.) Greuter & Burdet
L. peploides subsp. montevidensis (Spreng.) P.H. Raven
L. repens J.R. Forst.
Pontederiaceae
Eichhornia
E. crassipes (Mart.) Solms
Heteranthera
H. limosa (Sw.) Willd.
H. reniformis Ruiz & Pav.
H. rotundifolia (Kunth) Griseb.
Apiaceae
Hydrocotyle
H. ranunculoides L. f.
H. verticillata Thunb.
Salviniaceae
Salvinia
S. natans (L.) All.
The global spatial occurrences of 20 species were obtained from the Global Biodiversity Information Facility (GBIF, 2015). We tested the Iberian Peninsula occurrences with data showed by the Anthos Project (Anthos, 2015). The case of Ludwigia peploides subsp. montevidensis (Spreng.) P.H. Raven was checked in other additional sources (Verloove & Sánchez, 2008; Bou & Font, 2016). Records were considered from 1950 to the present to match the timeframe for the current climate data. In order to avoid underestimating the potential niche we counted all occurrences available for each species, showing the native and invasive ranges of species (Jiménez-Valverde & al., 2011).
We used the statistical software R (R Development Core Team, 2014) to clean data and removed duplicates, data without date, and erroneous occurrences in both taxonomic and geographic data. Furthermore, we also reduced the spatial autocorrelation of the data to not violate the assumption of independence (Heffner & al., 1996). Thus, the distance between data pairs was reduced to 10 km; the same distance was used for modeling the species’ potential distribution.
Predictor variables
The 19 bioclimatic layers and altitude —Digital Elevation Model, DEM— were taken from WorldClim-Global Climate Data (Hijmans & al., 2005; Worldclim, 2015). The resolution of the environmental layers used was 5 arc-min —~10 km at the equator.
Slope was derived from DEM layer using the software ArcGIS 9.3.1 (ESRI, 2008). The human footprint was considered a socio-economic factor that reflects the human influence on the territory following Sanderson & al. (2002). This authors used as proxies of this footprint several variables such as various human land uses, population density or distance to major roads, railways and rivers. The information was obtained from Socioeconomic Data and Applications Center (SEDAC, 2015) and its resolution is 30 arc-sec —~1 km.
The resolution of 22 variables (Table 2) was turned into 5 arc-min and was projected using the World Geodetic System 1984 projection. The spatial correlation between variables was analyzed by Raster package (Hijmans & van Etten, 2015). After obtaining the correlation tree, the variables were selected by a threshold limit of 0.5. In addition, to remove the linear combination between variables in the model, the Variance Inflation Factor —VIF— was calculated using the package HH, and taking 5 as limit value (Heiberger, 2015).
List and description of used variables.
Variable
Description
Bio 1
Annual mean temperature
Bio 2
Mean diurnal range [mean of monthly (max temp – min temp)]
Bio 3
Isothermality [(Bio 2 / Bio 7) * 100]
Bio 4
Temperature seasonality
Bio 5
Maximum temperature of warmest month
Bio 6
Minimum temperature of coldest month
Bio 7
Temperature annual range (Bio 5 – Bio 6)
Bio 8
Mean temperature of wettest quarter
Bio 9
Mean temperature of driest quarter
Bio 10
Mean temperature of warmest quarter
Bio 11
Mean temperature of coldest quarter
Bio 12
Annual precipitation
Bio 13
Precipitation of wettest month
Bio 14
Precipitation of driest month
Bio 15
Precipitation seasonality (coefficient of variation)
Bio 16
Precipitation of wettest quarter
Bio 17
Precipitation of driest quarter
Bio 18
Precipitation of warmest quarter
Bio 19
Precipitation of coldest quarter
DEM
Digital Elevation Model
Slope
Slope
HFP
Human Footprint
Species distribution modeling
We developed the species distribution models with the machine learning MaxEnt version 3.3.3.k (Phillips & al., 2006), which estimates species distribution by the principle of maximum entropy. This method was chosen because is one of the most effective species distribution model, and shows a high quality achievement with low sample sizes and moderate georeferencing errors (Elith & al., 2006; Wisz & al., 2008; Mateo & al., 2010).
The parameters employed for this study were taken from Phillips & al. (2006), Phillips & Dudík (2008), and Elith & al. (2011). Default parameters were convergence threshold = 0.00001, maximum iterations = 1,000, and prevalence = 0.5, multiple regularization —default is 1— was changed to 2.5 to reduce the probability of overfitting models following Elith & al. (2010). Models were fitted with the 70% occurrences data and the remaining 30% was used to evaluate the obtained models. Besides, we used 10-fold cross-validations to estimate the errors around the fitted functions and the predictive performance on the held-out data (Elith & al., 2011). We created 10,000 background points to simulate pseudo-absences (Phillips & Dudík, 2008; Elith & al., 2011). Likewise, we interpreted the logistic output as a habitat suitability map for each species. The model accuracy was estimated using the area under the receiving operating characteristic —ROC— curve —AUC—. According to it, the results within a value of 0.5 do not discriminate better than the random, while a model with a perfect discrimination would have an AUC of 1, and values bigger or equal than 0.7 correspond to the highest predictive models (Hosmer & Lemeshow, 2000). Finally, we calculated the AUC for each model and determined the average AUC for each set of 10 replicates (Barnes & al., 2014). 10th percentile training presence threshold was chosen because it shows a good ability to predict correctly the presence of invasive species (Pearson & al., 2007; Reshetnikov & Ficetola, 2011), representing the species distribution in suboptimal habitats (Kelly & al., 2014).
Invasion risk map
The invasion risks map was calculated by overlaying the 20 species distribution individual models (Aranda & Lobo, 2011; Fajardo & al., 2014) using the Geographic Information System ArcGIS 9.3.1 (ESRI, 2008). Thereby we obtained a cartography that reflects the cumulative risk of invasion, which represent the most favorable areas for colonization and spread for the studied species in the Iberian Peninsula.
RESULTS
A total of 8,892 records were used for modeling the global potential distribution of species. The number of records varied widely among species —Nymphaea mexicana Zucc. minimum global occurrence points: 46 and Azolla filiculoides Lam. maximum occurrence points: 1,617, after cleaning data—. Fig. 1 shows the number of records per decade and the accumulated number of records per decade, and Fig. 2 shows the current presences of studied species on the Iberian Peninsula.
Number of records per decade and accumulated number of records per decade of all the aquatic macrophytes studied in the Iberian Peninsula.
Potential distribution models for the selected species: a, Azolla filiculoides; b, Egeria densa; c, Eichhornia crassipes; d, Elodea canadensis; e, Heteranthera limosa; f, Heteranthera reniformis; g, Heteranthera rotundifolia; h, Hydrocotyle ranunculoides; i, Hydrocotyle verticillata; j, Lagarosiphon major; k, Lemna minuta; l, Lemna valdiviana; m, Ludwigia grandiflora; n, Ludwigia peploides subsp. montevidensis; o, Ludwigia repens; p, Myriophyllum aquaticum; q, Myriophyllum heterophyllum; r, Nymphaea mexicana; s, Pistia stratiotes; t, Salvinia natans. Darker areas correspond with higher suitability areas; red spots indicate the presence of occurrences of the studied species in the Iberian Peninsula —after data cleaning process.
The final factors included as predictors in MaxEnt were mean diurnal range —Bio 2—, temperature seasonality (Bio 4), annual precipitation —Bio 12—, precipitation seasonality —Bio 15—, precipitation in the driest quarter —Bio 17—, altitude, slope, and human footprint —HFP.
In Table 3 we show the main results for each studied species. The accuracy scores of models ranged between 0.918 and 0.981, which shows that our models provide a good performance (Hosmer & Lemeshow, 2000) indicating a better discrimination than random chance for the species analyzed (Phillips & al., 2006). The binomial test of omission showed statistical significance —p<0.001— for each of the 10 replicates (Phillips & al., 2006), supporting the reliability of the models. The use of 10th percentile training presence threshold allowed us to discriminate correctly the presence of non-native species (Pearson & al., 2007; Reshetnikov & Ficetola, 2011) in both optimal and suboptimal areas (Jiménez-Valverde & al., 2011; Kelly & al., 2014).
AUC values ± SD and percent contribution of each of the variables taken into account for the models. In bold the best factor in the potential distribution of each species.
Variables
Species
AUC ± SD
Bio 2
Bio 4
Bio 12
Bio 15
Bio 17
Altitude
Slope
HFP
A. filiculoides
0.923 ± 0.003
0.3
43.9
4.2
18.6
1.6
1.2
0.0
29.6
E. densa
0.956 ± 0.009
1.5
20.3
2.4
1.9
12.6
1.5
1.1
58.7
E. crassipes
0.918 ± 0.014
2.6
33.2
15.4
1.2
4.3
6.7
0.4
36.3
E. canadensis
0.919 ± 0.008
0.2
20.3
5.2
11.4
25.9
0.8
0.1
36.1
H. limosa
0.956 ± 0.016
18.4
15.6
15.4
8.4
8.5
2.7
1.7
29.4
H. ranunculoides
0.940 ± 0.014
1.9
26.7
5.4
5.0
1.2
4.0
0.9
55.0
H. reniformis
0.952 ± 0.009
1.9
30.5
21.9
2.5
5.1
0.2
3.8
34.0
H. rotundifolia
0.960 ± 0.012
26.6
11.0
21
5.5
9.1
2.4
2.2
22.3
H. verticillata
0.947 ± 0.011
9.8
49.4
1.4
6.7
3.2
9.3
2.0
18.3
L. major
0.971 ± 0.004
16.2
27.9
4.0
19.8
6.7
0.2
0.1
25.1
L. minuta
0.944 ± 0.007
4.3
19.5
4.0
28.0
7.7
3.5
0.0
32.9
L. valdiviana
0.932 ± 0.031
18.5
11.1
2.9
1.0
6.1
0.4
6.4
53.7
L. grandiflora
0.981 ± 0.005
1.1
22.4
1.0
15.0
16.7
6.9
0.5
36.4
L. peploides subsp. montevidensis
0.936 ± 0.014
6.7
32.0
2.1
2.9
4.2
3.8
1.2
47.2
L. repens
0.937 ± 0.029
11.8
21.5
1.3
0.7
1.2
3.5
2.6
57.4
M. aquaticum
0.948 ± 0.005
0.5
27.6
2.1
2.1
20.6
5.7
0.2
41.2
M. heterophyllum
0.973 ± 0.012
3.9
14.0
20.9
13.5
23.7
2.9
0.7
20.4
N. mexicana
0.967 ± 0.031
4.2
24.0
1.2
1.8
1.7
0.4
0.6
66.0
P. stratiotes
0.919 ±0.010
1.1
39.0
26.0
1.4
0.3
13.4
0.5
18.3
S. natans
0.966 ± 0.013
6.8
16.8
4.8
4.5
22.9
1.5
1.1
41.6
The best predictor of potential distribution for the majority of the species was the human footprint. In relation to Azolla filiculoides, Hydrocotyle verticillata Thunb., Lagarosiphon major (Ridl.) Moss ex Wager, and Pistia stratiotes L., the best predictor was the temperature seasonality; for Heteranthera rotundifolia (Kunth) Griseb. the mean diurnal range, and for Myriophyllum heterophyllum Michx. the precipitation in the driest quarter. Besides, for these species the human footprint was included among the three best predictors (Table 3).
The suitable habitat models for the invasion risk varied broadly between species (Fig. 2), showing a large favorable distribution for species as Azolla filiculoides, Egeria densa Planch., Elodea canadensis Michx., Lemna valdiviana Phil., Nymphaea mexicana, and Ludwigia repens J.R. Forst.
The combination of the 20 individual models is the risk map for non-native Iberian aquatic macrophytes (Fig. 3). It shows the suitability of presence of the species according to the factors selected in the model building. The most vulnerable areas coincide with the littoral fringe, the high population density sectors, and the large river basins.
Invasion risk map representing the risk suitability of 20 non-native aquatic macrophytes species in the Iberian Peninsula.
Fig. 4 shows the overlapping between the irrigated agricultural areas taken from European Environment Agency (2015) and the most vulnerable region in the invasion risk map.
Map showing the irrigated agricultural areas —black polygons— over suitable habitats for 20 non-native aquatic macrophyte species.
DISCUSSION
Our results show the first geographical representation of the potential invasion risk by non-native aquatic macrophytes in the Iberian Peninsula. The combination of both environmental and socio-economic factors allows us to identify those areas more susceptible to be invaded by non-native aquatic plants.
Large areas of the Iberian Peninsula were suitable to the invasion by different non-native aquatic macrophytes, like Azolla filiculoides, Egeria densa, Elodea canadensis, Lemna valdiviana, Ludwigia repens, Myriophyllum aquaticum (Vell.) Verdc., and Nymphaea mexicana (Fig. 2). Most of them are widely distributed in Europe, being Azolla filiculoides and Elodea canadensis the species present in more European countries (Hussner, 2012).
Temperature seasonality and precipitation in the driest quarter are key factors in the probability distribution of the studied species. This result is supported by the fact that the climatic characteristics of an area act as key elements for a successful colonization of non-native species (Thuiller & al., 2005; Broennimann & al., 2007). For instance, the temperature could limit the survival, growth, and reproduction in plants (Woodward & Willians, 1987), and the precipitation in the driest quarter is associated to water availability of water bodies (Reshetnikov & Ficetola, 2011), which acts as the principal factor for the persistence of aquatic plants communities. Similar results were obtained by others authors (Gallardo & Aldridge, 2013; Barnes & al., 2014; Kelly & al., 2014), implying that non-native aquatic macrophytes are able to tolerate a wide range of environmental conditions —v.gr., seasonality in Mediterranean environments— and extreme events. This ability benefits them versus native species (Rahel & Olden, 2008; Gallardo & Aldridge, 2013). Several authors (Pearson & Dawson, 2003; Broennimann & al., 2007; Walther & al., 2009) have suggested that shifts in climate could benefit non-native species, which often tolerate temperature and precipitation ranges broader than the native ones.
The human footprint was positively associated with the presence of all studied species. This association reflects the easiness these species have to establish in disturbed habitats (Chytrý & al., 2009; Kelly & al., 2014), due to the increased presence of introduction vectors and pathways like as channels, roads or railways by which these species can be introduced and the disturbances in land uses in the studied area by human activity (Catford & al., 2011; Gallardo & Aldridge, 2013).
For example, the increase of nutrients on watercourses and water bodies, which contributes to the growth of algal blooms and the rise of turbidity levels (Carter & Rybicki, 1990; Santamaría & al., 1996) is associated with human activities. It provokes the reduction of light and oxygen availability, stopping the growth of the submerged vegetation (Moss, 1990) but enhancing floating aquatic macrophytes (Egerston & al., 2004). The new ecological scheme will promote the establishment of non-native macrophytes, which are able to colonize degraded habitats, where native macrophytes are unable to survive (Quinn & al., 2001; Catford & Downes, 2010; Chappuis & al., 2011).
Areas under the highest risk of multiple invasions include large rivers basins, highly populated areas, and the coastline (Fig. 3). An important part of the areas for colonization and expansion of these non-native species coincide with territories with agricultural development increase over the last decades. From 1970, the number of records of non-native species in the Iberian Peninsula began to rise (Fig. 1). This period overlaps with the industrialization of agriculture —the Green Revolution— when traditional non-irrigated farming was transformed into huge irrigation areas (Ruiz & al., 2008) in the Iberian Peninsula.
In this period, the high dependence on agricultural chemicals has affected freshwater ecosystems (Galil & al., 2007). Hydrological alterations and the increase of dissolved nutrients, have contributed to the eutrophication of aquatic ecosystems (Chappuis & al., 2011; Quinn & al., 2011), and the intensive land use has favored sedimentation events (Allan, 2004). All these changes have facilitated the expansion of non-native aquatic macrophytes (Egertson & al., 2004; Chappuis & al., 2011; Quinn & al., 2011). Moreover, the increment of sedimentation events caused by an intensive land use also benefits submerged non-native species. Principal areas of irrigated agriculture in the Iberian Peninsula overlap with the most susceptible areas to be invaded by non-native macrophytes (Fig. 4). This phenomenon has been reported previously by García-Murillo & al. (2007) and Ruiz & al. (2008) for Azolla filiculoides and Eichhornia crassipes (Mart.) Solms expansion, respectively. Both studies support the hypothesis, together with ours, that the quick expansion of non-native macrophytes is due to the nutrients increase contributed by adjacent agricultural areas.
In addition, we also have observed that some areas predicted as being suitable (Fig. 3) were currently unoccupied —see Fig. 2, current presences of studied species—. This may be due to different causes: areas where species have been successfully eradicated —v.gr., Pistia stratiotes in neighborhood Doñana National Park, Southern Spain, as pointed up by García-Murillo & al. (2005)— or areas with geographical barriers or species interactions that limited its distributions —v.gr., Azolla filiculoides has not been detected in temporary ponds and marshes in Doñana National Park while the weevil Stenopelmus rufinasus Gyllenhal was present in samples, as pointed up by Florencio & al. (2015)—. Besides, they can also be areas where species have not been detected yet due to the lack of studies in these places, or because this species may have not been able to colonize these suitable areas yet (Liu & al., 2011) as a consequence of they are still in the early stages of the invasion process. These two last points are crucial for proper management and early control of non-native species.
Among the species studied in this work, we consider that the most harmful are Azolla filiculoides and Eichhornia crassipes, both present in the major part of the World, being the two more potentially invasive species in Europe and the Mediterranean basin (Hussner, 2012; Kriticos & Brunel, 2016). Their invasion capacity is due not only to climate tolerance and the adapting ability to eutrophic environments, but also to a high rate of vegetative reproduction that ensure the success of colonization in invaded habitats and a high competition with others species (Ruiz & al., 2008; Fernández-Zamudio & al., 2013).
In conclusion, our study, based on the global distribution of 20 non-native aquatic macrophyte species, contributes to the understanding of the distribution patterns of non-native aquatic macrophytes in the Iberian Peninsula, and it may be used as a base to develop useful tools to manage successfully the Iberian biodiversity in future conservation planning, and for the conservation and management of aquatic ecosystems in other lands. Species distribution models should not be a substitute for field work, but they are a first step that allows an early identification of the most vulnerable areas to implement more effective management efforts preventing biological invasions.
ACKNOWLEDGEMENTS
We are indebted to Professor Timothy H Keitt and Dr. Laura I. González, from the University of Texas, for their comments of the manuscript and for the English revision. Also, we thanks to associate editor and two anonymous reviewers which provided helpful recommendations that improved substantially our manuscript.
REFERENCESAguiarF.C.F.FerreiraM.T.Plant invasions in the rivers of the Iberian Peninsula, south-western Europe: A review201314711071119https://doi.org/10.1080/11263504.2013.861539Anthos2015Real Jardín Botánico de Madrid CSIC-Fundación Biodiversidad[http://www.anthos.es/]AllanJ.D.Landscapes and riverscapes: The influence of land use on stream ecosystems200435257284https://doi.org/10.1146/annurev.ecolsys.35.120202.110122ArandaS.C.LoboJ.M.How well does presence-only-based species distribution modelling predict assemblage diversity? A case study of the Tenerife flora2011343138https://doi.org/10.1111/j.1600-0587.2010.06134.xBarnesM.A.JerdeC.L.WittmannM.E.ChaddertonW.L.DingJ.ZhangJ.PurcellM.BudhathokiM.LodgeD.M.Geographic selection bias of occurrence data influences transferability of invasive Hydrilla verticillata distribution models2014425842593https://doi.org/10.1002/ece3.1120BouJ.FontJ.Situation in Catalonia of Ludwigia peploides (Onagraceae)2016805758BroennimannO.GuisanA.Predicting current and future biological invasions: both native and invaded ranges matter20084585589https://doi.org/10.1098/rsbl.2008.0254BroennimannO.TreierU.A.Muller-ScharerH.ThuillerW.PetersonA.T.GuisanA.Evidence of climatic niche shift during biological invasion200710701709https://doi.org/10.1111/j.1461-0248.2007.01060.xBrooksM.L.D’AntonioC.M.RichardsonD.M.GraceJ.B.KeeleyJ.E.DiTomasoJ.M.HobbsR.J.PellantM.PykeD.Effects of invasive alien plants on fire regimes200454677688https://doi.org/10.1641/0006-3568(2004)054[0677:EOIAPO]2.0.CO;2BrunduG.Plant invaders in European and Mediterranean inland waters: profiles, distribution, and threats20157466179https://doi.org/10.1007/s10750-014-1910-9CarterV.RybickiN.B.Light attenuation and submersed macrophyte distribution in the tidal Potomac River and estuary199013441452https://doi.org/10.2307/1351788CatfordJ.A.DownesB.J.Using multi-scale species distribution data to infer drivers of biological invasion in riparian wetlands2010162032https://doi.org/10.1111/j.1472-4642.2009.00631.xCatfordJ.A.VeskP.A.WhiteM.D.WintleB.A.Hotspots of plant invasion predicted by propagule pressure and ecosystem characteristics20111710991110https://doi.org/10.1111/j.1472-4642.2011.00794.xChambersP.A.LacoulP.MurphyK.J.ThomazS.M.Global diversity of aquatic macrophytes in freshwater2008595926https://doi.org/10.1007/s10750-007-9154-6ChappuisE.GaciaE.BallesterosE.Changes in aquatic macrophyte flora over the last century in Catalan water bodies (NE Spain)201195268277https://doi.org/10.1016/j.aquabot.2011.08.006ChappuisE.BallesterosE.GaciaE.Distribution and richness of aquatic plants across Europe and Mediterranean countries: patterns, environmental driving factors and comparison with total plant richness201223985997https://doi.org/10.1111/j.1654-1103.2012.01417.xChytrýM.PyšekP.WildJ.PinoJ.MaskellL.C.VilàM.European map of alien plant invasions based on the quantitative assessment across habitats20091598107https://doi.org/10.1111/j.1472-4642.2008.00515.xCirujanoS.MecoA.García-MurilloP.Chirino ArgentaM.2014Real Jardín Botánico de Madrid CSICMadridCollenB.WhittonF.DyerE.E.BaillieJ.E.M.CumberlidgeN.DarwallW.R.T.PollockC.RichmanN.I.SoulsbyA.M.BöhmM.Global patterns of freshwater species diversity, threat and endemism2014234051https://doi.org/10.1111/geb.12096CraftonR. E.Modeling invasion risk for coastal marine species utilizing environmental and transport vector data2015746349362https://doi.org/10.1007/s10750-014-2027-xDibbleE.D.ThomazS.M.PadialA.A.Spatial complexity measured at a multi-scale in three aquatic plant species200621239247https://doi.org/10.1080/02705060.2006.9664992EgertsonC.J.KopaskaJ.A.DowningJ.A.A century of change in macrophyte abundance and composition in response to agricultural eutrophication2004524145156https://doi.org/10.1023/B:HYDR.0000036129.40386.ceElithJ.LeathwickJ.R.Species distribution models: ecological explanation and prediction across space and time200940677697https://doi.org/10.1146/annurev.ecolsys.110308.120159ElithJ.KearneyM.PhillipsS.The art of modelling range-shifting species20101330342https://doi.org/10.1111/j.2041-210X.2010.00036.xElithJ.PhillipsS.J.HastieT.DudíkM.CheeY.E.YatesC.J.A statistical explanation of MaxEnt for ecologists2011174357https://doi.org/10.1111/j.1472-4642.2010.00725.xElithJ.GrahamC.H.AndersonR.P.DudikM.FerrierS.GuisanA.HijmansR.J.HuettmannF.LeathwickJ.R.LehmannA.LiJ.LohmannL.G.LoiselleB.A.ManionG.MoritzC.NakamuraM.NakazawaY.OvertonJ.M.PetersonA.T.PhillipsS.J.RichardsonK.Scachetti-PereiraR.SchapireR.E.SoberonJ.WilliamsS.WiszM.ZimmermannN.E.Novel methods improve prediction of species’ distributions from occurrence data200629129151https://doi.org/10.1111/j.2006.0906-7590.04596.xESRI2008Environmental Systems Research Institute, United States of AmericaEuropean Environment Agency2015[http://www.eea.europa]FajardoJ.LessmannJ.BonaccorsoE.DevenishC.MuñozJ.Combined use of systematic conservation planning, species distribution modeling, and connectivity analysis reveals severe conservation gaps in a megadiverse country (Peru)20149123https://doi.org/10.1371/journal.pone.0114367Fernández-ZamudioR.CirujanoS.Sánchez-CarrilloS.MecoA.García-MurilloP.Clonal reproduction of Azolla filiculoides Lam.: implications for invasiveness201332245252FlorencioM.Fernández-ZamudioR.BiltonD.T.Díaz-PaniaguaC.The exotic weevil Stenopelmus rufinasus Gyllenhal, 1835 (Coleoptera: Curculionidae) across a “host-free” pond network2015347984GalilB.S.NehringS.PanovV.2007W. Nentwig, SpringerBerlinGallardoB.AldridgeD.C.The dirty dozen: socio-economic factors amplify the invasion potential of 12 high-risk aquatic invasive species in Great Britain and Ireland201350757766https://doi.org/10.1111/1365-2664.12079GallardoB.ErreaM.P.AldridgeD.Application of bioclimatic models coupled with network analysis for risk assessment of the killer shrimp, Dikerogammarus villosus, in Great Britain20121412651278https://doi.org/10.1007/s10530-011-0154-0GallardoB.ClaveroM.SánchezM.I.VilàM.Global ecological impacts of invasive species in aquatic ecosystems201522113García-MurilloP.DanaE.D.RodríguezC.Pistia stratiotes L. (Araceae) Una planta acuática exótica en las proximidades del Parque Nacional de Doñana (SW España)200530235236García-MurilloP.Fernández-ZamudioR.CirujanoS.SousaA.EspinarJ.M.The invasion of Doñana National Park (SW Spain) by the mosquito fern (Azolla filiculoides Lam.)200726242250García-MurilloP.Fernández-ZamudioR.Díaz PaniaguaC.Las plantas de las lagunas temporales de Doñana2015Organismo Autónomo Parques NacionalesMadridGBIF2015[http://www.gbif.org/]HarrelS.DibbleE.Foraging efficiency of juvenile bluegill, Lepomis macrochirus, among different vegetated habitats200162441453https://doi.org/10.1023/A:1012259922727HeffnerR.A.ButlerM.J.ReillyC.K.Pseudoreplication revisited19967725582562https://doi.org/10.2307/2265754HeibergerR.M.2015[http://cran.r-project.org]HijmansR.J.Van EttenJ.2015[http://cran.r-project.org]HijmansR.J.CameronS.E.ParraJ.L.JonesP.G.JarvisA.Very high resolution interpolated climate surfaces for global land areas20052519651978https://doi.org/10.1002/joc.1276HosmerD.W.LemeshowS.2000John Wiley & Sons Inc.New Yorkhttps://doi.org/10.1002/0471722146HulmeP.E.Beyond control: wider implications for the management of biological invasions200643835847https://doi.org/10.1111/j.1365-2664.2006.01227.xHussnerA.Alien aquatic plant species in European countries201252297306https://doi.org/10.1111/j.1365-3180.2012.00926.xJiménez-ValverdeA.PetersonA.T.SoberónJ.OvertonJ.M.AragónP.LoboJ.M.Use of niche models in invasive species risk assessments20111327852797https://doi.org/10.1007/s10530-011-9963-4KellyR.LeachK.CameronA.MaggsC.A.ReidN.Combining global climate and regional landscape models to improve prediction of invasion risk201420111https://doi.org/10.1111/ddi.12194KriticosD.J.BrunelS.Assessing and managing the current and future pest risk from water hyacinth, (Eichhornia crassipes), an invasive aquatic plant threatening the environment and water security201611118https://doi.org/10.1371/journal.pone.0120054LesD.H.CrawfordD.J.KimballR.T.MoodyM.L.LandoltE.Biogeography of discontinuously distributed hydrophytes: A molecular appraisal of intercontinental disjunctions2003164917932https://doi.org/10.1086/378650LiuX.GuoZ.KeZ.WangS.LiY.Increasing potential risk of a global aquatic invader in Europe in contrast to other continents under future climate change20116111https://doi.org/10.1371/journal.pone.0018429LoweS.J.BrowneM.BoudjelasS.de PooterM.2004The Invasive Species Specialist Group (ISSG, SSC, IUCN)MackR.N.SimberloffD.Mark LonsdaleW.EvansH.CloutM.BazzazF.A.Biotic invasions: causes, epidemiology, global consequences, and control200010689710https://doi.org/10.1890/1051-0761(2000)010[0689:BICEGC]2.0.CO;2MateoR.G.CroatT.B.FelicísimoA.M.MuñozJ.Profile or group discriminative techniques? Generating reliable species distribution models using pseudo-absences and target-group absences from natural history collections2010168494https://doi.org/10.1111/j.1472-4642.2009.00617.xMolinaR.AparicioA.LavergneS.SlingsbyJ.ArroyoJ.Investigating the evolutionary assembly of a Mediterranean biodiversity hotspot: deep phylogenetic sgnal in the distribution of eudicots across elevational belts201542507518https://doi.org/10.1111/jbi.12398MossB.Engineering and biological approaches to the restoration from eutrophication of shallow lakes in which aquatic plant communities are important components1990200-201367377https://doi.org/10.1007/BF02530354MurphyK.J.Plant communities and plant diversity in softwater lakes of northern Europe200273287324https://doi.org/10.1016/S0304-3770(02)00028-1López-LópezP.MaioranoL.FalcucciA.BarbaE.BoitaniL.Hotspots of species richness, threat and endemism for terrestrial vertebrates in SW Europe201137399412https://doi.org/10.1016/j.actao.2011.05.004PearsonR.G.DawsonT.P.Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful?200312361371https://doi.org/10.1046/j.1466-822X.2003.00042.xPearsonR.G.RaxworthyC.J.NakamuraM.Townsend PetersonA.Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar200734102117https://doi.org/10.1111/j.1365-2699.2006.01594.xPetersonA.T.Predicting the geography of species’ invasions via ecological niche modeling20037841933https://doi.org/10.1086/378926PhillipsS.J.DudíkM.Modeling of species distributions with MaxEnt: new extensions and a comprehensive evaluation200831161175https://doi.org/10.1111/j.0906-7590.2008.5203.xPhillipsS.J.AndersonR.P.SchapireR.E.Maximum entropy modeling of species geographic distributions2006190231259https://doi.org/10.1016/j.ecolmodel.2005.03.026PimentelD.ZunigaR.MorrisonD.Update on the environmental and economic costs associated with alien-invasive species in the United States200552273288https://doi.org/10.1016/j.ecolecon.2004.10.002PyšekP.RichardsonD.M.Invasive species, environmental change and management, and health2010352555https://doi.org/10.1146/annurev-environ-033009-095548QuinnL.D.SchoolerS.S.van KlinkenR.D.Effects of land use and environment on alien and native macrophytes: lessons from a large-scale survey of Australian rivers201117132143https://doi.org/10.1111/j.1472-4642.2010.00726.xR Development Core Team2014Foundation for Statistical ComputingVienna, Austria[http://cran.r-project.org]RahelF.J.OldenJ.D.Assessing the effects of climate change on aquatic invasive species200822521533https://doi.org/10.1111/j.1523-1739.2008.00950.xRennieM.D.JacksonL.J.The influence of habitat complexity on littoral invertebrate distributions: patterns differ in shallow prairie lakes with and without fish20056220882099https://doi.org/10.1139/f05-123ReshetnikovA.FicetolaG.Potential range of the invasive fish rotan (Perccottus glenii) in the Holarctic20111329672980https://doi.org/10.1007/s10530-011-9982-1RicciardiA.KippR.Predicting the number of ecologically harmful exotic species in an aquatic system200814374380https://doi.org/10.1111/j.1472-4642.2007.00451.xRuizT.MartínE.LorenzoG.AlbanoE.MoránR.SánchezJ.M.The Water Hyacinth, Eichhornia crassipes: an invasive plant in the Guadiana River Basin (Spain)200834253https://doi.org/10.3391/ai.2008.3.1.8SalaO.E.ChapinF.S.ArmestoJ.J.BerlowE.BloomfieldJ.DirzoR.Huber-SanwaldE.HuennekeL.F.JacksonR.B.KinzigA.LeemansR.LodgeD.M.MooneyH.A.OesterheldM.PoffN.L.SykesM.T.WalkerB.H.WalkerM.WallD.H.Global biodiversity scenarios for the year 2100200028717701774https://doi.org/10.1126/science.287.5459.1770SandersonE.W.JaitehM.LevyM.A.RedfordK.H.WanneboA.V.WoolmerG.The human footprint and the last of the wild200252891904https://doi.org/10.1641/0006-3568(2002)052[0891:THFATL]2.0.CO;2SantamaríaL.Why are most aquatic plants widely distributed? Dispersal, clonal growth and small-scale heterogeneity in a stressful environment200223137154https://doi.org/10.1016/S1146-609X(02)01146-3SantamaríaL.MontesC.HootsmansM.J.M.Influence of environmental parameters on the biomass development of Ruppia drepanensis populations in Doñana National Park: the importance of conditions affecting the underwater light climate19965157180https://doi.org/10.1007/BF01995828SEDAC2015[http://sedac.ciesin.columbia.edu]SerranoL.Díaz PaniaguaC.Díaz PaniaguaC.Introducción2015Organismo Autónomo Parques NacionalesMadrid918StrayerD.L.DudgeonD.Freshwater biodiversity conservation: recent progress and future challenges201029344358https://doi.org/10.1899/08-171.1ThuillerW.RichardsonD.M.PyšekP.MidgleyG.F.HughesG.O.RougetM.Niche-based modelling as a tool for predicting the risk of alien plant invasions at a global scale20051122342250https://doi.org/10.1111/j.1365-2486.2005.001018.xVerlooveF.Sánchez-GullónE.New records of interesting xenophytes in the Iberian Peninsula200833147167WaltherG.R.RoquesA.HulmeP.E.SykesM.T.PysekP.KuhnI.ZobelM.BacherS.Botta-DukatZ.BugmannH.CzuczB.DauberJ.HicklerT.JarosikV.KenisM.KlotzS.MinchinD.MooraM.NentwigW.OttJ.PanovV.E.ReinekingB.RobinetC.SemenchenkoV.SolarzW.ThuillerW.VilaM.VohlandK.SetteleJ.Alien species in a warmer world: risks and opportunities200924686693https://doi.org/10.1016/j.tree.2009.06.008WillbyN.J.Managing invasive aquatic plants: problems and prospects200717659665https://doi.org/10.1002/aqc.913WilliamsS.GrosholzE.The invasive species challenge in estuarine and coastal environments: marrying management and science200831320https://doi.org/10.1007/s12237-007-9031-6WiszM.S.HijmansR.J.LiJ.PetersonA.T.GrahamC.H.GuisanA.NCEAS Predicting Species Distribution Working GroupEffects of sample size on the performance of species distribution models200814763773https://doi.org/10.1111/j.1472-4642.2008.00482.xWoodwardF.I.WilliamsB.G.Climate and plant distribution at global and local scales198769189197https://doi.org/10.1007/BF00038700Worldclim2015[http://www.worldclim.org]