Translator Disclaimer
30 November 2020 Predicting the Current and Future Distribution of the Invasive Weed Ageratina adenophora in the Chitwan–Annapurna Landscape, Nepal
Author Affiliations +

With increasing globalization, trade, and human movement, the rate of alien species introduction has increased all around the globe. In addition, climate change is thought to exacerbate the situation by allowing range expansion of invasive species into new areas. Predicting the distribution of invasive species under conditions of climate change is important for identifying susceptible areas of invasion and developing strategies for limiting their expansion. We used Maxent modeling to predict the distribution of one of the world's most aggressive invasive weeds, Ageratina adenophora (Sprengel) R. King and H. Robinson, in the Chitwan–Annapurna Landscape (CHAL) of Nepal under current conditions and 3 future climate change trajectories based on 3 representative concentration pathways (RCPs 2.6, 4.5, and 8.5) in 2 different time periods (2050 and 2070) using species occurrence data, and bioclimatic and topographic variables. Minimum temperature in the coldest month was the most important variable affecting the distribution of A. adenophora. About 38% (12,215 km2) of the CHAL area is climatically suitable for A. adenophora, with the Middle Mountain physiographic region being the most suitable one. A predicted increase in current suitable areas ranges from 1 to 2% under future climate scenarios (RCP 2.6 and RCP 8.5). All protected areas and 3 physiographic regions (Siwaliks, High Mountain, High Himalaya) are likely to gain climatically suitable areas in future climate scenarios. The upper elevational distribution limit of the weed is expected to expand by 31–48 m in future climate scenarios, suggesting that the weed will colonize additional areas at higher elevations in the future. In conclusion, our results showed that a vast area of CHAL is climatically suitable for A. adenophora. Expected further range expansion and upslope migration in the future make it essential to initiate effective management measures to prevent further negative impacts of this invasive plant.


Biological invasions, a major driver of global environmental changes, are posing serious threats to global biodiversity and ecosystem functioning (IPBES 2019). Climate change is likely to further amplify the risks of biological invasions (Walther et al 2009; Bradley, Wilcove, et al 2010). Biological invasions and climate change act synergistically, and this synergistic relation between the 2 parameters of global change has been identified as a major threat to biodiversity (Dukes and Mooney 1999; Walther et al 2009; Mainka and Howard 2010). The profound negative impacts of invasive plant species on the diversity of native species, soil dynamics, and ecosystem processes, which cause ecological and economic losses, are well known (Marbuah et al 2014; Villa and Hulme 2017; Castro-Diez et al 2019). These impacts are likely to be further exacerbated by climate change, enhancing traits that promote invasiveness and creating a more hospitable climate for invasive species to cross geographic barriers, thereby facilitating range expansion in new areas (Dukes and Mooney 1999; Stacowicz et al 2002; Walther et al 2009; Bradley, Blumenthal, et al 2010; Bellard et al 2013). In comparison to native plant species, parameters of global change, such as increased temperature and CO2 enrichment, enhance the performance of invasive species, imposing a threat of further spread (Liu et al 2017). Therefore, predicting the distribution of invasive weeds under climate change scenarios and identifying the areas potentially at risk are urgent needs for effective management planning to minimize ecological and economic impacts.

A first step to identify the risk of invasions is to use ecological niche models (ENMs) to predict suitable ecological niches for a species across a landscape. These relate documented presence records of the focal species with the environmental or spatial characteristics of the potential sites (Elith and Leathwick 2009; Franklin 2009). The niche concept is central to ENMs and is based on Hutchinson's (1957) concept of fundamental and realized niches (Araujo and Guisan 2006). It is highly likely that an invasive species at an early stage of invasion occupies only a small fraction (ie realized niche) of the fundamental niche in the introduced range, and there is always a risk of invasion in the unoccupied part of the fundamental niche (Soberon and Nakamura 2009). However, controversies persist over which facets of the niche are projected by ENMs (Araujo and Guisan 2006; McInerny and Etienne 2012). ENMs have been gaining popularity and are widely used by ecologists in invasive species risk assessments (Qin et al 2016; Suarez-Mota et al 2016; Wan et al 2017; Shrestha et al 2018; Thapa et al 2018). Among different ENMs, Maxent is one of the most popular species distribution modeling tools. This model uses presence-only records and has been commonly used in building habitat suitability maps for invasive species (Phillips et al 2006; Merow et al 2013; West et al 2016; Lamsal et al 2018).

Among 124 countries, Nepal has the third highest threat to agriculture sectors from invasive species (Paini et al 2016). To date, 179 species of flowering plants are naturalized, and, among them, 26 species are reported to be invasive in Nepal (Shrestha 2019; Shrestha, Budha, et al 2019). The crofton weed, Ageratina adenophora (Sprengel) R. King and H. Robinson (Asteraceae), is one of the most noxious invasive weeds in many parts of Asia, Oceania, and Africa. It has had serious ecological impacts on native biodiversity and caused enormous economic losses (Poudel et al 2019). It is ranked as the most problematic invasive weed in Nepal (Tiwari et al 2005).

Though many studies have been carried out to investigate the potential distribution of A. adenophora on broader spatial scales in Nepal (Shrestha and Shrestha 2019), China (Wang and Wang 2006; Wang et al 2017), the Himalayas (Lamsal et al 2018; Thapa et al 2018), and South Africa (Tererai and Wood 2014), there is a lack of such studies on smaller scales where management strategies are implemented. One of the most important landscapes in Nepal is the Chitwan–Annapurna Landscape (CHAL), located in central Nepal. CHAL harbors rich biodiversity due to its wide elevation gradient (200–8091 m above sea level [masl]), diverse topography, and climatic variations (subtropical to alpine) (WWF 2013; MFSC 2016). Invasion by alien species has already been recognized as a major threat to biodiversity in CHAL (WWF 2013), and A. adenophora is the invasive weed most prioritized by the local communities for management in natural ecosystems due to its negative impacts on biodiversity and livelihoods (Shrestha, Shrestha, et al 2019). Therefore, there is an urgent need to recognize potential areas of distribution of A. adenophora in CHAL under the current climate and identify areas at risk of being invaded by this weed under future climate scenarios.

In this study, we used the Maxent modeling tool to predict the current and future potential distribution of A. adenophora in CHAL using occurrence records from different sources. The objective was to prepare habitat suitability maps for the weed under current climatic conditions and future climate scenarios (RCP 2.6, RCP 4.5, and RCP 8.5 in the years 2050 and 2070) to identify the key environmental factors influencing its distribution and areas at risk of invasion. Information on its potential distribution will be very useful for the scientific community and managers in developing future monitoring and management strategies to prevent further expansion of the weed in this landscape.


Study area

