Anthropogenic activities and global climate change increase the risk of Chromolaena odorata invasion and habitat expansion. To predict its global distribution and habitat suitability under climate change, a random forest (RF) model was employed. The RF model, utilizing default parameters, analyzed species presence data and background information. The model revealed that the current spatial distribution of C. odorata covers 7,892,447 km². Predictions for 2061– 2080 indicate expansion of suitable habitat (42.59 and 46.30%), reduction of suitable habit (12.92 and 12.20%), and preservation of suitable habitat (87.08 and 87.80%) under the SSP (Shared Socio-economic Pathway) 2-4.5 and SSP5-8.5 scenarios, respectively, in comparison to the present distribution. Currently, C. odorata is predominantly found in South America, with limited presence in other continents. However, the data suggest that climate change will elevate the global invasion risk of C. odorata worldwide, particularly in Oceania, Africa, and Australia. Countries such as Gambia, Guinea-Bissau, and Lesotho, which currently have unsuitable habitats, are predicted to have highly suitable habitats with climate change, supporting the idea that global habitat expansion for C. odorata will occur due to climate change. This study indicates that proper management of C. odorata is crucial during the early invasion phase.

Species distribution modeling is widely used for evaluating invasion risk, and for prioritizing areas for the control and management of invasive species. However, selecting a modeling tool that accurately predicts species invasion risk requires a systematic approach. In this study, five species distribution models (SDMs), namely, artificial neural network (ANN), generalized linear model (GLM), multivariate adaptive regression splines (MARS), maximum entropy (MaxEnt), and random forest (RF), were performed and evaluated their model performance using the mean value of area under the curve (AUC), true skill statistics (TSS), and Kappa scores of 12 ecosystem disturbing alien plant species (EDAPS). The mean evaluation metric scores were highest in RF (AUC = 0.924 ± 0.058, TSS = 0.789 ± 0.109, Kappa = 0.671 ± 0.096, n = 12) and lowest in ANN. The ANOVA of AUC, TSS, and Kappa metrics revealed the RF model was significantly different from other SDMs and was therefore selected as the relatively best model. The potential distribution area and invasion risk for each EDAPS were quantified. Under the current climate conditions of South Korea, the average potential distribution area of EDAPS was estimated to be 13,062 km². However, in future climate change scenarios, the average percentage change of EDAPS distribution relative to the current climate was predicted to be increased over 219.93%. Furthermore, under the current climate, 0.16% of the area of the country was estimated to be under a very high risk of invasion, but this would increase to 60.43% by 2070. Invasion risk under the current climate conditions was highest in the northwestern, southern, and southeastern regions, and in densely populated cities, such as Seoul, Busan, and Daegu. By 2070, invasion risk was predicted to expand across the whole country except in the northeastern region. These results suggested that climate change induced the risk of EDAPS invasiveness, and SDMs could be valuable tools for alien and invasive plant species risk assessment.

The global climate change, including increases in temperature and precipitation, may exacerbate the invasion by P. hysterophorus. Here, MaxEnt modeling was performed to predict P. hysterophorus distribution worldwide and in South Korea under the current and future climate global climate changes, including increases in temperature and precipitation. Under the current climate, P. hysterophorus was estimated to occupy 91.26%, 83.26%, and 62.75% of the total land area of Australia, South America, and Oceania, respectively. However, under future climate scenarios, the habitat distribution of P. hysterophorus would show the greatest change in Europe (56.65%) and would extend up to 65°N by 2081–2100 in South Korea, P. hysterophorus currently potentially colonizing 2.24% of the land area, particularly in six administrative divisions. In the future, P. hysterophorus would spread rapidly, colonizing all administrative divisions, except Incheon, by 2081–2100. Additionally, the southern and central regions of South Korea showed greater habitat suitability than the northern region. These findings suggest that future climate change will increase P. hysterophorus distribution both globally and locally. Therefore, effective control and management strategies should be employed around the world and in South Korea to restrict the habitat expansion of P. hysterophorus.

Predicting the regions at risk of invasion from IIAPS is an integral horizon-scanning activity that plays a crucial role in preventing, controlling, and eradicating invasive species. Here, we quantify the spatial distribution area and invasion risk of IIAPS using a species distribution model under different levels of environmental change in South Korea. From the model predictions, the current average spatial extent of the 10 IIAPS is 33,948 km², and the individual spatial extents are estimated to change by −7% to 150% by 2050 and by −9% to 156% by 2070. The spatial invasion risk assessment shows that, currently, moderate-to-high invasion risk is limited to coastal areas and densely populated metropolitan cities (e.g., Seoul, Busan, and Gwangju), but that the area with this level of risk is expected to spread toward the central and northern regions of the country in the future, covering 86.21% of the total area of the country by 2070. These results demonstrate that the risk of invasion by IIAPS is estimated to enlarge across the whole country under future environmental changes. The modeling system provided in this study may contribute to the initial control and strategic management of IIAPS to maintain the dynamic ecosystems of South Korea.