Journal Information
Share
Share
Download PDF
More article options
Visits
266
Research Letters
Full text access
Available online 18 November 2024
Mapping the threat: projecting invasive plant distribution in the tropical Andes under climate change
Visits
266
Juan David González-Trujilloa,b,
Corresponding author
jdgonzalezt@unal.edu.co

Corresponding author.
, Magda R. Escobar-Albab, David E. Larab, Juan E. Carvajal-Cb
a Departamento de Biología, Facultad de Ciencias, Universidad Nacional de Colombia, Sede Bogotá, Colombia
b Escuela de Biología, Universidad Pedagógica y Tecnológica de Colombia, Tunja, Colombia
Highlights

  • Tropical mountains provide favorable habitats for invasive plant species.

  • Climate change is expected to open new space for the spread of invasive species.

  • Invasive species could gain between 1,000 and 10,000 km2 of more suitable areas.

  • The areas of greatest suitability are not necessarily at the highest elevations.

This item has received
Received 04 June 2024. Accepted 04 November 2024
Article information
Abstract
Full Text
Bibliography
Download PDF
Statistics
Figures (5)
Show moreShow less
Tables (1)
Table 1. List of modeled invasive plant species and their known impacts in the tropics.
Additional material (2)
Abstract

The high tropical mountains are renowned for their exceptional biodiversity, which plays a pivotal role in sustaining local communities and their economies. Nevertheless, these regions are confronted with considerable challenges, primarily due to the threat of invasive species, which is further compounded by the impacts of climate change. Here, we investigated the current and future distribution of invasive plant species in the tropical Andes, focusing on Páramos – ecosystems that provide essential water services to people. Species distribution modeling was used to assess the current and future suitability of 11 of the most harmful invasive species under low (SSP126), moderate (SSP370), and intense (SSP585) carbon emission scenarios. Our projections show that between 18% and 60% of the tropical mountain area (480,000 km2) is currently suitable for the establishment of at least one species. Additionally, projections indicate that this area could expand by 2–4 % by the mid-century, depending on the climate change scenario. This expansion could potentially make 500–3000 km2 of the Páramo extent suitable for up to five invasive species. Overall, our study provides relevant information for targeted management actions, such as manual removal and restoration, which are essential for limiting the expansion of invasive species to more suitable areas.

Keywords:
Invasion risk
Niche envelope
High mountain
Elevational change
Potential distribution
Full Text
Introduction

Tropical mountains are renowned for their high biodiversity both on land and in water. However, the increasing threat of invasive species, exacerbated by human activities and climate change, poses a major threat to this rich diversity (He et al., 2023), particularly at higher elevations (Ruiz et al., 2008). The growth of human populations at high altitudes has caused increased habitat fragmentation (He et al., 2023), enabling the establishment and expansion of invasive species (Fuentes-Lillo et al., 2023). Additionally, elevated temperatures facilitate the upward propagation of both native and invasive species (Iseli et al., 2023), which subsequently intensifies competition at the summit. Although these processes could lead to the degradation of native biodiversity, a comprehensive assessment of the potential threats posed by the climate-driven expansion of invasive species is currently lacking in tropical mountainous regions, such as the northern Andes.

The spread of invasive species can negatively impact ecosystems and human well-being. In tropical mountains, impacts range from displacement of native and cultivated plant species to soil degradation by altering natural water and nutrient cycles (Cárdenas López et al., 2017; Costa et al., 2019; Le Maitre et al., 2011; Valdez-Ramírez et al., 2020). Currently, invasive species such as Genista monspessulana ((L.) L.A.S. Johnson, 1962), Melinis minutiflora (P. Beauv., 1812), and Ulex europaeus (L., 1753) also increase the risk of fire in rural and urban landscapes by increasing fuel loads and fire corridors (Anderson and Anderson, 2009; D’Antonio et al., 2011; Pauchard et al., 2008), whereas others, such as Ricinus communis (L., 1753), corrode the integrity of roads and homes (Sandoval et al., 2022). Combined, these impacts could cost as much as USD 10 billion to governments in Andean countries (Heringer et al., 2021) and cause incalculable losses in terms of Nature Contributions to People (Roy et al., 2024). Thus, the future of biodiversity in the high mountains and the well-being of the people living in these areas (more than 150 people per square kilometer, Williams, 2011) are inextricably linked to the effective containment of invasive species.

Climate change could increase the invasion risk in high-elevation ecosystems such as the páramos above the Andean tree line. Projections indicate that changes in climate patterns will cause significant shifts in the distribution of native biodiversity at higher elevations (Peyre, 2022). Additionally, climate change may weaken historical environmental constraints, making it easier for non-native plants to establish themselves in the region (Valencia et al., 2013). Although both processes may create 'open spaces' for invasive species, it is currently unknown which species will use these newly accessible habitats and to what extent Páramo may be invaded. To address these knowledge gaps, we utilized species distribution modeling to identify the present-day and future zones of higher suitability for the establishment of the most harmful invasive plants currently found in the high Andean Mountains. In particular, we asked: (1) How extensive is the range of conditions currently suitable for invasive species in the tropical Andes? (2) How will range change under different climate change scenarios? (3) How much of the Páramo ecosystem is likely suitable for the establishment of invasive species? We investigated these questions in relation to three global circulation models and three climate change scenarios, each with low, mid, and high levels of carbon emission. As such, we aimed to capture the uncertainties in habitat suitability projections that may arise from different climate change projections.

Material and methods