CHAL is located in central Nepal and covers 19 districts. This landscape has a wide elevation gradient, ranging from 200 to 8091 masl, and covers an area of 32,057 km2 (WWF 2013). It spans 4 physiographic regions, namely, Siwalik, Middle Mountain, High Mountain, and High Himalaya. They have diverse climatic conditions, from subtropical in Siwalik to alpine in the High Himalaya, and a cold and dry climate in Trans-Himalayan regions. This geographic and climatic diversity shapes the habitat and environmental conditions for CHAL's rich biodiversity, which includes more than 104 species of mammals (Bhuju et al 2007), 500 species of birds (Baral and Inskipp 2005; Bhuju et al 2007), and 3430 species of plants, with high levels of endemism and genetic diversity (BPP 1995). Forests and grasslands are the main natural ecosystems, occupying 35.5% and 8.6% of the landscape, respectively, whereas 21.1% of the area is under agriculture (WWF 2013). The region has a population of 4.5 million people (CBS 2013). The average minimum and maximum temperatures are 5°C and 40°C, and the average annual rainfall ranges from 165 to 5244 mm (MFSC 2016). The landscape includes portions of 4 globally recognized ecoregions and comprises 3 national parks (Chitwan, Parsa, and Langtang) and 2 conservation areas (Annapurna and Manaslu) (MFSC 2016). Annapurna Conservation Area and Chitwan National Park are among the sites with a high number of visiting tourists (DNPWC 2018). Most of the lower and mid-hill forests of this region are at risk of fragmentation and conversion to other vegetation types due to climate change (Thapa et al 2015). In comparison to the eastern and western regions of Nepal, central Nepal, where CHAL is located, hosts higher numbers of naturalized plant species (Bhattarai et al 2014). The combination of diverse natural environment along with anthropogenic disturbances has made the region vulnerable to biological invasions (WWF 2013).

Study species

A. adenophora is native to Mexico but is established in 40 countries outside its native range (Poudel et al 2019). In Nepal, 26 invasive plant species are reported, 20 of which have been documented in the CHAL region, where A. adenophora is the most problematic weed in natural ecosystems, and its management is highly prioritized by local people (Tiwari et al 2005; Shrestha 2019; Shrestha, Shrestha, et al 2019). It was first reported in 1958, having been accidentally introduced from the eastern border of India to Nepal (Tiwari et al 2005). It covers a wide elevation gradient of 400–3280 masl in Nepal (Siwakoti et al 2016). The ability of A. adenophora to occupy a wide range of climatic habitats and spread rapidly can be attributed to its phenotypic plasticity, allelopathy, and ability to alter the soil microbial community to favor its further invasion (Poudel et al 2019). It has reduced forage supply, displaced native plant species, causing loss of biodiversity, and prevented forest regeneration in CHAL (WWF 2013; Shrestha, Shrestha, et al 2019).

Species occurrence data

In total, 686 occurrence points were collected from different sources. We noted 245 occurrence points from secondary sources (Siwakoti et al 2016; Shrestha and Shrestha 2019). The remaining 441 occurrence points were collected by the first author during field visits in 2016–2019. Road networks are the major conduit for dispersal of A. adenophora (Dong et al 2008); therefore, occurrence data were mainly collected through field surveys along roadsides and trekking routes (Figure 1). This sampling bias was addressed by spatial filtering of the data. Duplicate records of occurrence points were deleted and spatially thinned using the spThin package (Aiello-Lammens et al 2015) in the R software (version 3.4.4) (R Core Team 2017), so that only 1 single location occurred in each 1 km2 grid cell. Spatial filtering makes it possible to reduce overfitting to sampling bias in ENMs (Boria et al 2014). In total, 403 occurrence records obtained after filtering were used to build the models (Figure 1).


Study area, with elevation zone, major rivers, road networks, and occurrence locations of Ageratina adenophora in the Chitwan–Annapurna Landscape, Nepal. (Source: Survey Department, Government of Nepal)


Environmental and bioclimatic variables

In December 2018, we downloaded 19 grid-based bioclimatic variables that represent annual trends, seasonality, and extreme climatic conditions from the WorldClim database (version 1.4) (; Hijmans et al 2005) at a spatial resolution of 30 arc-seconds (∼1 km2) ( Appendix S1 (mred-40-02-15_s01.pdf), Supplemental material, Elevation was obtained from the Shuttle Radar Topographic Mission (SRTM) at 90 m spatial resolution. This was then resampled into 30 arc-second spatial resolution by using the nearest neighbor resampling technique in ArcGIS (version 10.3). Slope and aspect rasters of the study area were derived from the elevation data.

To predict climatically suitable areas in future climate scenarios, we chose projections from the Community Climate System Model (CCSM4) under b1 emission scenarios, which are based on the fifth phase of the Coupled Model Intercomparison Project5 (CMIP5) (Gent et al 2011). We selected 3 greenhouse gas (GHG) emission scenarios, also known as Representative Carbon Pathways (RCP 2.6, RCP 4.5, and RCP 8.5), for 2 different time periods (2050 and 2070) as adopted by the Intergovernmental Panel on Climate Change in its Fifth Assessment Report (AR5) (IPCC 2013). RCP 2.6, RCP 4.5, and RCP 8.5 represent the lowest, medium, and highest emission scenarios, corresponding to a 1.0°C, 1.4–1.8°C, and 2.0–3.7°C projected increase in global mean surface temperature, respectively (van Vuuren et al 2011; IPCC 2013).

The datasets were extracted for the study area (CHAL) using the Spatial Analyst Tool and the Extraction Tool in ArcGIS (version 10.3). These datasets were converted from raster format to ASCII files in ArcGIS. We repeated the procedure to prepare the predicted distribution maps for the 2 future climate scenarios (2050 and 2070).

Model building

Maxent (version 3.3.3) software was used to build the ENMs (Phillips and Dudik 2008). Maxent is one of the most commonly used habitat suitability modeling techniques. It uses presence-only data and is widely used for invasive species (Phillips et al 2006; Shrestha et al 2018; Maharjan et al 2019). Maxent works well with incomplete or limited data, so it can provide robust estimates of potentially suitable habitats for invasive species at small spatial scales (Jarnevich et al 2006; Jarnevich and Reynolds 2011; West et al 2016).

To reduce multicollinearity among predictor variables (19 bioclimatic and 3 topographic variables) and overfitting of the model, pairwise correlation analyses were performed in R, and highly correlated variables with a Pearson's correlation coefficient ≥0.8 were removed ( Appendix S2 (mred-40-02-15_s01.pdf), Supplemental material, (Merow et al 2013). While selecting one variable from each pair of highly correlated variables (r2 > 0.8), special attention was paid to existing biological and ecological insights into the species. Ultimately, 7 bioclimatic variables—isothermality (Bio 3), minimum temperature of the coldest month (Bio 6), temperature annual range (Bio 7), precipitation of the driest month (Bio 14), precipitation seasonality (Bio 15), precipitation of the warmest quarter (Bio 18), and precipitation of the coldest quarter (Bio 19)—and 2 topographic variables—aspect and slope—were used as predictors to build the habitat suitability model.

The Maxent model used 75% of the data for training and the remaining 25% for testing. We used a logistic format because it improves model calibration by estimating the probability of a species being present depending on environmental variables (Phillips and Dudik 2008). A convergence threshold of 10–5, a maximum iteration value of 5000, 15 replications with a replicated run type subsample, and 10,000 random background points were used to build the model. As a threshold rule, we chose tenth percentile training presence on the basis of the area under the curve (AUC) and true skill statistic (TSS). Tenth percentile training presence omits the 10% of localities or training presence records with the lowest predicted values (Radosavljevic and Anderson 2014) and is highly conservative in estimating species tolerance with respect to each climatic variable (Svenning et al 2008). The remaining parameters were kept at their default values.

We imported the Maxent output, which is continuous data with values ranging from 0 to 1, into ArcGIS (version 10.3) and classified the map using the Reclassify Tool into 2 classes, suitable habitat and unsuitable habitat, on the basis of a tenth percentile training presence logistic threshold. In this way, a binary habitat suitability map was created for the current and all future climate scenarios. We also calculated the climatically suitable area for present and future climate scenarios, as well as changes in suitable areas in terms of gain, loss, and stable areas in the future under all scenarios in ArcGIS. Data for physiographic regions and protected areas were then clipped to projected maps with suitable areas for current and future climate scenarios to calculate the suitable areas in these physiographic regions and protected areas. Changes in the upper and lower elevational distribution range under future climate scenarios in comparison to current were quantified using the Extraction Tool and digital elevation model (DEM) raster for all maps, current and future, in ArcGIS.

Model evaluation

Threshold-independent (area under the receiver operating characteristic [ROC] curve [AUC]) and threshold-dependent (TSS) measures of model accuracy were used to evaluate model performance (Fielding and Bell 1997; Allouche et al 2006; Franklin 2009). AUC values range from 0–1.0, with 0.5–0.7 considered low, 0.7–0.9 moderate, and >0.9 high (Swets 1988; Manel et al 2001). The TSS value ranges from –1 to +1, where +1 indicates a perfect agreement, and 0 or less indicates a performance no better than random (Allouche et al 2006). Marginal response curves were used to visually investigate the relationship between environmental variables (predictors) and the predicted index of habitat suitability of A. adenophora. The relative contribution of different predictor variables to the Maxent model was assessed by the variable percentage contribution and jackknife procedures (Elith et al 2011). The jackknife test of variable importance helps to identify those variables with important individual effects (Elith et al 2011). Two jackknife tests were taken into account. The jackknife test of regularized training gain shows the training gain of a variable when used in isolation and the training gain of a variable when omitted, and it compares these values to the training gain of all variables. Similarly, the jackknife test of AUC based on the AUC of test data shows the predictive performance of the variable when used in isolation and the predictive performance of the variable when omitted, and it compares these values with the AUC value when all variables are used (Phillips 2017).


Model performance and variable contribution

The current model for A. adenophora performed better than random, with a mean training AUC value of 0.85, mean test AUC value of 0.80, and a TSS value of 0.52 (Table 1). Responses of each predictor variable are shown in the response curve ( Appendix S3 (mred-40-02-15_s01.pdf), Supplemental material, Out of 9 predictor variables used for model building, minimum temperature of the coldest month (Bio 6) contributed the most (48.7%), followed by precipitation of the warmest quarter (Bio 18). Aspect had the lowest contribution (2.9%) (Table 2). Maxent's jackknife test of variable importance also showed that minimum temperature of the coldest month (Bio 6) had the highest training gain and AUC, followed by precipitation of the warmest quarter (Bio 18) when used in isolation (Figure 2). The response curve of the variable minimum temperature of the coldest month (Bio 6) showed that the probability that the weed will occur below 1°C was the lowest, increasing with increasing minimum temperature ( Appendix S3 (mred-40-02-15_s01.pdf), Supplemental material, In addition, the minimum temperature of the coldest month of about 96% of the occurrence points used in model building was above 0°C ( Appendix S4 (mred-40-02-15_s01.pdf), Supplemental material, Similarly, the response curve of precipitation of the warmest quarter (Bio 18) indicated that the probability the weed will occur increased with increasing precipitation in the warmest quarter above 500 mm ( Appendix S3 (mred-40-02-15_s01.pdf), Supplemental material,


Model evaluation matrices.



Relative contribution of the environmental variables to the Maxent model built for current climatic conditions.



Results of jackknife test of relative importance of predictor variables for A. adenophora for the current distribution. (A) Jackknife of regularized training gain. (B) Jackknife of AUC. Predictors used: Bio 3, isothermality; Bio 6, minimum temperature of coldest month; Bio 7, temperature annual range; Bio 14, precipitation of driest month; Bio 15, precipitation seasonality; Bio 18, precipitation of warmest quarter; Bio 19, precipitation of coldest quarter; slope; aspect.


Current potential distribution

Currently, 38% (12,215 km2) of the total area of CHAL is suitable for A. adenophora (Figure 3; Table 3). All districts of CHAL within the elevational range of 119–2824 masl had climatically suitable areas for this weed. The most suitable areas for A. adenophora were found to be in the Middle Mountain physiographic region (75%), followed by Siwalik (37%), High Mountain (29%), and High Himalaya (0.1%) (Table 4). High Himalaya is hardly suitable for the weed, with the fewest suitable areas found in the districts of Manang and Mustang. Nonetheless, all districts of the CHAL region were found to have climatically suitable areas. Because the Middle Mountain and High Mountain regions had the most suitable areas for the weed, we tried to observe the change in climatically suitable areas in 3 protected areas situated in these physiographic regions. Among the 3 protected areas—Annapurna Conservation Area, Langtang National Park, and Manaslu Conservation Area—Langtang National Park was predicted to have the highest percentage of area suitable for this weed (Figure 4).


Predicted suitable area for A. adenophora in the Chitwan–Annapurna Landscape, Nepal, under (A) current scenario; (B) RCP 2.6 for the year 2050; (C) RCP 2.6 for the year 2070; (D) RCP 4.5 for the year 2050; (E) RCP 4.5 for the year 2070; (F) RCP 8.5 for the year 2050; (G) RCP 8.5 for the year 2070. For future climate scenarios, likely stable, gain, and loss in areas are shown in bar graphs denoted by green, red, and yellow, respectively.



Predicted climatically suitable area for A. adenophora under current and future climate scenarios in the Chitwan–Annapurna Landscape, Nepal.



Predicted suitable area for A. adenophora in different physiographic regions of the Chitwan–Annapurna Landscape, Nepal.



Change in suitable areas of A. adenophora in different protected areas of the Chitwan–Annapurna Landscape, Nepal, under current and future climate scenarios.


Future invasion risk and change in habitat suitability

The predicted climatically suitable areas for A. adenophora would increase under RCP 2.6 for the year 2070 and RCP 4.5 for both the years 2050 and 2070 (Figure 3; Table 3). The highest increase (2%) in area of suitable habitat was predicted for the year 2070 under RCP 2.6 and 4.5. This gain in suitable areas was more prominent in districts like Lamjung, Gorkha, Dhading, Makwanpur, Chitwan, and Tanahun. However, an increase in radiative force (from RCP 4.5 to 8.5) would decrease climatically suitable areas for the weed in both the years 2050 and 2070 (Figure 3; Table 3).

Though the suitable area was predicted to decrease in extreme climate scenarios (RCP 8.5), the upper elevational distribution limit would expand by 31 m and 42 m for the years 2050 and 2070, respectively (Figure 5). Though a small decrease (24 m) in upper elevation limit was predicted under RCP 2.6 in 2050, the highest increase of 48 m was expected under RCP 2.6 for the year 2070 in comparison to current climatic conditions. However, under the medium emission scenario of RCP 4.5, the model predicted a contraction in the upper elevational limit for both years 2050 and 2070. In contrast, the lower elevation limit of A. adenophora would either remain stable or contract in future climate scenarios (Figure 5).


Change in upper and lower elevation range of A. adenophora in future climate scenarios compared to current conditions. (A) Year 2050. (B) Year 2070.


With climate change, all physiographic regions except the Middle Mountain region were expected to gain climatically suitable areas. Though a minimal loss in suitable areas was predicted for Middle Mountain, this region will still contain the most suitable areas for the weed, followed by High Mountain, Siwalik, and High Himalaya, in all future climate scenarios (Table 4). In 4 of the future climate scenarios—RCP 2.6 in 2050 and 2070, and RCP 8.5 in 2050 and 2070—Siwalik will gain climatically suitable areas. Except under RCP 8.5 for the year 2050, High Mountain will also gain suitable areas in all future climate scenarios. Furthermore, under RCP 4.5 and 8.5 for the year 2070, the climatically suitable area is also predicted to increase in High Himalaya. Among all physiographic regions, the percentage gain in suitable area was highest (2.4%) in Siwalik under RCP 8.5 for the year 2050 (Table 4).

Like under current climatic conditions, in future climate scenarios, Langtang National Park will have more suitable areas than the other 2 protected areas (Figure 4). Under RCP 2.6 and 4.5 for both 2050 and 2070, climatically suitable areas are predicted to increase in Annapurna Conservation Area, whereas under extreme climate scenarios, it will lose some suitable areas. For Langtang National Park, climatically suitable areas will increase under RCP 2.6 and 8.5 for both years, but not under RCP 4.5 (likewise for both years). For Manaslu Conservation Area, a remarkable gain in suitable areas is predicted only under RCP 8.5 for the year 2070.


This study is the first to predict current and future suitable habitat for A. adenophora in CHAL, Nepal. The model evaluation parameters (AUC and TSS) obtained for our models both lie within a range that confirms the robustness of the models (Table 1): AUC values above 0.8 and TSS values closer to 1 are considered to be acceptable (Swets 1988; Manel 2001; Allouche et al 2006). Though our study did not include a model transferability assessment, this would provide valuable information for model validation as well as model selection (Wenger and Olden 2012).

Climatic factors, such as temperature and precipitation, play a pivotal role in determining the pathways and success of plant invasions (Kathiresan and Gualbert 2016; Wang et al 2017). Our predictions suggest that the minimum temperature in the coldest month is the most influential factor for the distribution of A. adenophora. This is in line with the findings of other researchers (Wang and Wang 2006; Zhu et al 2007; Wang et al 2017; Lamsal et al 2018; Thapa et al 2018; Datta et al 2019). According to our model, the minimum temperature in winter (Bio 6) alone contributes about 49% to the habitat suitability model (Table 2), and the probability of occurrence of the weed decreased to almost 0 below the minimum winter temperature of 1°C ( Appendix S3 (mred-40-02-15_s01.pdf), Supplemental material, This indicates that the occurrence of this weed in high mountains is constrained by the lower minimum winter temperature, which is supported by observations in India, where the low temperature in winter limits the uppermost distribution range of A. adenophora (Datta et al 2017).

Our ENMs were built based on the occurrence data of invaded regions only, which indicates that the potential niche that we estimated is part of the fundamental niche of the species (Elith and Leathwick 2009; Soberon and Nakamura 2009). Moreover, factors such as biotic interactions and dispersal limitations also restrict the species from occupying its full potential niche (Soberon and Nakamura 2009). Thus, use of occurrence data from both native and invaded ranges would provide a more accurate estimate of the potential niche of the species (Jimenez-Valverde et al 2011).

In current climatic conditions, as well as in future climate scenarios, Middle Mountain is found to have more climatically suitable areas for A. adenophora compared to other physiographic regions. The elevation of the Middle Mountain region ranges from 1000 to 2500 masl (DHM 2017), which lies within the suitable range for the distribution of this weed (Wang and Wang 2006; Zhu et al 2007). A recent study that modeled the distribution of 24 invasive alien plants in Nepal also identified the Middle Mountain region as having a particularly large area of invasion hotspots, with suitable areas for the greatest number of species studied (Shrestha and Shrestha 2019). Displacement of native species, such as Artemisia indica and Urtica dioca, and reduction of the ground vegetation layer of Digitaria sp., Eragostris sp., and Imperata cylindrica by A. adenophora have already been observed in the Middle Mountain region (Tiwari et al 2005; Baral et al 2017). Therefore, this region requires the urgent attention of policymakers and land resource managers to implement effective management plans to prevent further spread of this weed. In contrast, the High Himalaya region is unsuitable for the weed because this region has a low minimum winter temperature (–10°C to 5°C) and low annual precipitation (400–1000 mm) (DHM 2017), which limit the distribution of the weed (Datta et al 2017). Most of the current areas predicted to be suitable for A. adenophora were found along road verges and river networks, which is in accordance with a study conducted in China (Wang and Wang 2006). These river and road networks enhance the rapid spread of A. adenophora (Dong et al 2008; Sang et al 2010). Therefore, control and monitoring efforts for the management of this weed should be focused on areas near rivers and roads.

There is growing evidence that climate change is likely to increase the risk of plant invasions, creating more suitable areas in the future (Bradley, Wilcove, et al 2010). Our projections also showed an increase in climatically suitable areas for A. adenophora in the future. The weed will gain suitable areas in 3 future climate scenarios: RCP 2.6 (2050) and RCP 4.5 (2050 and 2070). A recent study also reported that climatically suitable niches for A. adenophora would expand by 5.3% under RCP 6.0 in Nepal (Shrestha and Shrestha 2019). Studies in other parts of the world have also predicted the expansion of climatically suitable areas for A. adenophora in future climate scenarios (Wang and Wang 2006; Zhu et al 2007; Wang et al 2017; Lamsal et al 2018; Thapa et al 2018). However, a study conducted in Bhutan predicted a contraction of suitable areas for A. adenophora by 0.22% in 2050 (Thiney et al 2019). An experimental warming study conducted in China revealed that warming (2°C rise) increased biomass allocation and canopy cover of the weed, making it more stress tolerant (He et al 2012). Similarly, CO2 enrichment was also reported to increase the relative growth rate and biomass allocation of the weed (Lei et al 2012). Thus, an increase in suitable areas in future climate scenarios with increased temperature and CO2 concentration can be attributed to the innate and evolutionarily increased ecophysiological tolerances of A. adenophora favoring its growth (Blossey and Notzold 1995; He et al 2012; Lei et al 2012). Although, under extreme climate scenarios (RCP 8.5), a loss in total climatically suitable areas is predicted, the weed will still thrive and gain suitable areas in protected habitat and physiographic regions. Districts like Lamjung, Gorkha, Dhading, Makwanpur, Chitwan, and Tanahun are predicted to be vulnerable to further invasion by the weed due to climate change, so it is crucial to implement scientifically informed site-specific management policies, with the participation of local communities.

Despite its apomictic nature, with the associated evolutionary constraints, niche expansion has been observed in A. adenophora (Datta et al 2019). In addition, the weed is found to exhibit phenotypic plasticity that helps it to occupy a broader climatic niche (Zhao et al 2012). Our future climate models also demonstrate an expansion of the upper elevational distribution limit of the weed. Our results confirm findings in the Western Himalaya, where the weed was predicted to expand its upper elevational limit by 981 m compared to current climatic conditions (Thapa et al 2018). This indicates that the weed will spread toward cooler and drier regions in future. A similar trend was observed in a study of spatiotemporal patterns in China (Zhu et al 2007). Increased cold tolerance due to epigenetic modifications might help the weed to gain more suitable habitats in cooler and drier places at high elevations (Xie et al 2015). With climate change, all 3 protected areas will gain areas of potential suitable habitat. Though the weed has already been identified as the most problematic weed in Annapurna Conservation Area, impacting native diversity and livelihoods (Thapa and Maharjan 2014), no such studies have been carried out in Langtang National Park and Manaslu Conservation Area. Upward movement and colonization of A. adenophora due to recent climate changes have already been observed in Langtang National Park (Lamsal et al 2017). Our model also predicted that Langtang National Park has the highest proportion of climatically suitable areas in current climatic conditions and will continue to do so in future climate scenarios. A range shift of A. adenophora might threaten the habitat of 2 endangered animals, Ailurus fulgens (red panda) and Moschus chrysogaster (Himalayan musk deer), in Langtang National Park (Lamsal et al 2017). Thus, this information should act as a prompt for land managers, the scientific community, conservationists, and policymakers to adopt precautionary measures and formulate effective policies to prevent the further spread of this weed into new regions.

Taking climatic and topographic factors into consideration, our model predicted elevational range expansion, as well as an increase in suitable areas in future climate scenarios. However, other factors, such as biotic interactions, dispersal ability, demography, evolution, adaptation, and land-use change, also play key roles in determining the species range shift with climate change (Sinclair et al 2010; Urban et al 2016). Furthermore, long-term temporal predictions in climate change scenarios are associated with 2 other main errors. The first is the extrapolation of data beyond the training range to nonanalogue environmental conditions (climate change scenarios), which might make predictions unreliable (Fitzpatrick and Hargrove 2009). Another risk is that, as the relationship between climatic variables could change with time, the correlation structure of current and future climate variables might also change, thus reducing the certainty of models (Dormann et al 2013). To overcome these uncertainties and achieve more realistic predictions of species distribution, all factors, abiotic and biotic, that shape the distribution of invasive species should be incorporated in the species distribution models (Gonzalez-Salazar et al 2013; Leach et al 2016). Nevertheless, ENMs provide predictive information on species distribution required by vegetation managers and conservation practitioners for developing effective strategies to prevent further invasion (Peterson 2003). Thus, our study provides useful information about the current distribution of A. adenophora and identifies areas that may be at risk in the future on a local scale, demonstrating the urgent need for formulating effective management strategies to mitigate the impact of the weed.

Conclusions and management recommendations

Our study suggests that A. adenophora could spread further under future climate scenarios while retaining most of the currently suitable areas. Among the 4 physiographic regions in CHAL, the Middle Mountain region currently has the highest proportion and the High Himalaya region currently has the lowest proportion of climatically suitable areas for A. adenophora, and this is projected to continue under future projected climate scenarios. Similarly, Langtang National Park contains a higher percentage of areas suitable for the weed compared to the other 2 protected areas under current and future climate scenarios. All protected areas will gain additional suitable areas in future climate scenarios. Furthermore, it is predicted that the weed will expand its distribution range to higher elevations in future climate scenarios; this will amplify the consequences of climate change, which is already impacting these areas. For physiographic regions (Siwalik and High Himalaya) and protected areas (Manaslu Conservation Area) that have few suitable areas for invasion, regular inspection of habitats is needed to allow effective action to be taken in time to prevent further expansion of the weed. Management strategies for smaller and accessible invaded areas could either be mechanical control, for example, by hand pulling, or chemical control, by using herbicides (Parsons and Cuthbertson 2001; Di Tomaso et al 2013). Local communities in Nepal use A. adenophora for cattle bedding, composting, and also for making bio-briquettes (Tiwari et al 2005; Baral et al 2017; Shrestha, Shrestha, et al 2019); these cultural control methods can also be employed to manage the weeds, as long as precautions are taken to prevent seed dispersal. However, for regions like Middle Mountain and protected areas like Langtang National Park that have large areas at risk of invasion, effective management options, such as an integrated weed management approach, must be adopted and implemented along with regular monitoring of suitable habitats. By identifying areas that are potentially at risk in the future, our study constitutes a helpful resource for managers and policymakers to take appropriate and timely action to minimize the risk of invasion by A. adenophora associated with climate change.


This work was funded by the U.S. Agency for International Development (USAID) Bureau of Food Security under the Cooperative Agreement No. AID-OAA-L-15-00001 as part of the Feed the Future Innovation Lab for Integrated Pest Management. Any opinions, findings, conclusions, or recommendations expressed here are those of the authors alone. Some of the species location data used in this study were gathered by BBS during fieldwork supported by the International Foundation for Science (IFS, Sweden) and National Trust for Nature Conservation (Nepal). We express our gratitude to Dr DN Shah for his support during climatic data preparation. Thanks go to Ms Sara Hendery for English language editing in the manuscript. Comments and suggestions by anonymous reviewers greatly helped to improve the manuscript.



Aiello-Lammens ME, Boria RA, Radosavljevic A, Vilela B, Anderson RP . 2015. spThin: An R package for spatial thinning of species occurrence records for use in ecological niche models. Ecogeography 38:541–545. Google Scholar


Allouche O, Tsoar A, Kadmon R . 2006. Assessing the accuracy of species distribution models: Prevalence, kappa and true skill statistics (TSS). Journal of Applied Ecology 43:1223–1232. Google Scholar


Araujo MB, Guisan A . 2006. Five (or so) challenges for species distribution modeling. Journal of Biogeography 33:1677–1688. Google Scholar


Baral H, Inskipp C . 2005. Important Bird Areas in Nepal. Kathmandu, Nepal: Bird Conservation, Nepal. Google Scholar


Baral S, Adhikari A, Khanal R, Malla Y, Kunwar R, Basnyat B, Gauli K, Acharya RP . 2017. Invasion of alien plant species and their impact on different ecosystems of Panchase area, Nepal. Banko Janakari 27:31–42. Google Scholar


Bellard C, Thuiller W, Leroy B, Genovesi P, Bakkenes M, Courchamp F . 2013. Will climate change promote future invasions? Global Change Biology 19:3740–3748. Google Scholar


Bhattarai KR, Maren IE, Subedi SC . 2014. Biodiversity and invasibility: Distribution patterns of invasive plant species in the Himalayas, Nepal. Journal of Mountain Science 11:688–696. Google Scholar


Bhuju UR, Shakya PR, Basnet TB, Shrestha S . 2007. Nepal Biodiversity Resource Book: Protected Areas, Ramsar Sites, and World Heritage Sites. Kathmandu, Nepal: Ministry of Environment, Science and Technology and ICIMOD [International Centre for Integrated Mountain Development]. Google Scholar


Blossey B, Notzold R . 1995. Evolution of increased competitive ability in invasive non-indigenous plants: A hypothesis. Journal of Ecology 83:887–889. Google Scholar


Boria RA, Olson LE, Goodman SM, Anderson RP . 2014. Spatial filtering to reduce sampling bias can improve the performance of ecological niche models. Ecological Modelling 275:73–77. Google Scholar


BPP [Biodiversity Profile Project]. 1995. Biodiversity Profile of the High Himal and High Mountains Physiographic Zones. Publication No. 14. Kathmandu, Nepal: Department of National Parks and Wildlife Conservation, Ministry of Forests and Soil Conservation. Google Scholar


Bradley BA, Blumenthal DM, Wilcove DS, Ziska LH . 2010. Predicting plant invasions in an era of climate change. Trends in Ecology and Evolution 25:310–318. Google Scholar


Bradley BA, Wilcove DS, Oppenheimer M . 2010. Climate change increases risk of plant invasion in the Eastern United States. Biological Invasions 12:1855–1872. Google Scholar


Castro-Diez P, Vaz AS, Silva JS, van Loo M, Alonso A, Aponte C, Bayon A, Bellingham PJ, Chiuffo MC, DiManno N , et al. 2019. Global effects of non-native tree species on multiple ecosystem services. Biological Reviews 94:1477–1501. Google Scholar


CBS [Central Bureau of Statistics]. 2013. Preliminary Result of Census 2011. Kathmandu, Nepal: CBS. Google Scholar


Datta A, Kuhn I, Ahmad M, Michalski S, Auge H . 2017. Processes affecting altitudinal distribution of invasive Ageratina adenophora in western Himalaya: The role of local adaptation and the importance of different life-cycle stages. PLoS One 12:e0187708. Google Scholar


Datta A, Schweiger O, Kuhn I . 2019. Niche expansion of the invasive plant species Ageratina adenophora despite evolutionary constraints. Journal of Biogeography 46:1306–1315. Google Scholar


DHM [Department of Hydrology and Meteorology]. 2017. Observed Climate Trend Analysis in the Districts and Physiographic Regions of Nepal (1971–2014). Kathmandu, Nepal: DHM. Google Scholar


Di Tomaso JM, Kyser GB, Oneto SR, Wilson RG, Orloff SB, Anderson LW, Wright SD, Roncoroni JA, Miller TL, Prather TS , et al. 2013. Weed Control in Natural Areas in the Western United States. Berkeley, CA: Weed Research and Information Center, University of California. Google Scholar


DNPWC [Department of National Parks and Wildlife Conservation]. 2018. Annual Report. Babarmahal, Kathmandu: DNPWC. Google Scholar


Dong SK, Cui BS, Yang ZF, Liu SL, Liu J, Ding ZK, Zhu JJ, Yao WK, Wei GL . 2008. The role of road disturbance in the dispersal and spread of Ageratina adenophora along the Dian–Myanmar International Road. Weed Research 48:282–288. Google Scholar


Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carre G, Garcia Marquez JR, Gruber B, Lafourcade B, Leitao PJ , et al. 2013. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 36:27–46. Google Scholar


Dukes JS, Mooney HA . 1999. Does global change increase the success of biological invaders? Trends in Ecology and Evolution 14:135–139. Google Scholar


Elith J, Leathwick JR . 2009. Species distribution models: Ecological explanation and prediction across space and time. Annual Review of Ecology, Evolution and Systematics 40:677–697. Google Scholar


Elith J, Phillips SJ, Hastie T, Dudik M, Chee YE, Yates CJ . 2011. A statistical explanation of Maxent for ecologists. Diversity and Distributions 17:43–57. Google Scholar


Fielding AH, Bell JF . 1997. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environmental Conservation 24:38–49. Google Scholar


Fitzpatrick MC, Hargrove WW . 2009. The projection of species distribution models and the problem of non-analog climate. Biodiversity and Conservation 18:2255–2261. Google Scholar


Franklin J. 2009. Mapping Species DistributionsSpatial Inference and Prediction. Cambridge, United Kingdom: Cambridge University Press. Google Scholar


Gent PR, Danabasoglu G, Donner LJ, Holland MM, Hunke EC, Jayne SR, Lawrence DM, Neale RB, Rasch PJ, Vertenstein M , et al. 2011. The Community Climate System Model version 4. Journal of Climate 24:4973–4991. Google Scholar


Gonzalez-Salazar C, Stephens CR, Marquet PA . 2013. Comparing the relative contributions of biotic and abiotic factors as mediators of species distributions. Ecological Modelling 248:57–70. Google Scholar


He WM, Li JJ, Peng PH . 2012. A congeneric comparison shows that experimental warming enhances the growth of invasive Eupatorium adenophorum. PLoS One 7:e35681. Google Scholar


Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A . 2005. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25:1965–1978. Google Scholar


Hutchinson GE. 1957. Concluding remarks. Cold Spring Harbor Symposia on Quantitative Biology 22:145–199. Google Scholar


IPBES [Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services]. 2019. Summary for policymakers of the Global Assessment Report on Biodiversity and Ecosystem Services . Díaz S, Settele J, ES Brondízio ES, Ngo HT, Guèze M, Agard J, Arneth A, Balvanera P, Brauman KA, Butchart SHM, et al, editors. Bonn, Germany: IPBES Secretariat. Google Scholar


IPCC [Intergovernmental Panel on Climate Change]. 2013. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change . Stocker TF, Qin D, Plattner GK, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM, editors. Cambridge, United Kingdom, and New York, NY: Cambridge University Press. Google Scholar


Jarnevich CS, Reynolds LV . 2011. Challenges of predicting the potential distribution of a slow invading invader: A habitat suitability map for an invasive riparian tree. Biological Invasions 13:153–163. Google Scholar


Jarnevich CS, Stohlgren TJ, Barnett D, Kartesz J . 2006. Filling in the gaps: Modeling native species richness and invasions using spatially incomplete data. Diversity and Distributions 12:511–520. Google Scholar


Jimenez-Valverde A, Peterson AT, Soberon J, Overton JM, Aragon P, Lobo JM . 2011. Use of niche models in invasive species risk assessments. Biological Invasions 13:2785–2797. Google Scholar


Kathiresan R, Gualbert G . 2016. Impact of climate change on the invasive traits of weeds. Weed Biology and Management 16:59–66. Google Scholar


Lamsal P, Kumar L, Atreya K . 2017. Historical evidence of climatic variability and changes, and its effect on high-altitude regions: Insights from Rara and Langtang, Nepal. International Journal of Sustainable Development and World Ecology 24:471–484. Google Scholar


Lamsal P, Kumar S, Aryal A, Atreya K . 2018. Invasive alien plant species dynamics in the Himalayan region under climate change. Ambio 47:697–710. Google Scholar


Leach K, Montgomery WI, Reid N . 2016. Modelling the influence of biotic factors on species distribution patterns. Ecological Modelling 337:96–106. Google Scholar


Lei YB, Wang WB, Feng YL, Zheng YL, Gong HD . 2012. Synergistic interactions of CO2 enrichment and nitrogen deposition promote growth and ecophysiological advantages of invading Eupatorium adenophorum in Southwest China. Planta 236:1205–1213. Google Scholar


Liu Y, Oduor AMO, Zhang Z, Manea A, Tooth IM, Leishman MR, Xu X, van Kleunen M . 2017. Do invasive alien plants benefit more from global environmental change than native plants? Global Change Biology 23:3363–3370. Google Scholar


Maharjan S, Shrestha BB, Joshi MD, Devkota A, Muniappan R, Adiga A, Jha PK . 2019. Predicting suitable habitat of an invasive weed Parthenium hysterophorus under future climate scenarios in Chitwan Annapurna Landscape, Nepal. Journal of Mountain Science 16:2243–2256. Google Scholar


Mainka SA, Howard GW . 2010. Climate change and invasive species: Double jeopardy. Integrative Zoology 5:102–111. Google Scholar


Manel S, Ceri Williams H, Ormerod SJ . 2001. Evaluating presence–absence models in ecology: The need to account for prevalence. Journal of Applied Ecology 38:921–931. Google Scholar


Marbuah G, Gren IM, McKie B . 2014. Economics of harmful invasive species: A review. Diversity 6:500–523. Google Scholar


McInerny GJ, Etienne RS . 2012. Ditch the niche—Is the niche a useful concept in ecology or species distribution modelling? Journal of Biogeography 39:2096–2102. Google Scholar


Merow C, Smith MJ, Silander JA . 2013. A practical guide to MaxEnt for modeling species' distributions: What it does, and why inputs and settings matter. Ecography 36:1058–1069. Google Scholar


MFSC [Ministry of Forest and Soil Conservation]. 2016. Conservation Landscapes of Nepal. Singha Durbar, Kathmandu, Nepal: MFSC. Google Scholar


Paini DR, Sheppard AW, Cook DC, De Barro PJ, Worner SP, Thomas MB . 2016. Global threat to agriculture from invasive species. Proceedings of the National Academy of Sciences of the United States of America 113:7575–7579. Google Scholar


Parsons WT, Cuthbertson EG . 2001. Noxious Weeds of Australia. 2nd edition (1st edition 1992). Collingwood, Australia: CSIRO Publishing. Google Scholar


Peterson AT. 2003. Predicting the geography of species invasion via ecological niche modeling. Quarterly Review of Biology 78:419–433. Google Scholar


Phillips SJ. 2017. A Brief Tutorial on Maxent. New York, NY: American Museum of Natural History.; accessed on 16 September 2017. Google Scholar


Phillips SJ, Anderson RP, Schapire RE . 2006. Maximum entropy modeling of species geographic distributions. Ecological Modelling 190:231–259. Google Scholar


Phillips SJ, Dudik M . 2008. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecogeography 31:161–175. Google Scholar


Poudel AS, Jha PK, Shrestha BB, Muniappan R . 2019. Biology and management of the invasive weed Ageratina adenophora (Asteraceae): Current state of knowledge and future research needs. Weed Research 59:79–92. Google Scholar


Qin Z, Zhang JE, DiTommaso A, Wang RL, Liang KM . 2016. Predicting the potential distribution of Lantana camara L under RCP scenarios using ISI-MIP models. Climatic Change 134:193–208. Google Scholar


R Core Team. 2017. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.; accessed on 16 September 2017. Google Scholar


Radosavljevic A, Anderson RP . 2014. Making better Maxent models of species distributions: Complexity, overfitting and evaluation. Journal of Biogeography 41:629–643. Google Scholar


Sang W, Zhu L, Axmacher JC . 2010. Invasion pattern of Eupatorium adenophorum Spreng in southern China. Biological Invasions 12:1721–1730. Google Scholar


Shrestha BB. 2019. Management of invasive alien plants in Nepal: Current practices and future prospects. In : Garkoti S, Van Bloem S, Fule P, Semwal R, editors. Tropical Ecosystems: Structure, Functions and Challenges in the Face of Global Change . Singapore: Springer, pp 45–68. Google Scholar


Shrestha BB, Budha PB, Pagad S, Wong LJ . 2019. Global Register of Introduced and Invasive Species—Nepal. ISSG [Invasive Species Specialist Group] Checklist Dataset. Copenhagen, Denmark: Global Biodiversity Information Facility.; accessed on 8 November 2019. Google Scholar


Shrestha BB, Shrestha UB, Sharma KP, Thapa-Parajuli RB, Devkota A, Siwakoti M . 2019. Community perception and prioritization of invasive alien plants in Chitwan–Annapurna Landscape, Nepal. Journal of Environmental Management 229:38–47. Google Scholar


Shrestha UB, Sharma KP, Devkota A, Siwakoti M, Shrestha BB . 2018. Potential impact of climate change on the distribution of six invasive alien plants in Nepal. Ecological Indicators 95:99–107. Google Scholar


Shrestha UB, Shrestha BB . 2019. Climate change amplifies plant invasion hotspots in Nepal. Diversity and Distributions 25:1599–1612. Google Scholar


Sinclair SJ, White MD, Newell GR . 2010. How useful are species distribution models for managing biodiversity under future climates? Ecology and Society 15:art8. Google Scholar


Siwakoti M, Shrestha BB, Devkota A. 2016. Assessment of the effects of climate change on distribution of invasive alien plant species in Nepal. In : Bhuju DR, McClaughlin K, Sijapati J, Devkota BD, Shrestha N, Ghimire GP, Neupane PK, editors. Building Knowledge for Climate Resilience in Nepal: Research Brief. Kathmandu, Nepal: Nepal Academy of Science and Technology, pp 66–72. Google Scholar


Soberon J, Nakamura M. 2009. Niches and distributional areas: Concepts, methods, and assumptions. Proceedings of the National Academy of Sciences of the United States of America 106:19644–19650. Google Scholar


Stacowicz JJ, Terwin JR, Whitlatch RB, Osman W . 2002. Linking climate change and biological invasions: Ocean warming facilitates nonindigenous species invasion. Proceedings of the National Academy of Sciences of the United States of America 99:15497–15500. Google Scholar


Suarez-Mota ME, Oritz E, Villasenor JL, Espinosa-Garcia FJ . 2016. Ecological niche modeling of invasive plant species according to invasion status and management needs: The case of Chromolaena odorata (Asteraceae) in South Africa. Polish Journal of Ecology 64:369–383. Google Scholar


Svenning JC, Normand S, Kageyama M . 2008. Glacial refugia of temperate trees in Europe: Insights from species distribution modeling. Journal of Ecology 96:1117–1127. Google Scholar


Swets JA. 1988. Measuring the accuracy of diagnostic systems. Science 240:1285–1293. Google Scholar


Tererai F, Wood AR . 2014. On the present and potential distribution of Ageratina adenophora (Asteraceae) in South Africa. South African Journal of Botany 95:152–158. Google Scholar


Thapa GJ, Wikramanayake E, Forrest J . 2015. Climate-Change Impacts on the Biodiversity of the Terai Arc Landscape and the Chitwan–Annapurna Landscape. Kathmandu, Nepal: Hariyo Ban, WWF [World Wide Fund for Nature]. Google Scholar


Thapa N, Maharjan M . 2014. Invasive alien species: Threats and challenges for biodiversity conservation—A case study of Annapurna Conservation Area, Nepal. In : Thapa GJ, Subedi N, Pandey MR, Thapa SK, Chapagain NR, Rana A, editors. Proceedings of the International Conference on Invasive Alien Species Management. Kathmandu, Nepal: National Trust for Nature Conservation, pp 18–22. Google Scholar


Thapa S, Chitale V, Rijal SJ, Bisht N, Shrestha BB . 2018. Understanding the dynamics in distribution of invasive alien plant species under predicted climate change in Western Himalaya. PLoS One 13:e0195752. Google Scholar


Thiney U, Banterng P, Gonkhamdee S, Katawatin R . 2019. Distributions of alien invasive weeds under climate change scenarios in mountainous Bhutan. Agronomy 9:442. Google Scholar


Tiwari S, Adhikari B, Siwakoti M, Adhikari B, Subedi K . 2005. An Inventory and Assessment of Invasive Alien Species of Nepal. Kathmandu, Nepal: IUCN [International Union for Conservation of Nature]. Google Scholar


Urban MC, Bocedi G, Hendry AP, Mihoub JB, Peer G, Singer A, Bridle JR, Crozier LG, De Meester L, Godsoe W , et al. 2016. Improving the forecast for biodiversity under climate change. Science 353:aad8466. Google Scholar


van Vuuren DP, Edmonds J, Kainuma M, Riahi K, Thomson A, Hibbard K, Hurtt GC, Kram T, Krey V, Lamarque JF , et al. 2011. The representative concentration pathways: An overview. Climatic Change 109:5–31. Google Scholar


Villa M, Hulme PE . 2017. Non-native species, ecosystem services, and human well-being. In : Vila M, Hulme PE, editors. Impact of Biological Invasions on Ecosystem Services. Bern, Switzerland: Springer, pp 1–14. Google Scholar


Walther GR, Roques A, Hulme PE, Sykes MT, Pysek P, Kuhn I, Zobel M . 2009. Alien species in a warmer world: Risks and opportunities. Trends in Ecology and Evolution 24:686–693. Google Scholar


Wan JZ, Wang CJ, Tan JF, Yu FH . 2017. Climatic niche divergence and habitat suitability of eight alien invasive weeds in China under climate change. Ecology and Evolution 7:1541–1552. Google Scholar


Wang C, Lin H, Feng Q, Jin C, Cao A, He L . 2017. A new strategy for the prevention and control of Eupatorium adenophorum under climate change in China. Sustainability 9:2037. Google Scholar


Wang R, Wang YZ . 2006. Invasion dynamics and potential spread of the invasive alien plant species Ageratina adenophora (Asteraceae) in China. Diversity and Distribution 12:397–408. Google Scholar


Wenger SJ, Olden JD . 2012. Assessing transferability of ecological models: An underappreciated aspect of statistical validation. Methods in Ecology and Evolution 3:260–267. Google Scholar


West AM, Kumar S, Brown CS, Stohlgren TJ, Bromberg J . 2016. Field validation of an invasive species Maxent model. Ecological Informatics 36:126–134. Google Scholar


WWF [World Wide Fund for Nature]. 2013. Chitwan Annapurna Landscape: A Rapid Assessment Hariyo Ban Program. Kathmandu, Nepal: WWF. Google Scholar


Xie HJ, Li H, Liu D, Dai WM, He JY, Lin S, Duan H, Liu LL, Chen SG, Song XL , et al. 2015. ICE1 demethylation drives the range expansion of a plant invader through cold tolerance divergence. Molecular Ecology 24:835–850. Google Scholar


Zhao X, Liu W, Zhou M . 2012. Lack of local adaptation of invasive crofton weed (Ageratina adenophora) in different climatic areas of Yunnan Province, China. Journal of Plant Ecology 6:316–322. Google Scholar


Zhu L, Sun OJ, Sang W, Li Z, Ma K . 2007. Predicting the spatial distribution of invasive plant species (Eupatorium adenophorum) in China. Landscape Ecology 22:1143–1154.; assessed on 16 December 2018. Google Scholar


Supplemental material

 APPENDIX S1 (mred-40-02-15_s01.pdf) Bioclimatic variables used for modeling suitable habitat of Ageratina adenophora.

 APPENDIX S2 (mred-40-02-15_s01.pdf) Correlation matrix of 19 bioclimatic and 3 topographic variables.

 APPENDIX S3 (mred-40-02-15_s01.pdf) Marginal response curves of 2 predictor variables.

 APPENDIX S4 (mred-40-02-15_s01.pdf) Scatter plot of occurrence points and minimum temperature of the coldest month for current climate.

Found at:

© 2020 Poudel et al. This open access article is licensed under a Creative Commons Attribution 4.0 International License ( Please credit the authors and the full source.
Anju Sharma Poudel, Bharat Babu Shrestha, Mohan Dev Joshi, Rangaswamy Muniappan, Abhijin Adiga, Srinivasan Venkatramanan, and Pramod Kumar Jha "Predicting the Current and Future Distribution of the Invasive Weed Ageratina adenophora in the Chitwan–Annapurna Landscape, Nepal," Mountain Research and Development 40(2), R61-R71, (30 November 2020).
Received: 1 March 2020; Accepted: 1 June 2020; Published: 30 November 2020

Back to Top