Our modelling approach consists of four stages (Fig. 1): (A) preparing data of species occurrences and environmental predictors, (B) training three algorithms to predict the distribution of each species in its native region, (C) validating the usefulness of the models to project the distribution of invasive species in the tropical Andes, and (D) employing the trained models to project the potential distribution of invasive species for the mid-century. All modelling stages were conducted using R ver. 4.2.1 and implemented in accordance with a standardized and replicable methodological pipeline, which has been successfully employed with other biological groups in situations of limited data (e.g., Gouvêa et al., 2024). R scripts to reproduce this pipeline are available in a permanent repository on GitHub (https://github.com/jorgeassis/speciesDistributionModelling).

Fig. 1.

Flowchart showing the step-by-step procedure for projecting the current and future distribution of invasive plant species in the tropical Andes. (A) Data preparation included (1) the compilation and curation of species occurrences in their native regions; (2) the implementation of a Principal Component Analysis (PCA) performed on the environmental variables in the native region; (3) the application of a kernel-based filter, which splits the environmental space into two subspaces associated with either the environmental conditions more suitable for the species (in yellow) or those associated with less or unsuitable environmental conditions; and (4) uniform sampling of pseudo-across a sampling grid of a chosen resolution superimposed on the 2-dimensional environmental space. (B) Model training and evaluation included calibration of three algorithms and evaluation of their individual and ensemble performance using a set of metrics. (C) Model validation included evaluation of model transferability (i.e., performance outside the training region) using occurrences registered in the tropical Andes. (D) Ensemble prediction involved using the ensemble of three algorithms to project the distribution of each invasive plant species under three scenarios of future climate change.

(0.58MB).

Below, we provide a detailed description of the steps and data involved in each stage of the modelling approach.

  • A

    Data preparation

Occurrence dataset. Occurrence data for eleven plant invasive species (Table 1) in the Tropical Andes were compiled from various databases, including the Global Biodiversity Information Facility (GBIF, retrieved in November 2023), the Catalog of Plants of Colombia, Tropicos (www.tropicos.org), and scientific publications (e.g., Sandoval et al., 2022). To ensure data quality, we retained only records from the Tropical Andes and regions of origin for each species (Table 1) with taxonomic verification and GPS errors below 1 km. Furthermore, the retained records were uniformly sampled within a 2-dimensional environmental space created using principal component analysis (PCA) with all environmental predictors (Fig. 1A). The 'uniformSampling' function from the USE R package (Da Re et al., 2023) was used for this purpose. This environmentally driven thinning approach addresses the biases that stem from uneven sampling across the entire environmental range, as highlighted by Castellanos et al., 2019 and Da Re et al., 2023. By employing PCA to capture over 75% of the total environmental variation and thinning occurrences within this space, we were able to achieve a balanced representation across the environmental gradient, thereby reducing the oversampling of accessible areas and undersampling of remote regions. As such, this approach enhances the reliability of our distribution models, particularly regarding the projection of species distributions in future climate change scenarios (Varela et al., 2014). The final number of records used for calibrating and evaluating species distribution models is described in Table S1.

Table 1.

List of modeled invasive plant species and their known impacts in the tropics.

Species  Biogeographic origin  Impacts  References 
Acacia decurrens(J.C.Wendl.) Willd., 1813  Australian  Create regular layers of litter on the ground, inhibit or redirect successional processes, inhibit the expression of seed banks, and limit resource supply, leading to displacement of native plants and animals and increasing the frequency of fires.  (Cárdenas López et al., 2017; Le Maitre et al., 2011
Acacia melanoxylonR.Br., 1813  Australian  Alter the structure and function of their ecosystems, thereby displacing their native flora. It also causes soil erosion and alters hydrological cycles, negatively affecting agriculture.  (Kumschick and Jansen, 2023; Le Maitre et al., 2011
Arundo donaxL., 1753.  Holarctic  Alter the natural vegetation structure, outcompete native plant species and diminish the diversity and abundance of animals such as arthropods and birds. It also drives out soil, fuels forest fires, displaces native species, and increases the invasion of ticks that affect livestock.  (Cárdenas López et al., 2017; Girotto et al., 2021; Lambert et al., 2010
Genista monspessulana(L.) L.A.S. Johnson, 1962  Holarctic  Alter fire regime and nutrient cycling, displace native species and decrease native diversity by forming dense monospecific stands. It also facilitates the establishment of other invasive species and produces seeds that are toxic to livestock and humans.  (Cárdenas López et al., 2017; Herrera et al., 2016; Pauchard et al., 2008
Hedychium coronariumJ. Koenig, 1783  Indo-Malesian  Alter hydrological and nutrient cycles in soil. It forms thickets that suppress the successional and regeneration processes of native species, thus affecting the native flora and crops.  (Cárdenas López et al., 2017; Costa et al., 2019
Melinis minutifloraP. Beauv., 1812  African  Increases the occurrence of fires, displaces native species, and alters soil properties and decomposition. It also inhibits the growth of native species.  (Cárdenas López et al., 2017; Nogueira et al., 2019; Sandoval et al., 2022
Pteridium aquilinum(L.) Kuhn, 1834  Holarctic  Alter vegetation success processes affect crops and cause livestock poisoning. It also produces acids that inhibit root growth in native and cultivated species.  (Berget et al., 2015; Cárdenas López et al., 2017; Valdez-Ramírez et al., 2020
Ricinus communisL., 1753  African  Alter vegetation success processes affect crops and cause livestock poisoning. It also produces acids that inhibit root growth in native and cultivated species.  (Cárdenas López et al., 2017; Sandoval et al., 2022; Silva and Fabricante, 2022
Senecio madagascariensisPoir., 1804  African  Alter soil nutrient cycles, damage to agricultural crops, and outcompete native species. It also contains substances that are toxic to both animals and humans.  (Wijayabandara et al., 2021
Thunbergia alata(Lindl.) T. Anderson, 1866  African  Displace native species and reduce habitat heterogeneity, thereby affecting the structure and function of native ecosystems.  (Cárdenas López et al., 2017; Quijano-Abril et al., 2021
Ulex europeausL., 1753  Holarctic  Dry soil and increase the occurrence of fires. Inhibits vegetative growth, including pastures in agricultural and livestock lands.  (Anderson and Anderson, 2009; Cárdenas López et al., 2017

Pseudo-absence dataset. Given that the available species distribution data contain presence-only records, we adopted the approach of Da Re and colleagues (2023) to generate a sample of pseudo-absences. This involved three steps: recovering the core area of the species' bioclimatic niche within the environmental space using records within species native ranges; using a kernel-based filter to split the environmental space into subspaces associated with suitable and less suitable conditions; and uniformly sampling pseudo-absences across subspaces representing less suitable conditions. All three steps were implemented using the 'paSampling' function from the USE R package (Da Re et al., 2023). The native ranges of the species were delineated based on the extent of the phylogeographic region in which the evolutionary lineage of each species originated. This approach allows for the consideration of all potentially accessible areas from the lineage's center of origin and the species' dispersal potential and associated barriers, as proposed by Rojas‐Soto et al. (2024). The boundaries of each phylogeographic region were extracted from a study on floristic regions of the world by Liu et al. (2023).

Environmental Predictors. Fifteen environmental predictors were chosen based on established ecological relationships (see Table S1). These predictors were derived from various sources and encompassed climate (Brun et al., 2022), soil characteristics (Lembrechts et al., 2022; Poggio et al., 2021), and human-mediated effects (e.g., livestock density, Gilbert et al., 2018). All environmental layers were averaged to a grid scale of 1 × 1 km grid cells and reprojected to the ESPG:4326 coordinate system using the functions ‘resample’ and ‘project’ from the terra R package (Hijmans, 2024). The Tropical Andes Mountain range, spanning from central Peru to northern Colombia and Venezuela, was utilized for the projections. The Huancabamba depression, located at 6 °S–79 °W, serves as the southern limit as it denotes the transition between tropical and subtropical climates (Weigend, 2002).

  • B

    Model training and evaluation

For each species under consideration, we employed three distinct machine learning algorithms to develop species distribution models (SDMs) (Fig. 1B). The algorithms utilized were Adaptive Boosting (AdaBoost) (Krause-Jensen et al., 2020), Boosted Regression Trees (BRT) (Elith et al., 2008), and Extreme Gradient Boosting (XGBoost) (Chen and Guestrin, 2016). These algorithms are recognized for their effectiveness in the presence of limited data (Barbet‐Massin et al., 2012) as well as their capacity to capture the intricate interplay between predictors and response variables and their adaptability to altered hyperparameters (Elith et al., 2008). Each algorithm has a distinct advantage. AdaBoost is an effective method for handling imbalanced datasets, BRT can capture intricate species-environment interactions, and XGBoost enhances computational efficiency while reducing overfitting through regularization. Accordingly, within an ensemble forecasting framework (Araujo and New, 2007), these models utilize their complementary strengths, thereby enhancing their transferability and facilitating reliable species distribution predictions within and outside their native regions (Elith et al., 2008). AdaBoost models were implemented using the 'mboost' function from the mboost package (Hofner et al., 2014); BRT models were implemented using the 'gbm' function from the gbm package (Greenwell et al., 2019); and XGBoost models were implemented using the 'xgboost' function from the xgboost package (Chen and Guestrin, 2016).

Following Gouvêa et al. (2024), we used a 10-fold cross-validation framework and the grid search method to identify the optimal combination of hyperparameter values for the three algorithms. The grid for BRT changes with respect to the learning rate (0.1, 0.01, and 0.001), tree complexity (1–4), and the number of trees (50–1000, step 50). The grid for AdaBoost changes with respect to the number of interactions (50–250, step 50), degrees of freedom (1–12), and shrinkage (0.25–1, step 0.25). The grid for XGBoost changes with respect to gamma (0–5, step 1), interaction depth (1–4), shrinkage (0.1−0.5, step 0.1), and the number of rounds (10–100, step 10). We also ensured that the models had monotonic responses to positive or negative constraints (Table S2). This helped to reduce overfitting and align with the expected relationships between species occurrences and environmental predictors (Hofner et al., 2011). Prior to modeling, we examined collinearity among predictors using the variance inflation factor (VIF) metric and selected a reduced set of uncorrelated predictors through a stepwise procedure to avoid multicollinearity (Dormann et al., 2013). The threshold for selecting uncorrelated predictors was set at a VIF < 3.

The performance of the individual models and their weighted average ensemble (Araújo et al., 2006) in predicting the distribution of each species within each species' native range (see Table 1) was evaluated using the Boyce index, area under the receiver operating characteristic curve (AUC), and total score (TSS) (Boyce et al., 2002; Fielding and Bell, 1997). The values of these metrics range from −1 to +1 for the Boyce index and from 0 to 1 for the AUC and TSS. A positive Boyce index above 0, or AUC and TSS values above 0.5, indicate that the model's predictions outperform random expectations, with values approaching 1 indicating strong agreement between the predicted and observed distribution of each species (Hirzel et al., 2006).

  • C

    Model validation

As the models were trained (and evaluated) with data from regions where the species are native, it is essential to ascertain their transferability to regions where the species are considered invasive (Fig. 1C). In other words, the models should be validated by determining whether the projected distribution in the Andes chain aligns with documented records. We validated the transferability of the ensemble model using the Boyce index (Boyce et al., 2002). One advantage of the Boyce index is that it is not dependent on threshold-based metrics such as TSS and AUC. The Boyce index assesses the performance of a model across the entire suitability gradient, thereby obviating the need for a fixed cutoff. This makes it particularly well suited for situations in which reliable absence data are lacking (Hirzel et al., 2006). Moreover, because it captures the relationship between predicted suitability and observed occurrences across all suitability levels, the Boyce index provides a more comprehensive assessment of model performance and reliability (Leroy et al., 2018).

The presence data of each species were used to validate the transferability of the model to the tropical Andes. As with the occurrence dataset (see above), occurrence records across the Tropical Andes were uniformly sampled using the 'uniformsampling' function of the USE R package (Da Re et al., 2023). This was performed to prevent the inclusion of an excessive number of occurrences derived from sites with similar environmental characteristics, which would otherwise inflate metric values. The numbers of occurrences used to validate the distribution models for each species are presented in Table S1.

  • D

    Ensemble forecasting

For each species, we combined the outcome of the different models within an ensemble forecasting framework (Araujo and New, 2007) to produce maps of future (2040–2070) habitat suitability (Fig. 1D). To assess the uncertainties in projections of habitat suitability that may result from various scenarios of future climate change, we generated projections based on three shared socio-economic pathways (SSP) scenarios: SSP126, SSP370, and SSP585. These scenarios correspond to the lowest (best), mid, and highest (worst) greenhouse gas emission scenarios (Riahi et al., 2017). Each SSP scenario was produced using three different global circulation models (GFDL-ESM [Dunne et al., 2020], MPI-ESM [Gutjahr et al., 2019], IPSL-CM6A [Boucher et al., 2020]). Distribution maps are available in the Zenodo repository (DOI: 10.5281/zenodo.10907376).

Habitat suitability maps for the present and future were transformed into binary maps using a threshold that optimizes TSS (sensitivity and specificity) (Liu et al., 2016). To measure the distribution area of each species, we summed all cells where its presence was predicted. Similarly, we determined the elevational ranges of each species based on the altitude of the cells in which they were predicted, obtained from a digital elevation model (DEM). To determine whether the species' elevation range will increase or decrease in the future, we used Cliff's delta. This effect size measures the probability that a selected value from one group will be greater than a selected value from the other group, ranging from −1 to +1, with values closer to zero indicating no overlap between the two groups. Positive values indicated an increase in the elevational range of the species, while negative values indicated a decrease. The extent of elevational change was determined using the thresholds outlined in Romano et al. (2006): D < 0.147 is considered 'negligible,' D < 0.33 is considered 'small,' D < 0.474 is considered 'medium,' and anything beyond that is considered 'large’.

ResultsPerformance of species distribution models

Ensembles exhibited a greater mean performance compared to the mean of the individual models, as shown by a mean difference of TSS = +0.16, AUC = +0.13, and Boyce index = +0.04 (see Table S3). Furthermore, ensemble models demonstrated high accuracy in predicting species occurrence in the Andes. The TSS, AUC, and Boyce indices were positive and greater than 0.5, as shown in Table S4, which validated their usefulness in mapping the distribution of invasive species. The only exceptions were A. decurrens, which exhibited the poorest predictive performance (Boyce index = −0.5), and A. melanoxylon and S. madagascariensis (Boyce index close to 0, Table S4). As a result, these three species were excluded from further analysis.

How extensive is the range of conditions currently suitable for invasive species in the tropical Andes?

In the tropical Andes, H. coronarium has the largest suitable area among the invasive species, covering approximately 62.6% of the total area (300,000 out of 480,000 km2), followed by M. minutiflora (55.4%), P. aquilinum (45.9%), T. alata (46.7%), R. communis (42.9%), A. donax (25.4%), U. europaeus (18.4%), and G. monspessulana (15.1%). Most species can occupy multiple altitudinal zones, ranging from sea level to 4500 meters. However, the highest climatic suitability is projected to be confined to specific elevation boundaries (Fig. 2). Three of the eight species were most suited to environments above 2000 m, with U. europaeus and G. monspessulana being the two species restricted to high altitudes (2000–4000 meters). The other five species were better suited to environmental conditions at lower elevations, particularly below 2000 m.

Fig. 2.

Elevational ranges of invasive species in the Andes. Histograms display the projected amount of suitable area in intervals of 100 meters of elevation. Black vertical bars indicate elevation ranges where species have been recorded in the Andean region.

(0.35MB).

According to the models, the suitability of the environment for invasive species was determined by both climate and anthropogenic variables (see Supplementary Material 2). These variables included precipitation in the driest quarter (BIO17), which explained between 15% and 45% of species distribution; mean diurnal range of surface temperature (BIO2), which explained between 17% and 50%; soil temperature (SBIO2), which explained between 7% and 31%; mean growing season temperature (gst), which explained between 11% and 22%; livestock density, which explained between 10% and 31%; soil pH, which explained between 17% and 20%; and organic carbon content, which explained between 10% and 31%.

How will range change under different climate change scenarios?

The projected reduction in suitable areas is estimated to reach up to three times the amount of suitable area gained by invasive species by the mid-century, irrespective of the climate change scenario (Fig. 3a). On average, invasive species are projected to lose approximately 16,000 km2 of suitable area under the most benign scenario (SSP126). However, this area is projected to increase by a factor of 1.3 and 1.6, under the more carbon-intensive scenarios (SSP370 and 585, respectively). It is predicted that species such as U. europaeus, M. minutiflora, H. coronarium, and T. alata would lose the greatest suitable habitat (as indicated by the green bars in Fig. 3a), particularly in the southern Colombian Andes region, where the three mountain ranges converge (see Supplementary Material 2).

Fig. 3.

Projected range shifts of invasive species in the tropical mountains by mid-century. (a) Bar plots summarize the amount of suitable area projected to increase ('expansion' - purple) and decrease ('contraction' - green) under the three socio-economic pathways (SSPs) evaluated. (b) Box-and-whisker plots represent expected changes in median elevation for each species under each SSP. Each bar represents an average of the areas predicted by the three circulation models. (b) Box-and-whisker plots represent the expected changes in mean elevation of each species under each SSP by merging the projections of the three global circulation models.

(0.42MB).

In terms of the suitable area gained, species are expected to gain an average of approximately 6,100 km2 of suitable area under a benign climate change scenario and between 11,300 and 11,200 km2 under the other two carbon-intensive scenarios. It is anticipated that M. minutiflora, P. aquilinum, R. communes, and T. alata will experience the greatest increase in suitable areas, particularly within Colombian Andes. With regard to elevation, it should be noted that not all species gain suitable areas within the same elevation range (Fig. 3b). For example, H. coronarium and P. aquilinum will expand below their current median ranges, whereas A. donax and G. monspessulana will expand above their current median ranges. T. alata and U. europeaus were also projected to expand above their current mean elevations, but the magnitude of Cliff’s delta was negligible (See Table S4).

How much of the Páramo ecosystem is likely suitable for the establishment of invasive species?

The extent of suitable habitats for invasive species is estimated to cover up to 10% of the total area of the Páramo, which is approximately 3,000 km². However, this value varies depending on the species, the SSPs and the global circulation models (Fig. 4). H. coronarium and P. aquilinum have the smallest increase in suitable area, averaging less than 110 km2 in each SSP. The remaining species may have an increase in suitable area ranging from 500 km2 to 1200 km2 depending on the SSP (Fig. 4a). Global circulation models and SSPs are not entirely consistent in their predictions of areas with the highest future suitability (see Fig. 4b and c). All three models indicate that there will be an expansion of the suitable area in the western Andes under both low- and high-carbon scenarios (yellow and red areas in Fig. 4b and c). However, GFDL-ESM predicted a large expansion of suitable areas in the southern Andes between Ecuador and Colombia, which was only partially predicted by MPI-ESM and not fully predicted by IPSL-CM6A.

Fig. 4.

Projected shifts in the suitable area for invasive species in the Andean Páramos by mid-century. (A) The mean and standard deviation values computed for each species with predictions of the three algorithms and global circulation models are represented by points and bars. (B and C) Areas projected to gain suitability for one or more invasive species by the mid of century under SSP126 and SSP370.

(0.34MB).
Discussion

Tropical mountains provide a favorable environment for the potential establishment of invasive plant species. Currently, 18–60% of the region is estimated to have suitable conditions for at least one of the studied species. Our models indicated that suitability is higher in areas with higher dry season rainfall, stable temperatures throughout the year, specific soil properties (such as high pH and organic carbon), and/or dense livestock populations. However, as the climate changes, the spatial distribution of these conditions may change. By the mid-21st century, under the most benign scenario, invasive species could lose between 2,000 and 28,000 km2 of suitable environmental areas and gain between 1,000 and 10,000 km2 of more suitable areas. Although this result may seem counterintuitive, it is consistent with the expected effects of rising temperature on mountain plant species (Dirnböck et al., 2011). Specifically, rising temperatures are expected to affect the probability of species to establish new populations at lower elevations, particularly for cold-adapted species, resulting in upward shifts in their range (Alexander et al., 2018). Such upward shifts and range contractions have already been documented for U. europaeus in Colombia (Ángel-Vallejo et al., 2024), with a limited seed germination range of 0–15 °C (Udo et al., 2016) identified as the limiting factor, confining the species to the cooler highlands of the tropical Andes.

Our projections can be used to guide future conservation and management plans. Areas predicted to be highly suitable for invasive species could be targeted for conservation or management actions, such as population removal, blockage of invasion pathways, and active/passive restoration of ecosystems at a high risk of invasion (Roy et al., 2024). The implementation of such measures may prove pivotal in areas exhibiting a higher degree of suitability for more than one species (illustrated by the red cells in Figs. 4b and c), and notably in areas encompassing páramo ecosystems. The establishment and propagation of invasive species into páramos has the potential to pose a substantial risk to the region's water regulation services. Grasses such as A. donax and M. minutiflora have the potential to alter the soil structure and hinder the regeneration of native species (Zenni et al., 2022), which could ultimately impact the system's ability to retain water. Furthermore, the establishment of species, such as G. monspessulana and U. europaeus, combined with the projected longer droughts in the area (González‐Trujillo et al., 2024), may increase soil impoverishment and fire hazards (Galappaththi et al., 2023), ultimately compromising the paramos' capacity for water regulation and provisioning.

It is important to emphasize that given the life histories of the eight species studied, both modeling and management efforts should be species-specific and integrate multiple strategies. The suitability of the habitat for each species is determined not only by distinct abiotic factors (Supp. Mat. 1) but also by varying human activities that contribute to the establishment and spread of these species (Roy et al., 2024). Consequently, effective invasion control must go beyond managing and protecting the most suitable areas; it must also address key factors, such as invasion pathways. For several species, human-facilitated spread is driven by economic and cultural factors. For example, R. communis and M. minutiflora have been introduced for agricultural purposes, while T. Alata, A. donax, and H. coronarium have been spread for ornamental use. Thus, managing invasion pathways alongside the designation of priority protection areas is critical for effective control.

A more comprehensive approach to invasion control can be found in the restoration of native ecosystems. Restoration can function as a filter, impeding the establishment of invasive species, while facilitating the colonization of native species (Bakker and Wilson, 2004). Besides, local efforts to restore high Andean forests and páramos have shown promising results in controlling the spread of exotic grasses, gorses, and shrubs when combined with proactive management actions, such as the mechanical removal of species, prescribed burns, and artificial shading (Castiblanco-Álvarez, 2012; Gómez-Ruiz et al., 2013). The challenge ahead is to raise awareness of the pressing issue among stakeholders to increase the resources available for testing and scaling up the most feasible combinations of restoration and active management efforts.

Although our projections identify potential management priority areas until the mid-century, it is important to acknowledge the inherent limitations of our modelling approach. First, it should be noted that our models are conservative in nature, and that habitat suitability is predicted based on the assumption of no niche shift, which may not be applicable to all species (Atwater et al., 2017). For example, recent evidence indicates that U. europaeus has expanded its climate niche and colonized areas with climates distinct from its native range (Christina et al., 2020). This indicates that areas deemed suitable for invasion may be more expansive than our estimates suggest. Second, despite the careful selection of predictors, the analysis did not include important factors such as groundwater levels, spatial distribution of man-made corridors, and incidence of biotic interactions. The inclusion of data on these predictors could enhance the precision of our projections; however, they are currently unavailable at a global level or are not projected for future climate scenarios. For instance, the distribution of hydrophyte species such as A. donax is constrained to areas with shallow groundwater (Quinn and Holt, 2008), whereas species such as R. communis are primarily dispersed along man-made corridors such as roads and urban areas (Sandoval et al., 2022).

In the absence of more detailed knowledge, our projections provide the first comprehensive assessment of the suitability range of invasive species in the tropical Andes under different climatic scenarios. This information provides an overview of the magnitude of the problem and can help prioritize areas for implementing ecosystem-based management actions, such as restoration. In addition, these findings can inform future theoretical and empirical research on the environmental factors contributing to the spread of invasive species in the tropics, which has historically received less attention when studying biological invasions (Chiu et al., 2023).

CRediT authorship contribution statement

JDGT conceived the study. Species data collection and analyses were performed by JDGT. The first draft of the manuscript was written by JDGT, and all authors commented on subsequent versions of the manuscript. All authors read and approved the final manuscript.

Statements and declarations

The authors have no relevant financial or non-financial interests to disclose.

Acknowledgments

This research is funded by the Colombian Ministry of Science (Minciencias) through the project “DE LA MONTAÑA NO VENGO: MAPEO Y MODELIZACIÓN COMUNITARIAS DE LAS ESPECIES INVASORAS EN LA ALTA MONTAÑA TROPICAL - ESTANCIAS CON PROPÓSITO, PROGRAMA ONDAS, EN LOS MUNICIPIOS DE TUNJA, ARCABUCO, CIÉNEGA, DUITAMA, PESCA, SAMACÁ, SIACHOQUE Y GÜICÁN. SGI 3603”

Appendix A
Supplementary data

The following are Supplementary data to this article:

References
[Alexander et al., 2018]
J.M. Alexander, L. Chalmandrier, J. Lenoir, T.I. Burgess, F. Essl, S. Haider, C. Kueffer, K. McDougall, A. Milbau, M.A. Nuñez, A. Pauchard, W. Rabitsch, L.J. Rew, N.J. Sanders, L. Pellissier.
Lags in the response of mountain plant communities to climate change.
Global Change Biol., 24 (2018), pp. 563-579
[Anderson and Anderson, 2009]
S.A.J. Anderson, W.R. Anderson.
Predicting the elevated dead fine fuel moisture content in gorse (Ulex europaeus L.) shrub fuels.
Can. J. For. Res., 39 (2009), pp. 2355-2368
[Ángel-Vallejo et al., 2024]
M.C. Ángel-Vallejo, N. Aguirre-Acosta, G.T. Rodríguez-Rey, E.J. García-Marín, L.M. Álvarez-Mejía, C. Feuillet-Hurtado.
Distribution models in invasive plants with climatic niche expansion: a case study of Ulex europaeus L. in Colombian Andes.
Biol. Invasions, 26 (2024), pp. 1919-1930
[Araujo and New, 2007]
M. Araujo, M. New.
Ensemble forecasting of species distributions.
Trends in Ecol. Evol., 22 (2007), pp. 42-47
[Araújo et al., 2006]
M.B. Araújo, W. Thuiller, R.G. Pearson.
Climate warming and the decline of amphibians and reptiles in Europe.
J. Biogeogr., 33 (2006), pp. 1712-1728
[Atwater et al., 2017]
D.Z. Atwater, C. Ervine, J.N. Barney.
Climatic niche shifts are common in introduced plants.
Nat. Ecol. Evol., 2 (2017), pp. 34-43
[Bakker and Wilson, 2004]
J.D. Bakker, S.D. Wilson.
Using ecological restoration to constrain biological invasion.
J. Appl. Ecol., 41 (2004), pp. 1058-1064
[Barbet‐Massin et al., 2012]
M. Barbet‐Massin, F. Jiguet, C.H. Albert, W. Thuiller.
Selecting pseudo‐absences for species distribution models: how, where and how many?.
Methods Ecol. Evol., 3 (2012), pp. 327-338
[Berget et al., 2015]
C. Berget, E. Duran, D.B. Bray.
Participatory Restoration of Degraded Agricultural Areas Invaded by Bracken Fern (Pteridium aquilinum) and Conservation in the Chinantla Region, Oaxaca, Mexico.
Hum. Ecol., 43 (2015), pp. 547-558
[Boucher et al., 2020]
O. Boucher, J. Servonnat, A.L. Albright, O. Aumont, Y. Balkanski, V. Bastrikov, S. Bekki, R. Bonnet, S. Bony, L. Bopp, P. Braconnot, P. Brockmann, P. Cadule, A. Caubel, F. Cheruy, F. Codron, A. Cozic, D. Cugnet, F. D’Andrea, P. Davini, C. De Lavergne, S. Denvil, J. Deshayes, M. Devilliers, A. Ducharne, J. Dufresne, E. Dupont, C. Éthé, L. Fairhead, L. Falletti, S. Flavoni, M. Foujols, S. Gardoll, G. Gastineau, J. Ghattas, J. Grandpeix, B. Guenet, L. Guez, E. Guilyardi, M. Guimberteau, D. Hauglustaine, F. Hourdin, A. Idelkadi, S. Joussaume, M. Kageyama, M. Khodri, G. Krinner, N. Lebas, G. Levavasseur, C. Lévy, L. Li, F. Lott, T. Lurton, S. Luyssaert, G. Madec, J. Madeleine, F. Maignan, M. Marchand, O. Marti, L. Mellul, Y. Meurdesoif, J. Mignot, I. Musat, C. Ottlé, P. Peylin, Y. Planton, J. Polcher, C. Rio, N. Rochetin, C. Rousset, P. Sepulchre, A. Sima, D. Swingedouw, R. Thiéblemont, A.K. Traore, M. Vancoppenolle, J. Vial, J. Vialard, N. Viovy, N. Vuichard.
Presentation and evaluation of the IPSL‐CM6A‐LR climate model.
J. Adv. Model. Earth Syst., 12 (2020),
[Boyce et al., 2002]
M.S. Boyce, P.R. Vernier, S.E. Nielsen, F.K.A. Schmiegelow.
Evaluating resource selection functions.
Ecol. Modell., 157 (2002), pp. 281-300
[Brun et al., 2022]
P. Brun, N.E. Zimmermann, C. Hari, L. Pellissier, D.N. Karger.
Global climate-related predictors at kilometer resolution for the past and future.
Earth Syst. Sci. Data, 14 (2022), pp. 5573-5603
[Cárdenas López et al., 2017]
D. Cárdenas López, M.P. Baptiste, N. Castaño Arboleda, Instituto de Investigación de Recursos Biológicos Alexander von Humboldt.
Plantas exóticas con alto potencial de invasión en Colombia.
Instituto Humboldt Colombia, (2017),
[Castellanos et al., 2019]
A.A. Castellanos, J.W. Huntley, G. Voelker, A.M. Lawing.
Environmental filtering improves ecological niche models across multiple scales.
Methods Ecol. Evol., 10 (2019), pp. 481-492
[Castiblanco-Álvarez, 2012]
F. Castiblanco-Álvarez.
Control de pastos exóticos mediante sombreado artificial y reubicación de especies nativas como estrategias para la restauración ecológica del páramo andino (PNN Chingaza-Colombia).
UNiversidad Nacional de Colombia, (2012),
[Chen and Guestrin, 2016]
T. Chen, C. Guestrin.
XGBoost: A Scalable Tree Boosting System, in: Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining.
Presented at the KDD’ 16: The 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, ACM, (2016), pp. 785-794 http://dx.doi.org/10.1145/2939672.2939785
[Chiu et al., 2023]
J.H. Chiu, K.Y. Chong, S.K.Y. Lum, D.A. Wardle.
Trends in the direction of global plant invasion biology research over the past two decades.
Ecol. Evol., 13 (2023), pp. e9690
[Christina et al., 2020]
M. Christina, F. Limbada, A. Atlan.
Climatic niche shift of an invasive shrub (Ulex europaeus): a global scale comparison in native and introduced regions.
J. Plant Ecol., 13 (2020), pp. 42-50
[Costa et al., 2019]
R.O. Costa, A.F. Batisteli, E.L. Gaeta Espindola, D.M. Da Silva Matos.
Invasive Hedychium coronarium inhibits native seedling growth through belowground competition.
[D’Antonio et al., 2011]
C.M. D’Antonio, R.F. Hughes, J.T. Tunison.
Long-term impacts of invasive grasses and subsequent fire in seasonally dry Hawaiian woodlands.
Ecol. Appl., 21 (2011), pp. 1617-1628
[Da Re et al., 2023]
D. Da Re, E. Tordoni, J. Lenoir, J.J. Lembrechts, S.O. Vanwambeke, D. Rocchini, M. Bazzichetto.
USE it: Uniformly sampling pseudo‐absences within the environmental space for applications in habitat suitability models.
Methods Ecol. Evol., 14 (2023), pp. 2873-2887
[Dirnböck et al., 2011]
T. Dirnböck, F. Essl, W. Rabitsch.
Disproportional risk for habitat loss of high‐altitude endemic species under climate change.
Global Change Biol., 17 (2011), pp. 990-996
[Dormann et al., 2013]
C.F. Dormann, J.S. Elith Bacher, et al.
Collinearity: a review of methods to deal with it and a simulation study evaluating their performance.
[Dunne et al., 2020]
J.P. Dunne, L.W. Horowitz, A.J. Adcroft, P. Ginoux, I.M. Held, J.G. John, J.P. Krasting, S. Malyshev, V. Naik, F. Paulot, E. Shevliakova, C.A. Stock, N. Zadeh, V. Balaji, C. Blanton, K.A. Dunne, C. Dupuis, J. Durachta, R. Dussin, P.P.G. Gauthier, S.M. Griffies, H. Guo, R.W. Hallberg, M. Harrison, J. He, W. Hurlin, C. McHugh, R. Menzel, P.C.D. Milly, S. Nikonov, D.J. Paynter, J. Ploshay, A. Radhakrishnan, K. Rand, B.G. Reichl, T. Robinson, D.M. Schwarzkopf, L.T. Sentman, S. Underwood, H. Vahlenkamp, M. Winton, A.T. Wittenberg, B. Wyman, Y. Zeng, M. Zhao.
The GFDL Earth System Model Version 4.1 (GFDL‐ESM 4.1): overall coupled model description and simulation characteristics.
J. Adv. Model. Earth Syst, (2020), pp. 12
[Elith et al., 2008]
J. Elith, J.R. Leathwick, T. Hastie.
A working guide to boosted regression trees.
J. Animal Ecol., 77 (2008), pp. 802-813
[Fielding and Bell, 1997]
A.H. Fielding, J.F. Bell.
A review of methods for the assessment of prediction errors in conservation presence/absence models.
Envir. Conserv., 24 (1997), pp. 38-49
[Fuentes-Lillo et al., 2023]
E. Fuentes-Lillo, J.J. Lembrechts, A. Barros, V. Aschero, R.O. Bustamante, L.A. Cavieres, J. Clavel, I. Herrera, A. Jiménez, P. Tecco, P.E. Hulme, M.A. Núñez, R. Rozzi, R.A. García, D. Simberloff, I. Nijs, A. Pauchard.
Going up the Andes: patterns and drivers of non-native plant invasions across latitudinal and elevational gradients.
[Galappaththi et al., 2023]
H.S.S.D. Galappaththi, W.A.P.P. De Silva, A. Clavijo Mccormick.
A mini-review on the impact of common gorse in its introduced ranges.
Trop. Ecol., 64 (2023), pp. 1-25
[Gilbert et al., 2018]
M. Gilbert, G. Nicolas, G. Cinardi, T.P. Van Boeckel, S.O. Vanwambeke, G.R.W. Wint, T.P. Robinson.
Global distribution data for cattle, buffaloes, horses, sheep, goats, pigs, chickens and ducks in 2010.
[Girotto et al., 2021]
L. Girotto, A.C. Franco, C.V. Nunez, S.C.C. Oliveira, M.C. Scheffer De Souza, M.T. Fachin-Espinar, C.S. Ferreira.
Phytotoxicity and allelopathic potential of extracts from rhizomes and leaves of Arundo donax, an invasive grass in neotropical savannas.
Not. Bot. Horti. Agrobo., 49 (2021), pp. 12440
[Gómez-Ruiz et al., 2013]
P.A. Gómez-Ruiz, R. Lindig-Cisneros, O. Vargas-Ríos.
Facilitation among plants: a strategy for the ecological restoration of the high-andean forest (Bogotá, D.C.—Colombia).
Ecol. Eng., 57 (2013), pp. 267-275
[González‐Trujillo et al., 2024]
J.D. González‐Trujillo, D. Alagador, P. González‐del‐Pliego, M.B. Araújo.
Exposure of protected areas in Central America to extreme weather events.
Conserv. Biol., (2024),
[Gouvêa et al., 2024]
L. Gouvêa, E. Fragkopoulou, T. Legrand, E.A. Serrão, J. Assis.
Range map data of marine ecosystem structuring species under global climate change.
[Greenwell et al., 2019]
B. Greenwell, B. Boehmke, J. Cunningham, G. Developers, M.B. Greenwell.
Package ‘gbm.’.
R package version 2, (2019),
[Gutjahr et al., 2019]
O. Gutjahr, D. Putrasahan, K. Lohmann, J.H. Jungclaus, J.-S. Von Storch, N. Brüggemann, H. Haak, A. Stössel.
Max Planck Institute Earth System Model (MPI-ESM1.2) for the High-Resolution Model Intercomparison Project (HighResMIP).
Geosci. Model Dev., 12 (2019), pp. 3241-3281
[He et al., 2023]
X. He, A.D. Ziegler, P.R. Elsen, Y. Feng, J.C.A. Baker, S. Liang, J. Holden, D.V. Spracklen, Z. Zeng.
Accelerating global mountain forest loss threatens biodiversity hotspots.
One Earth, 6 (2023), pp. 303-315
[Heringer et al., 2021]
G. Heringer, E. Angulo, L. Ballesteros-Mejia, C. Capinha, F. Courchamp, C. Diagne, V.G. Duboscq-Carra, M.A. Nuñez, R.D. Zenni.
The economic costs of biological invasions in Central and South America: a first regional assessment.
[Herrera et al., 2016]
I. Herrera, E. Goncalves, A. Pauchard, R.O. Bustamante.
Manual de plantas invasoras de Sudamérica, Primera edición. ed.
IEB Chile, Instituto de Ecología y Biodiversidad, (2016),
[Hijmans, 2024]
R.J. Hijmans.
Terra: Spatial Data Analysis.
(2024),
[Hirzel et al., 2006]
A.H. Hirzel, G. Le Lay, V. Helfer, C. Randin, A. Guisan.
Evaluating the ability of habitat suitability models to predict species presences.
Ecol. Modell., 199 (2006), pp. 142-152
[Hofner et al., 2011]
B. Hofner, J. Müller, T. Hothorn.
Monotonicity-constrained species distribution models.
Ecology, 92 (2011), pp. 1895-1901
[Hofner et al., 2014]
B. Hofner, A. Mayr, N. Robinzonov, M. Schmid.
Model-based boosting in R: a hands-on tutorial using the R package mboost.
Comput. Stat., 29 (2014), pp. 3-35
[Iseli et al., 2023]
E. Iseli, C. Chisholm, J. Lenoir, S. Haider, T. Seipel, A. Barros, A.L. Hargreaves, P. Kardol, J.J. Lembrechts, K. McDougall, I. Rashid, S.B. Rumpf, J.R. Arévalo, L. Cavieres, C. Daehler, P.A. Dar, B. Endress, G. Jakobs, A. Jiménez, C. Küffer, M. Mihoc, A. Milbau, J.W. Morgan, B.J. Naylor, A. Pauchard, A. Ratier Backes, Z.A. Reshi, L.J. Rew, D. Righetti, J.M. Shannon, G. Valencia, N. Walsh, G.T. Wright, J.M. Alexander.
Rapid upwards spread of non-native plants in mountains across continents.
Nat. Ecol. Evol., 7 (2023), pp. 405-413
[Krause-Jensen et al., 2020]
D. Krause-Jensen, P. Archambault, J. Assis, I. Bartsch, K. Bischof, K. Filbee-Dexter, K.H. Dunton, O. Maximova, S.B. Ragnarsdóttir, M.K. Sejr, U. Simakova, V. Spiridonov, S. Wegeberg, M.H.S. Winding, C.M. Duarte.
Imprint of climate change on Pan-Arctic marine vegetation.
[Kumschick and Jansen, 2023]
S. Kumschick, C. Jansen.
Evidence-Based Impact Assessment for Naturalized and Invasive Australian Acacia Species.
[Lambert et al., 2010]
A.M. Lambert, T.L. Dudley, K. Saltonstall.
Ecology and Impacts of the Large-Statured Invasive Grasses Arundo donax and Phragmites australis in North America.
Invasive Plant Sci. Manag., 3 (2010), pp. 489-494
[Le Maitre et al., 2011]
D.C. Le Maitre, M. Gaertner, E. Marchante, E. Ens, P.M. Holmes, A. Pauchard, P.J. O’Farrell, A.M. Rogers, R. Blanchard, J. Blignaut, D.M. Richardson.
Impacts of invasive Australian acacias: implications for management and restoration.
Divers. Distrib., 17 (2011), pp. 1015-1029
[Lembrechts et al., 2022]
J.J. Lembrechts, J. van den Hoogen, J. Aalto, M.B. Ashcroft, P. De Frenne, J. Kemppinen, M. Kopeckỳ, M. Luoto, I.M. Maclean, T.W. Crowther, et al.
Global maps of soil temperature.
Global Change Biol., 28 (2022), pp. 3110-3144
[Leroy et al., 2018]
B. Leroy, R. Delsol, B. Hugueny, C.N. Meynard, C. Barhoumi, M. Barbet‐Massin, C. Bellard.
Without quality presence–absence data, discrimination metrics such as TSS can be misleading measures of model performance.
J. Biogeogr., 45 (2018), pp. 1994-2002
[Liu et al., 2016]
C. Liu, G. Newell, M. White.
On the selection of thresholds for predicting species occurrence with presence‐only data.
Ecol. Evol., 6 (2016), pp. 337-348
[Liu et al., 2023]
Y. Liu, X. Xu, D. Dimitrov, L. Pellissier, M.K. Borregaard, N. Shrestha, X. Su, A. Luo, N.E. Zimmermann, C. Rahbek, Z. Wang.
An updated floristic map of the world.
Nat. Commun., 14 (2023), pp. 2990
[Nogueira et al., 2019]
C.B. Nogueira, E. Menéndez, M.H. Ramírez-Bahena, E. Velázquez, Á. Peix, P.F. Mateos, M.R. Scotti.
The N-fixing legume Periandra mediterranea constrains the invasion of an exotic grass (Melinis minutiflora P. Beauv) by altering soil N cycling.
Sci. Rep., 9 (2019), pp. 11033
[Pauchard et al., 2008]
A. Pauchard, R.A. García, E. Peña, C. González, L.A. Cavieres, R.O. Bustamante.
Positive feedbacks between plant invasions and fire regimes: Teline monspessulana (L.) K. Koch (Fabaceae) in central Chile.
Biol Invasions, 10 (2008), pp. 547-553
[Peyre, 2022]
G. Peyre.
What does the future hold for páramo plants? A modelling approach.
Front. Ecol. Evol., 10 (2022),
[Poggio et al., 2021]
L. Poggio, L.M. De Sousa, N.H. Batjes, G.B.M. Heuvelink, B. Kempen, E. Ribeiro, D. Rossiter.
SoilGrids 2.0: producing soil information for the globe with quantified spatial uncertainty.
SOIL, 7 (2021), pp. 217-240
[Quijano-Abril et al., 2021]
M.A. Quijano-Abril, M.D.L.Á. Castaño López, D. Marín-Henao, D.A. Sánchez-Gómez, J.M. Rojas-Villa, J. Sierra-Escobar.
Rasgos funcionales de la especie invasora Thunbergia alata (Acanthaceae) y su importancia en la adaptación a bosques andinos.
[Quinn and Holt, 2008]
L.D. Quinn, J.S. Holt.
Ecological correlates of invasion by Arundo donax in three southern California riparian habitats.
Biol. Invasions, 10 (2008), pp. 591-601
[Riahi et al., 2017]
K. Riahi, D.P. Van Vuuren, E. Kriegler, J. Edmonds, B.C. O’Neill, S. Fujimori, N. Bauer, K. Calvin, R. Dellink, O. Fricko, W. Lutz, A. Popp, J.C. Cuaresma, S. Kc, M. Leimbach, L. Jiang, T. Kram, S. Rao, J. Emmerling, K. Ebi, T. Hasegawa, P. Havlik, F. Humpenöder, L.A. Da Silva, S. Smith, E. Stehfest, V. Bosetti, J. Eom, D. Gernaat, T. Masui, J. Rogelj, J. Strefler, L. Drouet, V. Krey, G. Luderer, M. Harmsen, K. Takahashi, L. Baumstark, J.C. Doelman, M. Kainuma, Z. Klimont, G. Marangoni, H. Lotze-Campen, M. Obersteiner, A. Tabeau, M. Tavoni.
The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: an overview.
Global Environ. Change, 42 (2017), pp. 153-168
[Rojas‐Soto et al., 2024]
O. Rojas‐Soto, J.S. Forero‐Rodríguez, A. Galindo‐Cruz, C. Mota‐Vargas, K.D. Parra‐Henao, A. Peña‐Peniche, J. Piña‐Torres, K. Rojas‐Herrera, J.D. Sánchez‐Rodríguez, F.A. Toro‐Cardona, C.D. Trinidad‐Domínguez.
Calibration areas in ecological niche and species distribution modelling: unravelling approaches and concepts.
J. Biogeogr., 51 (2024), pp. 1416-1428
[Romano et al., 2006]
J. Romano, J.D. Kromrey, J. Coraggio, J. Skowronek.
Appropriate statistics for ordinal level data: Should we really be using t-test and Cohen’sd for evaluating group differences on the NSSE and other surveys.
in: Annual Meeting of the Florida Association of Institutional Research, (2006), pp. 34
[Roy et al., 2024]
H.E. Roy, A. Pauchard, P. Stoett, T. Renard Truong.
IPBES invasive alien species assessment: full report.
[Ruiz et al., 2008]
D. Ruiz, H.A. Moreno, M.E. Gutiérrez, P.A. Zapata.
Changing climate and endangered high mountain ecosystems in Colombia.
Sci. Total Environ., 398 (2008), pp. 122-132
[Sandoval et al., 2022]
L.M. Sandoval, O. Vargas Ríos, J.D. González-Trujillo.
Invasion of plants in towns and roadsides around the Iguaque massif (Boyacá, Colombia).
[Silva and Fabricante, 2022]
F.O.D. Silva, J.R. Fabricante.
Impacts of the biological invasion by Ricinus communis L. on the native biota of the Atlantic Forest, Aracaju, Sergipe State, Brazil.
[Udo et al., 2016]
N. Udo, M. Tarayre, A. Atlan.
Evolution of germination strategy in the invasive species Ulex europaeus.
[Valdez-Ramírez et al., 2020]
C. Valdez-Ramírez, S.I. Levy-Tacher, N.S. León-Martínez, D.A. Navarrete-Gutiérrez, Á.I. Ortiz-Ceballos.
Cambios químicos y biológicos del suelo provocados por Pteridium aquilinum (L.) Kuhn en áreas de influencia de la reserva de la biosfera de Calakmul, Campeche.
Terra, 38 (2020), pp. 289-300
[Valencia et al., 2013]
J. Valencia, L. Lassaletta, E. Velázquez, J.M. Nicolau, A. Gómez‐Sal.
Factors Controlling Compositional Changes in a Northern Andean Páramo (La Rusia, Colombia).
[Varela et al., 2014]
S. Varela, R.P. Anderson, R. García‐Valdés, F. Fernández-González.
Environmental filters reduce the effects of sampling bias and improve predictions of ecological niche models.
Ecography, 37 (2014), pp. 1084-1091
[Weigend, 2002]
M. Weigend.
Observations on the Biogeography of the Amotape-Huancabamba Zone in Northern Peru.
[Wijayabandara et al., 2021]
K. Wijayabandara, S. Campbell, J. Vitelli, A. Shabbir, S. Adkins.
Review of the biology, distribution, and management of the invasive fireweed (Senecio madagascariensis Poir).
[Williams, 2011]
J.N. Williams.
Human Population and the Hotspots Revisited: A 2010 Assessment.
Eds.), Biodiversity Hotspots., pp. 61-81 http://dx.doi.org/10.1007/978-3-642-20992-5_4
[Zenni et al., 2022]
R.D. Zenni, I. Herrera, M. De Sá Dechoum, S.R. Ziller, A.C.L. De Matos, C.I. Núñez, M.A. Núñez, A. Pauchard.
Plant Invasions in South America.
Global Plant Invasions, pp. 187-208 http://dx.doi.org/10.1007/978-3-030-89684-3_9
Copyright © 2024. Associação Brasileira de Ciência Ecológica e Conservação
Download PDF
Perspectives in Ecology and Conservation
Article options
Tools
Supplemental materials