Mapping impacts of alien species on biodiversity in the face of climate change

Mapping impacts of alien species on biodiversity in the face of climate change | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Mapping impacts of alien species on biodiversity in the face of climate change Sabrina Kumschick, Lysandre Journiac, Océane Boulesnane-Genguant, and 3 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4437291/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 18 Dec, 2024 Read the published version in Biological Invasions → Version 1 posted 4 You are reading this latest preprint version Abstract Risk maps are a useful tool to prioritise sites for management and allocate resources where they are most needed as they can show us where impacts of biological invasions are most likely to happen or expected to be largest. Given the pace of global changes, we need to understand not only the risks under current conditions, but future risks taking into account aspects such as climate change. In this study, we use Australian acacias alien to South Africa as a case study to model their potential distribution under climate change scenarios, and map their potential impacts and the uncertainty related to variability across socio-economic pathways and climatic models, to help us set better priorities for the future. We found that although climatic suitability and therefore the potential high risk area is predicted to decrease under climate change, the highly vulnerable fynbos biome remains an area with high projected impacts. We found that potential impacts would decrease under climate change within protected areas of the Western Cape. However, in this province, national parks and other areas with high levels of protection would experience greater impacts than private protected areas of lower protection status. This can help us prioritise management actions and aid the development of suitable plans to protect biodiversity not only under current, but also future climate conditions. EICAT biological invasions Australian Acacia SDM protected area risk map Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Introduction Biological invasions and climate change are two of the main drivers of global biodiversity loss (Master and Norgrove 2010; Johnson et al. 2017; IPBES 2023 ). Impacts of biological invasions are manifold, including changes to ecosystem functions such as fire regimes and hydrology and changes to biodiversity including extirpations and extinctions of native species (Vilà et al. 2011 ; Smith 2020 ; Vilà et al. 2024 ; IPBES 2023 ). Similarly, climate change can have severe impacts on native species, for example increasing their risk of extinctions (e.g., Thomas et al. 2004 ; Manes et al. 2021 ) and causing range changes (e.g., Lenoir et al. 2008 ; Doak et al. 2010). While these drivers can act in isolation, climate change can also affect the invasion of alien species as climatic conditions play an important role in the growth, survival and distribution of species (Walther et al. 2009; Master and Norgrove 2010). Climate change can not only lead to changes in native species distributions, but also those of alien species (Bellard et al. 2013 ). Member states of the Convention on Biological Diversity (CBD) are mandated to improve the protection of natural resources under the Kunming-Montreal Global Biodiversity Framework (GBF). Target 1 asks to minimize the loss of areas with high biodiversity. To achieve this, we need to understand the pressures on these areas from drivers such as climate change and biological invasions, and map them to prioritise control. Risk maps have been used for various purposes in invasion science. For example, they can be useful to improve detection of new incursions or individuals of species with limited populations (Kaplan et al. 2014 ). They have also been used to map the potential impacts to assess sites experiencing the highest current or potential impact, translating to sites at risk of invasion impacts (Nentwig et al. 2010 ). However, as species’ distributions are expected to change with climate change, considering only current climatic conditions in such risk maps does not prepare us well for future risk scenarios. Therefore, to improve our ability to manage biological invasions, as stipulated under Target 6 of the GBF, we need to understand how biological invasions and their impacts could change under climate change. Simultaneously, this addresses Target 8, which aims to minimize the effect of climate change on biodiversity. Still, climate change is not routinely included in risk assessments for alien taxa (but see Marchioro and Krechemer 2021 for some insect pests). Correlative species distribution models (hereafter SDM) are popular in various fields of ecology and are increasingly used in invasion science (Guisan et al. 2014 ). They are statistical methods which correlate known occurrences of species with environmental variables and predict a species' potential distribution in other regions over space and time (Guisan and Zimmermann 2000 ). The advantages of SDMs are manifold as they are easy to implement and can cover large geographic areas. Therefore, they are widely used to produce maps of potential invasion by delineating areas based on climatic suitability for a species (e.g. Bradley et al. 2010 ; Jiménez-Valverde et al. 2011 ; Tingley et al. 2017 ). Despite known limitations (see for example Hui 2023 ), they have shown potential for accurate prediction of alien species spread (Barbet-Massin et al. 2018 ) and response of biological invasions to climate change (Sheppard et al. 2014 ) if properly calibrated. In this study, we use Australian Acacia species alien to South Africa as a case study to assess i) the current distribution of the species and their potential impacts, ii) the potential distribution under current climate, and iii) the projected future distribution under climate change scenarios. By overlaying the outputs of the models with the species potential impacts, we produce risk maps for current and projected future climatic conditions and assess potential impacts for protected areas of the Western Cape Province in South Africa. Methods We detail below how we selected Australian Acacia species and fitted SDM models per species that were then used to project the habitat suitability across South Africa per time horizon, accounting for various climatic scenarios. Then, we explain how we aggregated an impact score from the species-wise SDM projections and the documented impacts per species. Species selection We selected, as a case study, Australian Acacia species ( Acacia Mill. sensu stricto – synonym: Acacia subg. Phyllodineae (DC.) Seringe) which are alien in South Africa. This includes 33 species (Magona et al. 2018 ; Jansen and Kumschick 2022 ; Kumschick and Jansen 2023 ; Supplementary Material Appendix S4.1), of which three species ( Acacia crassiuscula , A. acuminata and A. koa ) had to be excluded from all analysis because they had less than 10 presence records after the filtering steps detailed below. Acacias are among the most highly impacting alien taxa in South Africa with a diverse range of impacts on native biodiversity (le Maître et al. 2011; Kumschick & Jansen 2023 ). Occurrence data We gathered distribution data from the native and introduced ranges of the 30 Acacia species, and used it along with bioclimatic variables to build individual ensemble models which were then projected based on current and future climate in South Africa. For each species, worldwide occurrence data were obtained from the Global Biodiversity Information Facility (GBIF, http://www.gbif.org/ ) using the rgbif package (Chamberlain et al. 2022 ). For this extraction, we kept all records which had complete coordinates, which were tagged as “presence”, and which were either human or machine observations or coming from the literature. No time limit was specified. Then we filtered occurrences to reduce the uncertainty arising from questionable records and only kept records which came from a list of 65 trusted GBIF datasets. This list was previously developed for the global distribution of Acacia species (Botella et al. 2023 ) and is composed of datasets which include at least one documented step of taxonomic verification in the recording of the occurrences. The documented native and introduced countries of each species were also obtained from that study. Besides, we extracted the global land cover raster (approx. 300 m resolution) from the European Space Agency GlobCover product. We excluded records located outside the native and introduced countries of the species (obtained from Botella et al. 2023 ) or whose 4 km x 4 km cell was located in urban areas based on the GlobCover 2009 global land cover raster (ESA and UCLouvain GlobCover 2009 Project, last accessed 29/01/2024). The latter step was implemented to exclude planted trees. More precisely, we excluded cells classified as “Artificial surfaces and associated areas (Urban areas > 50%)” by GlobCover. We then resampled our remaining records to keep only up to one record per 4 km × 4 km cell for each species to reduce the oversampling of some areas (Boria et al. 2014 ). We obtained a total of 37,940 filtered records (Figure S4 .1 & S4.2). The number of filtered records were imbalanced across the 30 considered species, with 27 for the least sampled ( A. adunca ) and 3,738 for the most sampled ( A. dealbata ). Climate variables To estimate the climatic suitability for each Acacia species from their occurrence data, we extracted four bioclimatic variables from the Climatologies at High Resolution for Earth’s Land Surface Areas database (CHELSA, Karger et al. 2017 ). We chose the mean annual air temperature (bio1), mean diurnal air temperature range (bio2), annual precipitation (bio12) and precipitation seasonality (bio15) as environmental variables as they are likely to affect plant productivity and survival across the landscape (Mod et al., 2016 ), and the correlation between each of these variables was moderate (Pearson’s | r | < 0.55). These data are available globally at a 4 km × 4 km resolution for the period 1981–2010. Then, we predicted the suitable climatic range of each species per time period in South Africa. The “current” South African climate was derived from the 1981–2010 period and the “future” climate was decomposed in two horizons: mid-century (2041–2070) and end of the century (2071–2100). For both time horizons, we extracted the four climatic variables based five Global Climate Models from CHELSA (“gfdl-esm4”, “ipsl-cm6a”, “mpi-esm1”, “mri-esm2”, “ukesm1”) for three shared socio-economic pathway scenarios (“ssp1.26”, “ssp3.70” and “ssp5.85”) for the South African region: ssp1.26: scenario with low greenhouse gas (GHG) emission; ssp3.70: scenario with high GHG emission; ssp5.85: scenario with very high GHG emission (Fig. 1 and Supplementary Material, Appendix S1). Species distribution models Pseudo-absences selection For each species we drew ten datasets of pseudo-absences composed of the same number of pseudo-absences as presences. For species with less than 100 presences, we increased the number of pseudo-absences to 100 in each dataset, to reduce the variance (over-fitting) of each model. These pseudo-absences were drawn in the countries being part of the native and introduced ranges of the species assuming that the species likely had an opportunity to establish at pseudo-absence locations. Following recommendations of Barbet-Massin et al. ( 2012 ), pseudo-absences were selected based on the surface range envelope method (Busby 1991 ), by randomly selecting pseudo-absences outside the usual climatic range of the species occurrences (Thuiller et al. 2009 ). Given that pseudo-absences were drawn in the native and introduced ranges but outside of the climatic envelope where the species occurrences were recorded (PA.sre.quant = 0), we can assume false absences to be rare. Therefore, the final model predictions will approximate the probability of presence conditional to the species presence in the area, i.e. an establishment probability. Model implementation For each species we built a species distribution model using functions and workflow of the biomod2 package (Thuiller et al. 2009 ). We applied a random forest classifier (Liaw and Wiener 2002 ) on the presence/pseudo-absence datasets and the four climatic variables to decipher which climatic conditions are suitable for each species. Each forest was composed of 500 trees, where each branch tests two randomly picked variables (mtry = 2), and we impose a minimum of five sampled presence/pseudo-absences per terminal node. We implemented a step of cross-validation for which the presence/pseudo-absence dataset was divided into four blocks along the longitude axis with presences equally balanced between each block (Wenger and Olden 2012 ). Then, three blocks were used for model training and one for testing. This approach allows for spatial and environmental independence between training and testing sets. Kappa and TSS metrics were computed each time on the testing dataset. As we had ten pseudo-absence datasets and four blocks, we computed 40 random forest models for each species. For each species we kept only the random forest models with a Kappa and a TSS over 0.7, indicating a robust predictive capacity, and which we refer to as sub-models. We then derived for each species an ensemble of its sub-models that was a “committee averaging” of the binarized responses given by the sub-models (Thuiller et al. 2009 ). We used a probability threshold of 0.5 to binarize each submodel response, corresponding to a majority vote among trees of the corresponding random forest. Committee averaged values are the proportion of sub-models predicting the species as present. Here, we considered that committee averaged values can be regarded as a proxy for the probability of a species to establish. Thus, we obtained one projection of establishment probability for South Africa under the current climate and 15 scenarios (five global climate models x three socio-economic pathways) for each time horizon (2041–2070 and 2071–2100) as represented in Fig. 1 . So, for each time horizon and each spatial cell, we averaged the 15 establishment probabilities of the 15 scenarios to obtain the final establishment probability of the species in that cell. This final establishment probability gives an equal weight to the 15 plausible scenarios. Finally, for each time horizon, we produced a potential Acacia species richness map summing the 30 species establishment probability maps, and an associated uncertainty map. The uncertainty map shows standard deviation of the potential richness estimate based on the variability standard deviation of the predicted establishment per species and its variability across the 15 scenarios. Risk maps As we aimed to map the accumulated risk of Acacia species potentially causing impacts in different regions in South Africa, we firstly derived a potential impact score per species based on its documented impacts. Data on the impacts of the Acacia species were taken from Jansen and Kumschick ( 2022 ). They assigned categories of Minimal Concern to Massive to each species, based on the impacts reported in the literature and following the International Union for the Conservation of Nature (IUCN) standard for impact classification, the Environmental Impact Classification for Alien Taxa, EICAT (Blackburn et al. 2014 ; Hawkins et al. 2015 ; IUCN 2020a , b ; Volery et al. 2020 ). We used the maximum impact category per species as a precautionary measure, which is the recommended approach used by the IUCN EICAT Standard (IUCN 2020a ; Kumschick et al. 2024 ), noting that other calculation methods are possible depending on the aim of the study (see also Kumschick et al. 2024 , Boulesnane-Genguant et al. in prep). We transformed the impact categories into numerical values as follows: Minimal Concern (MC) = 0, Minor (MN) = 1, Moderate (MO) = 2, Major (MR) = 3, Massive (MV) = 4. 19 species classified as Data Deficient (DD), i.e. no data on impact was found in a comprehensive literature search, and they were assigned a score of 0 (Supplementary Material Appendix S4.1). As a prediction of local impact risk per species, we multiplied the predicted establishment probability derived from the SDMs in each grid cell by the species impact score. Finally, we summed up the local impact risks across species in each grid cell to obtain the impact risk map (see also Nentwig et al. 2010 ). Hence, each species occurring adds to the risk of the potential impacts actually manifesting at that site, again using a precautionary approach. Our maps therefore do not represent sites where impacts are definitely going to occur, as the potential establishment of each species alone does not necessarily happen and, even if it does, it does not necessarily result in the worst impact documented for that species. Other factors such as the recipient community and the abundance of the alien species play important roles in the manifestation of impacts. Note also that even if data deficient species are present, they do not affect the impact risk value, as for species whose maximum impact is MC. This might underestimate the impact of some species which actually cause higher impacts but which have not been studied. Changes in impact risk in protected areas To assess the specific change in risk across time on protected areas, we assessed the potential changes of impact risk for protected areas in the Western Cape Province of South Africa. For each protected area, we calculated the mean impact risk score and compared it under current and future (2071–2100) projected climate. To see whether differences could be found based on the protected area status, we also calculated the mean impact risk score per protected area type. In the Western Cape, protected areas are classified as follows (with increasing level of protection): private nature reserves, provincial nature reserves and national parks. We obtained protected areas data from the South African National Biodiversity Institute ( www.bgis.sanbi.org ). Results Firstly, we present a method to map the potential cumulative impacts of alien species under climate change projections (Fig. 1 , Supplementary Material, Appendix S1). It combines the recorded impacts alien species have had with the projected species richness. Applying this method to Acacia species in South Africa, we find that across all species currently present in the country, the suitable area declines over time and only the south coast in the Eastern Cape province and parts of the Western Cape province still remain suitable for a considerable number of species (Fig. 2 ). Large parts of the north eastern (NE) part of the country are predicted to become unsuitable for most wattles in future. Indeed, while the current projected impact score lies between 15 and 25 in a relatively large part of the NE area, it is inferior to 5 in nearly all of the area for the end of the century (Fig. 2 ). Yet, this decrease should be interpreted with caution, as the standard deviation of the end-of-century impact score, due to the variability across socio-economic pathways and climate models, is around 2 in the NE area that currently shows a high risk (Figure S3 ). Comparing the current impact risk map to a risk map based on actual records for Acacia species (Fig. 3 ), the general pattern looks similar to potential impacts under current climate, with a less continuous distribution along the east coast, and a marked region of high risk in the Gauteng province. The latter is probably due to the fact that in the map showing impacts based on recorded occurrences, urban areas were not excluded, and much of the Gauteng province would fall under this category (Fig. 3 ). Furthermore, a positive sampling bias in this highly populated area might be responsible for higher detection rates of the species present and explain the relatively high impact risk when based on actual occurrences. [insert Fig. 2 here] [insert Fig. 3 here] Although species richness of Acacias is not projected to be as high in the southwestern part of the country as along the south coast, the projected risks are similar in both areas (Fig. 2 ). That shows that the areas around Cape Town and Hermanus, although not suitable for as many species, might remain suitable for the high impacting species. The impact score in this area remains most likely high compared to other areas even when considering its standard deviation (Figure S3 ). We found that potential impacts would decrease under climate change within protected areas of the Western Cape. Our results showed a lower impact risk with climate change for protected areas along the coast line, almost no change for northern protected areas and a higher impact risk in few isolated protected areas (Fig. 4 ). However, the impact risk differed according to the protection status. The impact risk remains high in areas of high protection status (such as national parks) and low for private protected areas (with lower protection status) (Fig. 5 ). Discussion This study presents, to the best of our knowledge, a first attempt at combining impacts of alien taxa with their projected distribution under climate change to produce risk maps for future climatic scenarios. Our results indicate that the projected suitable areas for alien Acacia species in South Africa are reduced under climate change (Fig. 2 ). This is similar to a study looking at A. mangium and A. auriculiformis in Brazil, where they found a shift in potential distribution and a reduction in suitable area (Heringer et al. 2019 ). More generally, this is in line with Bellard et al. ( 2018 ) who showed in a review of 71 papers covering 423 alien species that climate change is more frequently projected to contribute to a decrease in range size than an increase. However, to assume that the problem of Acacia invasions is going to sort itself out over time might be a bit too optimistic. For once, the sites at highest risk are also the regions where the highly vulnerable fynbos biome is located. This biome is already under pressure from the effects of climate change, including increased temperature and drought, and therefore more vulnerable to other pressures (Slingsby et al. 2017 ; Skowno et al. 2021 ). Furthermore, Acacia species which are not currently alien in the country (not yet introduced) were not modelled in this study. That means that other species not in our dataset might behave differently from what we present here. Therefore, this does not exclude the possibility of other Acacia species, and other alien species in general, becoming more problematic in future (see also Sheppard et al. 2016 ). Furthermore, Acacias and trees in general are long lived species which can result in a long lag between reduced climate suitability and decline in populations or observed reduction in impacts (Kowarik 1995 ; Robeck et al. 2024 ). Besides, perennial and predominantly allogame plants like Acacia species often exhibit lag phases of several decades before the beginning of an invasion (Robeck et al. 2024 ). Hence, several Acacia species could still be in a lag phase preceding an invasion in South Africa, and this invasion debt was not accounted for in the trends presented here (e.g., Rouget et al. 2016 ). As with any model, there are uncertainties linked to the projections. SDMs are based on several strong assumptions (Guisan et al. 2017 ; Hui 2023 ). First, occurrence records should reflect the true performance of the species. However, sampling biases are present in most if not all record databases (Beck et al. 2014 ). Moreover, low sample size can highly impact the performance of SDMs and data available may not be sufficient to fully inform the models (Wisz et al. 2008 ; Stockwell et al. 2002). To address these shortcomings, we designed a conservative methodology to keep only the most trustworthy records and limit oversampling in some areas. We also kept records from both the native and the alien range (Broennimann and Guisan 2008 ). Second, the species’ performance should respond directly to the variation of the selected predictors. In this study we chose four climate variables for their known link with plant species survival and development (Mod et al. 2016 ) and their availability in both fine spatial and temporal scales. Moreover, Sheppard et al. ( 2014 ), showed that predictions of similar SDMs of the response of three invasive plant species to climate change was highly correlated with field experiments. However, SDMs are correlative and may yield incorrect estimates of habitat suitability if climatic variables are correlated to other unknown variables in the training area (Guisan et al. 2017 ; Jarnevich et al. 2015 ). Future projections from our fitted models of climatic suitability could be biased by potential confounding factors of climate, as we do not take into account other factors which determine the success of alien species in a new region, such as dispersal capabilities and biotic interactions. Our model fitting procedure implicitly assumed that a species had the opportunity to colonise a large part of any country where it was recorded. Regarding interactions, Australian Acacia species often have competitive advantages over native plant species and tend to become dominant among plant communities, especially after disturbance (Morris et al. 2011 ), suggesting that the spatial extent of the realised niche would not be strongly restricted by competitive interactions compared to the potential niche. Yet, other interactions, such as mutualistic and trophic interactions with soil fungi (Birnbaum et al. 2018 ) may contribute to constrain the actual range of Acacias and act as confounding factors of climate in fitted SDMs, inducing bias in future projections. Furthermore, factors such as topography could be taken into account to improve the models (Bradley and Mustard 2006 ). Thus, our model predictions could be improved by considering non-climatic drivers (such as soil composition) if these variables were available at a fine resolution. Third, the species’ distribution, represented by recorded occurrences, should be stable and fill any available niche in the study environment. Several studies have shown violations of the niche conservatism hypothesis during invasion with niche shifts between the native niche and the introduced niche (e.g. Parravicini et al. 2015 ; Broennimann et al. 2007 ; Guisan et al. 2014 ). Moreover, predictions based on extrapolations on a new territory and with future climatic conditions may not be robust because the data used for model parameterization cannot represent all conditions in the extrapolated region (Elith and Leathwick 2009 ; Barbet-Massin et al. 2010 ; Sinclair et al. 2010 ). Thus, it must be kept in mind that our suitability and risk maps may be underestimated and should not be taken as a prediction of true future species richness and impacts. Fourth, some future climatic conditions may not have an analog amongst the historical climates of the study area (Williams & Jackson 2007 ). Yet, SDMs cannot predict exactly how species will respond to conditions that were not used for calibration (Pouteau et al. 2021 ). This could partly explain why most alien species (Bellard et al. 2018 ) including Acacias alien to South Africa, are projected to experience a decrease in the size of their potential range according to our current knowledge. Future work should consider the identification of novel climates so as to avoid putting too much confidence in climates with no current analogues. Caution is also advised when interpreting the risk maps including the sum of potential impacts of Acacias. Firstly, there are many ways to aggregate impacts, both, within species (calculating one impact value taking into account all impact records for the species) and across species (calculating an impact score for a site where several alien species are present) (Boulesnane-Genguant et al. in prep). Some of the most prominent methods to get one impact value per species have been to sum scores (e.g., Nentwig et al. 2016 ), to calculate a mean value (e.g., Rumlerova et al. 2016), and to take a maximum value (e.g., Blackburn et al. 2014 ) (see also Kumschick et al. 2024 ). They each come with underlying assumptions which can affect the results. The maximum score per species, as used here, was chosen for range shifting species in the US to anticipate if any high impacting invaders are likely to arrive under climate change (Rockwell-Postel et al. 2020 ). However, the applications of scores aggregated across species are scarce (for an example, see Nentwig et al. 2010 ). Furthermore, climate change can not only affect the potential distribution of species, but also modify their invasion behaviour and impacts (e.g., Le Maitre et al. 2020 ). Changes to fire regimes could have large impacts on alien species and native ecosystems, especially in the fynbos biome of South Africa where fires have been increasing in frequency and intensity due to alien invasions and climate change (Le Maitre et al. 2020 ; Slingsby et al. 2017 ). Furthermore, CO 2 concentration could lead to woody plant densification, which is already shown for native woody plants in some southern African habitats (Skowno et al. 2017 ). In the fynbos biome, increased CO 2 could also favour the alien wattles as they are nitrogen fixers as opposed to the native flora which are adapted to low nitrogen conditions (Richardson et al. 2014). In other examples, the synchronisation of the flowering period of native and alien plants may favour the latter, through increased interactions with pollinators to the detriment of the pollination of native species. This is the case in New Zealand with Calluna vulgaris whose greater phenological plasticity compared to the native species Dracophyllum subulatum means that it can reproduce more easily in areas with a high floral density (Giejsztowt et al. 2020 ). Morphological responses to climate change can also increase the competitiveness of alien species. For example, milder winter temperatures in China facilitate the survival of the water hyacinth Pontederia crassipes , and also allow it to develop a greater biomass, forming denser foliage that excludes submerged native plants the following season (You et al. 2013 ). In terrestrial environments, rising temperatures can also encourage the development of alien plants, reducing the availability of water for native plants, which are then at a disadvantage when it comes to coping with dry spells. This is the case of Tamarix spp. introduced in the United States, which develops greater capacity to capture and use water resources than native riparian species under the effect of drought (Hellmann et al. 2008 ). Given all these potential interactions between climate change and impacts of alien species, what we are showing here is not an accurate representation of sites with future impacts of alien Acacia species in South Africa. However, risk maps like these of sites where impacting species could occur under climate change projections can be valuable in helping us prioritise sites for future protection from invasion impacts. Despite the potential shortcomings of the models presented here, our study contributes to our understanding of the impacts of climate change on the risk of alien Acacia invasions, including their impacts. Furthermore, it can aid the prioritisation of clearing actions for alien Acacia species in South Africa and feed into strategies for protected area management and expansion (e.g., Department of Environmental Affairs 2016 ). Protected areas are set up to safeguard biodiversity and ecosystem services for the future. We show that for the Western Cape province, most protected areas should experience less impact risk in future. However, the areas with increased risk are of lowest protection, which makes them more vulnerable to negative changes and biodiversity loss. Elevating the protection of these sites and controlling harmful Acacias can help us reduce the loss of high biodiversity areas, as stipulated in the GBF Target 1. Understanding where some of the most damaging invaders might occur under climate change is important so our management is tackling not only current, but also potential future problems. Declarations Competing Interests The authors have no relevant financial or non-financial interests to disclose. Funding SK acknowledges the support of the Centre for Invasion Biology (CIB) at Stellenbosch University and the South African Department of Forestry, Fisheries and the Environment (DFFE). SK also acknowledges funding from the B3 project. B3 (Biodiversity Building Blocks for Policy) receives funding from the European Union’s Horizon Europe Research and Innovation Programme (through grant no .584 101059592). LJ was supported by the OUSNEXT project funded by the Structure fédérative de recherche BioST of the Université de La Réunion (PI: RP). OBG was supported by the Centre de Coopération Internationale en Recherche Agronomique pour le Développement (CIRAD), Fonds Vert and by the Université de La Réunion. This work was co-funded by the European Union: Agricultural Fund for Rural Development (EAFRD), by the Conseil Départemental de La Réunion, and by the Centre de Coopération Internationale en Recherche Agronomique pour le Développement (CIRAD). Author contributions SK and MR conceptualised and designed the study. CB provided distribution data. CB, LJ, OBJ and MR performed the data analyses. The first draft of the manuscript was written by SK and all authors contributed to the writing of the submitted draft and read and approved the final manuscript. Acknowledgements We thank Cally Jansen for data collection on impacts of alien Acacia species, and the participants of the workshop on “Managing invasions in protected areas” held in Stellenbosch from 3–5 April 2024 for valuable comments received on an earlier draft of the manuscript. SK acknowledges the support of the Centre for Invasion Biology (CIB) at Stellenbosch University and the South African Department of Forestry, Fisheries and the Environment (DFFE), noting that this manuscript does not necessarily reflect the opinions of the DFFE or its employees. SK also acknowledges funding from the B3 project. B3 (Biodiversity Building Blocks for Policy) receives funding from the European Union’s Horizon Europe Research and Innovation Programme (through grant no .584 101059592). The views and opinions expressed in the present article are, however, those of the authors only and do not necessarily reflect those of the European Union or the European Commission. Neither the EU nor the EC can be held responsible for them. LJ was supported by the OUSNEXT project funded by the Structure fédérative de recherche BioST of the Université de La Réunion (PI: RP). OBG was supported by the Centre de Coopération Internationale en Recherche Agronomique pour le Développement (CIRAD), Fonds Vert and by the Université de La Réunion. This work was co-funded by the European Union: Agricultural Fund for Rural Development (EAFRD), by the Conseil Départemental de La Réunion, and by the Centre de Coopération Internationale en Recherche Agronomique pour le Développement (CIRAD). References Barbet-Massin M, Rome Q, Villemant C, Courchamp F (2018) Can species distribution models really predict the expansion of invasive species? PLoS ONE 13(3):e0193085. https://doi.org/10.1371/journal.pone.0193085 Barbet-Massin M, Thuiller W, Jiguet F (2010) How much do we overestimate future local extinction rates when restricting the range of occurrence data in climate suitability models? Ecography 33:878–886 Barbet‐Massin M, Jiguet F, Albert CH, Thuiller W (2012) Selecting pseudo‐absences for species distribution models: How, where and how many?. Methods in Ecology and Evolution 3(2):327-338 Beck J, Böller M, Erhardt A, Schwanghart W (2014) Spatial bias in the GBIF database and its effect on modeling species’ geographic distributions. Ecological Informatics 19:10–15. https://doi.org/10.1016/j.ecoinf.2013.11.002 Bellard C, Thuiller W, Leroy B, Genovesi P, Bakkenes M, Courchamp F (2013) Will climate change promote future invasions?. Global Change Biology 19(12):3740-3748 Bellard C, Jeschke JM, Leroy B, Mace GM (2018) Insights from modeling studies on how climate change affects invasive alien species geography. Ecology and Evolution 8:5688–5700. https://doi.org/10.1002/ece3.4098 Birnbaum C, Morald TK, Tibbett M, Bennett RG, Standish RJ (2018) Effect of plant root symbionts on performance of native woody species in competition with an invasive grass in multispecies microcosms. Ecology and Evolution 8(17):8652-8664. Blackburn TM, Essl F, Evans T, Hulme PE, Jeschke JM, Kühn I, Kumschick S, Mrugała A, Marková Z, Nentwig W, Pergl J, Pyšek P, Rabitsch W, Ricciardi A, Richardson DM, Sendek A, Vilà M, Wilson JRU, Winter M, Genovesi P, Bacher S (2014) A unified classification of alien species based on the magnitude of their environmental impacts. PLoS Biology 12(5): e1001850. DOI: 10.1371/journal.pbio.1001850 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 Botella C, Marchante H, Celesti-Grapow L, Brundu G, Geerts S, Ramirez-Albores JE, González-Moreno P, Ritter M, Richardson DM (2023) The global distribution of Acacia. Wattles: Australian Acacia Species Around the World, 131. https://doi.org/10.1079/9781800622197.0009 Bradley BA, Wilcove DS, Oppenheimer M (2010) Climate change increases risk of plant invasion in the Eastern United States. Biological Invasions 12:1855–1872 Bradley BA, Mustard JF (2006) Characterizing the landscape dynamics of an invasive plant and risk of invasion using remote sensing. Ecological Applications 16(3):1132-1147 Broennimann O, Guisan A (2008) Predicting current and future biological invasions: both native and invaded ranges matter. Biol Lett 4:585–589 Broennimann O, Treier UA, Müller-Schärer H, Thuiller W, Peterson AT, Guisan A (2007) Evidence of climatic niche shift during biological invasion. Ecol Lett 10:701–709 Busby JR (1991) BIOCLIM-a bioclimate analysis and prediction system. Plant Protection Quarterly 6(1):8–9 Chamberlain S, Barve V, Mcglinn D, Oldoni D, Desmet P, Geffert L, Ram K (2022) rgbif: Interface to the Global Biodiversity Information Facility API. R package version 3.7.2. https://CRAN.R-project.org/package=rgbif Department of Environmental Affairs (2016) National protected areas expansion strategy for South Africa 2016. Department of Environmental Affairs, Pretoria, South Africa. Doak D, Morris W (2010) Demographic compensation and tipping points in climate-induced range shifts. Nature 467:959–962. doi:10.1038/nature09439 Elith J, Leathwick J (2009) Species distribution models: Ecological explanation and prediction across space and time. Evolution and Systematics 40:677–697 Giejsztowt J, Classen AT, Deslippe JR (2020) Climate change and invasion may synergistically affect native plant reproduction. Ecology 101(1):e02913 Guisan A, Petitpierre B, Broennimann, et al (2014) Unifying niche shift studies: insights from biological invasions. Trends in Ecology and Evolution 29:260–269. https://doi.org/10.1016/j.tree.2014.02.009 Guisan A, Thuiller W, Zimmermann NE (2017) Habitat suitability and distribution models: with applications in R. Cambridge University Press, Cambridge Guisan A, Zimmermann NE (2000) Predictive habitat distribution models in ecology. Ecological Modelling 135(2):147–186. https://doi.org/10.1016/S0304-3800(00)00354-9 Hawkins CL, Bacher S, Essl F, Hulme PE, Jeschke JM, Kühn I, Kumschick S, Nentwig W, Pergl J, Pyšek P, Rabitsch W, Richardson DM, Vilà M, Wilson JRU, Genovesi P, Blackburn TM (2015) Framework and guidelines for implementing the proposed IUCN Environmental Impact Classification for Alien Taxa (EICAT). Diversity and Distributions 21(11):1360-1363 Hellmann JJ, Byers JE, Bierwagen BG, Dukes JS (2008) Five potential consequences of climate change for invasive species. Conservation Biology 22(3):534-543 Heringer G, Bueno ML, Meira-Neto JA, Matos FA, Neri AV (2019) Can Acacia mangium and Acacia auriculiformis hinder restoration efforts in the Brazilian Atlantic Forest under current and future climate conditions?. Biological Invasions 21(9):2949-2962 Hui C (2023) The dos and don’ts for predicting invasion dynamics with species distribution models. Biological Invasions 25:947–953. https://doi.org/10.1007/s10530-022-02976-3 IPBES (2023) Summary for policymakers of the thematic assessment report on invasive alien species and their control of the Intergovernmental science-policy Platform on Biodiversity and Ecosystem Services. Roy HE, Pauchard A, Stoett P, Renard Truong T, Bacher S, Galil BS, Hulme PE, Ikeda T, Sankaran KV, McGeoch MA, Meyerson LA, Nuñez MA, Ordonez A, Rahlao SJ, Schwindt E, Seebens H, Sheppard AW, Vandvik V (eds.). IPBES secretariat, Bonn, Germany. https://doi.org/10.5281/zenodo.7430692 IUCN (2020a) IUCN EICAT Categories and Criteria. The Environmental Impact Classification for Alien Taxa, 1st edn. IUCN, Gland, Switzerland and Cambridge, UK. https://doi.org/10.2305/IUCN.CH.2020.05.en IUCN (2020b). Guidelines for using the IUCN Environmental Impact Classification for Alien Taxa (EICAT) Categories and Criteria: First edition. Version 1.1. Gland, Switzerland and Cambridge, UK. IUCN. Jansen C, Kumschick S (2022) A global impact assessment of Acacia species introduced to South Africa. Biological Invasions 24(1):175-187. https://doi.org/10.1007/s10530-021-02642-0 Jarnevich CS, Stohlgren TJ, Kumar S, Morisette JT, Holcombe TR (2015) Caveats for correlative species distribution modeling. Ecol Inform 29:6–15 Jiménez-Valverde A, Peterson AT, Soberón J, Overton JM, Aragón P, Lobo JM (2011) Use of niche models in invasive species risk assessments. Biol Invasions 13:2785–2797 Kaplan H, van Niekerk A, Le Roux JJ, Richardson DM, Wilson JR (2014) Incorporating risk mapping at multiple spatial scales into eradication management plans. Biological Invasions 16:691-703 Karger DN, Conrad O, Böhner J, Kawohl T, Kreft H, Soria-Auza RW, Zimmermann NE, Linder HP, Kessler M (2017) Climatologies at high resolution for the earth’s land surface areas. Scientific Data 4(1):1-20 Kowarik I (1995) Time lags in biological invasions with regard to the success and failure of alien species. Plant invasions: general aspects and special problems, 15-38. Kumschick S, Jansen C (2023) Evidence-based impact assessment for naturalized and invasive Australian Acacia species. Wattles: Australian Acacia species around the world (eds. D.M. Richardson, J.J. Le Roux & E.M. Marchante). CABI, Wallingford, pp. 359-381. Kumschick S, Bertolino S, Blackburn TM, Brundu G, Costello KE, de Groot M, Evans T, Gallardo B, Genovesi P, Govender T, Jeschke JM, Lapin K, Measey GJ, Novoa A, Nunse AL, Probert AF, Pyšek P, Preda C, Rabitsch W, Roy HE, Smith KG, Tricarico E, Vilà M, Vimercati G, Bacher S (2024) Using the IUCN Environmental Impact Classification for Alien Taxa (EICAT) to inform decision-making. Conservation Biology 38. DOI:10.1111/cobi.14214 Le Maitre DC, Blignaut JN, Clulow A, Dzikiti S, Everson CS, Görgens AH, Gush MB (2020) Impacts of plant invasions on terrestrial water flows in South Africa. In Biological Invasions in South Africa (pp. 431-457). Cham: Springer International Publishing. Le Maitre DC, Gaertner M, Marchante E et al. (2011) Impacts of invasive Australian acacias: implications for management and restoration. Diversity and Distributions 17:1015–1029 Lenoir J et al. (2008) A significant upward shift in plant species optimum elevation during the 20th century. Science 320:1768–1771 Liaw A, Wiener M (2002) Classification and regression by randomForest. R news 2(3):18-22 Magona N, Richardson DM, le Roux JJ, Kritzinger-Klopper S, Wilson JRU (2018) Even well-studied groups of alien species might be poorly inventoried: Australian Acacia species in South Africa as a case study. NeoBiota 39:1-29. https://doi.org/10.3897/neobiota.39.23135 Manes S et al. (2021) Endemism increases species’ climate change risk in areas of global biodiversity importance. Biological Conservation 257. https://doi.org/10.1016/j.biocon.2021.109070. Marchioro CA, Krechemer FS (2021) Prevention is better than cure: Integrating habitat suitability and invasion threat to assess global biological invasion risk by insect pests under climate change. Pest Management Science 77(10):4510-4520 Mod HK, Scherrer D, Luoto M, Guisan A (2016) What we use is not what we know: environmental predictors in plant distribution models. Journal of Vegetation Science 27(6):1308-1322 Morris TL, Esler KJ, Barger NN, Jacobs SM, Cramer MD (2011) Ecophysiological traits associated with the competitive ability of invasive Australian acacias. Diversity and Distributions 17(5):898-910 Nentwig W, Bacher S, Pyšek P, Vilà M, Kumschick S (2016) The Generic Impact Scoring System (GISS): a standardized tool to quantify the impacts of alien species. Environmental Monitoring and Assessment 188:315. doi:10.1371/journal.pbio.1001850. Nentwig W, Kühnel E, Bacher S (2010) A generic impact‐scoring system applied to alien mammals in Europe. Conservation Biology 24(1):302-311 Parravicini V, Azzurro E, Kulbicki M, Belmaker J (2015) Niche shift can impair the ability to predict invasion risk in the marine realm: an illustration using Mediterranean fish invaders. Ecol Lett 18:246–253 Pouteau R, Biurrun I, Brunel C, Chytrý M, Dawson W, Essl F, Fristoe T, Haveman R, Hobohm C, Jansen F, Kreft H, Lenoir J, Lenzner B, Meyer C, Moeslund JE, Pergl J, Pyšek P, Svenning J-C, Thuiller W, Weigelt P, Wohlgemuth T, Yang Q, van Kleunen M (2021) Potential alien ranges of European plants will shrink in the future, but less so for already naturalized than for not yet naturalized species. Diversity and Distributions 27:2063–2076. https://doi.org/10.1111/ddi.13378 Richardson DM, Allsopp N, D’Antonio CM et al. (1994) Plant invasions – the role of mutualisms. Biol Rev 75:65–93. https://doi.org/10.1111/j.1469-185X.1999.tb00041.x Robeck P, Essl F, van Kleunen M, Pyšek P, Pergl J, Weigelt P, Mesgaran MB (2024) Invading plants remain undetected in a lag phase while they explore suitable climates. Nature Ecology & Evolution, 1-12. Rockwell-Postel M, Laginhas BB, Bradley BA (2020) Supporting proactive management in the context of climate change: prioritizing range-shifting invasive plants based on impact. Biological Invasions 22:2371-2383 Rouget M, Robertson MP, Wilson JRU, Hui C, Essl F, Renteria JL, Richardson DM (2016) Invasion debt - quantifying future biological invasions. Diversity and Distributions 22(4):445-456 Rumlerová Z, Vilà M, Pergl J, Nentwig W, Pyšek P (2016) Scoring environmental and socioeconomic impacts of alien plants invasive in Europe. Biological Invasions 18:3697-3711 Sheppard CS, Burns BR, Stanley MC (2014) Predicting plant invasions under climate change: are species distribution models validated by field trials? Glob Change Biol 20:2800–2814 Sheppard CS, Burns BR, Stanley MC (2016) Future-proofing weed management for the effects of climate change: is New Zealand underestimating the risk of increased plant invasions?. New Zealand Journal of Ecology 40(3):398-405 Sinclair SJ, White MD, Newell GR (2010) How useful are species distribution models for managing biodiversity under future climates? Ecology and Society 15(1):8. Skowno AL, Thompson MW, Hiestermann J et al (2017) Woodland expansion in South African grassy biomes based on satellite observations (1990-2013): general patterns and potential drivers. Glob Chang Biol 23:2358–2369. https://doi.org/10.1111/gcb.13529 Skowno AL, Jewitt D, Slingsby JA (2021) Rates and patterns of habitat loss across South Africa's vegetation biomes. South African Journal of Science 117(1-2):1-5 Slingsby JA, Merow C, Aiello-Lammens M, Allsopp N, Hall S, Kilroy Mollmann H, Turner R, Wilson AM, Silander Jr JA (2017) Intensifying postfire weather and biological invasion drive species loss in a Mediterranean-type biodiversity hotspot. Proceedings of the National Academy of Sciences 114(18):4697-4702 Smith K (2020) The IUCN Red List and invasive alien species: An analysis of impacts on threatened species and extinctions. International Union for the Conservation of Nature. Stockwell DRB, Peterson AT (2002) Effects of sample size on accuracy of species distribution models. Ecol Modell 148:1–13. Thomas CD, Cameron A, Green RE, Bakkenes M, Beaumont LJ, Collingham YC, Erasmus BF, De Siqueira MF, Grainger A, Hannah L, Hughes L, et al. (2004) Extinction risk from climate change. Nature 427(6970):145-148 Thuiller W, Lafourcade B, Engler R, Araújo MB (2009) BIOMOD–a platform for ensemble forecasting of species distributions. Ecography 32(3):369-373 Tingley R, García-Díaz P, Arantes CRR, Cassey P (2017) Integrating transport pressure data and species distribution models to estimate invasion risk for alien stowaways. Ecography Vilà M, Espinar JL, Hejda M, Hulme PE, Jarošík V, Maron JL, Pergl J, Schaffner U, Sun Y, Pyšek P (2011) Ecological impacts of invasive alien plants: a meta‐analysis of their effects on species, communities and ecosystems. Ecology Letters 14(7):702-708 Vilà M, Trillo A, Castro-Díez P, Gallardo B, Bacher S (2024)  Field studies of the ecological impacts of invasive plants in Europe. NeoBiota 90:139-159 Volery L, Blackburn TM, Bertolino S, Evans T, Genovesi P, Kumschick S, Roy HE, Smith KG, Bacher S (2020) Improving the Environmental Impact Classification for Alien Taxa (EICAT): a summary of revisions to the framework and guidelines. NeoBiota 62:547–567 Wenger SJ, Olden JD (2012) Assessing transferability of ecological models: an underappreciated aspect of statistical validation. Methods in Ecology and Evolution 3(2):260-267 Williams JW, Jackson ST (2007) Novel climates, no-analog communities, and ecological surprises. Frontiers in Ecology and the Environment 5(9):475–482. https://doi.org/10.1890/070037 Wisz MS, Hijmans RJ, Li J, Peterson AT, Graham CH, Guisan A (2008) Predicting Species Distributions Working Group. Effects of sample size on the performance of species distribution models. Diversity and Distributions 14:763–773 You W, Yu D, Xie D, Yu L (2013) Overwintering survival and regrowth of the invasive plant Eichhornia crassipes are enhanced by experimental warming in winter. Aquatic Biology 19:45-53 Supplementary Files SupplementaryMaterial.docx appS2finalImpactmapsformulas.pdf Cite Share Download PDF Status: Published Journal Publication published 18 Dec, 2024 Read the published version in Biological Invasions → Version 1 posted Reviewers agreed at journal 13 Jun, 2024 Reviewers invited by journal 26 May, 2024 Editor assigned by journal 21 May, 2024 First submitted to journal 17 May, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4437291","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":306902120,"identity":"873f8fd6-4243-4f6f-9562-7681a8868a32","order_by":0,"name":"Sabrina Kumschick","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABCElEQVRIie2RsUoDQRCG5zi4aiHtBsF7hQkLaxP2We5YMM0pqWwUTHU2IdvmScRy4SBpTuxkIU2usQohEBACIs4FtNpgYmWxX7MDy8f/DwMQCPxLGEAGgO0Y04Od05Xu6BgFfhQa0P4ipA/Ps03zpAS+zpvtcKiEeKnlEu4UffllrK/1NK+1xEUhzqaopXTFBcJMQzTO/AoUAvIy7pOSxAxtXzomOSSWavqV1Kxa5Z6UwdteEaYm5dNC0ln6l3H7lIqKZbJVJOVKHpUWGD9QzK01KXPRbXdhqAV3lzc8n2jGuT8lNVdVtCtve5PFoNmyD9Uzpnrkm3d1nhp/ygGy73sFAoFA4E98AdcgVMQeccraAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0001-8034-5831","institution":"Stellenbosch University Faculty of Science","correspondingAuthor":true,"prefix":"","firstName":"Sabrina","middleName":"","lastName":"Kumschick","suffix":""},{"id":306902121,"identity":"48e648d4-c8ac-4365-bfaf-122a2d423447","order_by":1,"name":"Lysandre Journiac","email":"","orcid":"","institution":"Muséum National d'Histoire Naturelle: Museum National d'Histoire Naturelle","correspondingAuthor":false,"prefix":"","firstName":"Lysandre","middleName":"","lastName":"Journiac","suffix":""},{"id":306902122,"identity":"ea747222-a802-42fe-87f3-e208c2476469","order_by":2,"name":"Océane Boulesnane-Genguant","email":"","orcid":"","institution":"Reunion University: Universite de la Reunion","correspondingAuthor":false,"prefix":"","firstName":"Océane","middleName":"","lastName":"Boulesnane-Genguant","suffix":""},{"id":306902123,"identity":"9dd1473d-5f7c-44e8-af6b-81a193d7c37e","order_by":3,"name":"Christophe Botella","email":"","orcid":"","institution":"Montpellier University: Universite de Montpellier","correspondingAuthor":false,"prefix":"","firstName":"Christophe","middleName":"","lastName":"Botella","suffix":""},{"id":306902124,"identity":"1a214a39-253b-4cc2-a90e-d1b847e8776d","order_by":4,"name":"Robin Pouteau","email":"","orcid":"","institution":"AMAP: Botanique et Modelisation de l'Architecture des Plantes et des Vegetations","correspondingAuthor":false,"prefix":"","firstName":"Robin","middleName":"","lastName":"Pouteau","suffix":""},{"id":306902125,"identity":"1de06cd3-e7fe-402e-bb2e-ced3fa50fa90","order_by":5,"name":"Mathieu Rouget","email":"","orcid":"","institution":"Centre de Coopération Internationale en Recherche Agricole pour le Développement Systèmes biologiques: CIRAD Departement Systemes biologiques","correspondingAuthor":false,"prefix":"","firstName":"Mathieu","middleName":"","lastName":"Rouget","suffix":""}],"badges":[],"createdAt":"2024-05-17 14:17:44","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4437291/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4437291/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s10530-024-03490-4","type":"published","date":"2024-12-18T15:57:42+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":58151657,"identity":"03006192-d2c5-48aa-befc-c6cae3f7d65d","added_by":"auto","created_at":"2024-06-11 20:12:23","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":110704,"visible":true,"origin":"","legend":"\u003cp\u003eStep by step explanation of the development of the models and maps for future climate projections. Five climate models were selected, and three socio-economic pathway scenarios (ssp) modelled. For details on the models and variables selected, see text.\u003c/p\u003e","description":"","filename":"1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4437291/v1/5afaa53374fb20476114fa40.jpg"},{"id":58151658,"identity":"76b8bc45-fd3d-4ed8-9b7d-2dfa921f9a1e","added_by":"auto","created_at":"2024-06-11 20:12:23","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":100656,"visible":true,"origin":"","legend":"\u003cp\u003eCumulative probability of presence ranging from zero (no species projected to establish) to thirty (all species with the highest probability to establish) of alien acacias (left) and projected impacts (right) based on the maximum recorded impacts of the species present over three time horizons: Top row: Current climate (1981-2010), middle row: mid-century (2041-2070), bottom row: end of century (2071-2100). The cumulative probability of presence on the left is the sum of single species establishment probability. Each species probability is a mean across combinations of global climate models and socio-economic scenarios. The right hand column shows the sum of species maximum impact weighted by their establishment probability.\u003c/p\u003e","description":"","filename":"2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4437291/v1/9ee44d2408417f2d5b4e95f3.jpg"},{"id":58151659,"identity":"c282d84f-2c86-48c2-afae-a4f261713bed","added_by":"auto","created_at":"2024-06-11 20:12:23","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":75234,"visible":true,"origin":"","legend":"\u003cp\u003eSum of maximum recorded impacts across acacia species based on EICAT records in Jansen et al. (2023) occurring per Quarter Degree Square (QDS; 25 km x 25 km) grid cell. This map is based on actual occurrence records of acacia species in South Africa taken from the database of Botella et al. (2023), with sum of impacts overlayed. Contrary to the projected maps in Figure 2, no records were excluded (including urban areas). The darker the red color, the higher the impact sum of the species occurring in these cells.\u003c/p\u003e","description":"","filename":"3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4437291/v1/782e019e2a998fdfbbaaac95.jpg"},{"id":58151662,"identity":"72d875d2-ddc1-4cdf-93f2-955d6df138a5","added_by":"auto","created_at":"2024-06-11 20:12:23","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":97616,"visible":true,"origin":"","legend":"\u003cp\u003eProjected changes in impact risk of Acacia species under future climate for protected areas of the Western Cape province. For each protected area, we compared the mean impact risk under current climate and end-of-century (2071-2100) climate. Negative values (in blues) indicate a reduction in impact risk, while positive values (in red) indicate an increase in impact risk.\u003c/p\u003e","description":"","filename":"4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4437291/v1/f97a6dc62bee713f70b7916e.jpg"},{"id":58151663,"identity":"b0845a83-2684-4c20-9140-b9cb57c6dfac","added_by":"auto","created_at":"2024-06-11 20:12:23","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":58545,"visible":true,"origin":"","legend":"\u003cp\u003eChanges in impact risk of Acacia species in the Western Cape province of South Africa according to protected areas status. Protected areas are listed in increasing order of protection status (from left to right). For each type, impact risk of Acacia species have been summarised under current and end-of-century (2071-2100) climate.\u003c/p\u003e","description":"","filename":"5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4437291/v1/c6a381d0c9da4c046bb5e07d.jpg"},{"id":72201846,"identity":"60bbe405-62f0-4d21-9dee-629731281c39","added_by":"auto","created_at":"2024-12-23 16:11:06","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":868009,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4437291/v1/e217404f-0f78-46f9-8a20-9637744c696f.pdf"},{"id":58151660,"identity":"367bd6c3-3c6f-4187-b423-ae58cf8796fc","added_by":"auto","created_at":"2024-06-11 20:12:23","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":757336,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryMaterial.docx","url":"https://assets-eu.researchsquare.com/files/rs-4437291/v1/196e83b65a033e49eff1f37e.docx"},{"id":58151661,"identity":"c836e989-3d3e-48ac-9c67-a334d8bbb03e","added_by":"auto","created_at":"2024-06-11 20:12:23","extension":"pdf","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":79029,"visible":true,"origin":"","legend":"","description":"","filename":"appS2finalImpactmapsformulas.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4437291/v1/7df74cc2c9737d3bb9def9c8.pdf"}],"financialInterests":"","formattedTitle":"Mapping impacts of alien species on biodiversity in the face of climate change","fulltext":[{"header":"Introduction","content":"\u003cp\u003eBiological invasions and climate change are two of the main drivers of global biodiversity loss (Master and Norgrove 2010; Johnson et al. 2017; IPBES \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Impacts of biological invasions are manifold, including changes to ecosystem functions such as fire regimes and hydrology and changes to biodiversity including extirpations and extinctions of native species (Vil\u0026agrave; et al. \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Smith \u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Vil\u0026agrave; et al. \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2024\u003c/span\u003e; IPBES \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Similarly, climate change can have severe impacts on native species, for example increasing their risk of extinctions (e.g., Thomas et al. \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e2004\u003c/span\u003e; Manes et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) and causing range changes (e.g., Lenoir et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Doak et al. 2010). While these drivers can act in isolation, climate change can also affect the invasion of alien species as climatic conditions play an important role in the growth, survival and distribution of species (Walther et al. 2009; Master and Norgrove 2010). Climate change can not only lead to changes in native species distributions, but also those of alien species (Bellard et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2013\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eMember states of the Convention on Biological Diversity (CBD) are mandated to improve the protection of natural resources under the Kunming-Montreal Global Biodiversity Framework (GBF). Target 1 asks to minimize the loss of areas with high biodiversity. To achieve this, we need to understand the pressures on these areas from drivers such as climate change and biological invasions, and map them to prioritise control. Risk maps have been used for various purposes in invasion science. For example, they can be useful to improve detection of new incursions or individuals of species with limited populations (Kaplan et al. \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). They have also been used to map the potential impacts to assess sites experiencing the highest current or potential impact, translating to sites at risk of invasion impacts (Nentwig et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). However, as species\u0026rsquo; distributions are expected to change with climate change, considering only current climatic conditions in such risk maps does not prepare us well for future risk scenarios. Therefore, to improve our ability to manage biological invasions, as stipulated under Target 6 of the GBF, we need to understand how biological invasions and their impacts could change under climate change. Simultaneously, this addresses Target 8, which aims to minimize the effect of climate change on biodiversity. Still, climate change is not routinely included in risk assessments for alien taxa (but see Marchioro and Krechemer \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2021\u003c/span\u003e for some insect pests).\u003c/p\u003e \u003cp\u003eCorrelative species distribution models (hereafter SDM) are popular in various fields of ecology and are increasingly used in invasion science (Guisan et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). They are statistical methods which correlate known occurrences of species with environmental variables and predict a species' potential distribution in other regions over space and time (Guisan and Zimmermann \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2000\u003c/span\u003e). The advantages of SDMs are manifold as they are easy to implement and can cover large geographic areas. Therefore, they are widely used to produce maps of potential invasion by delineating areas based on climatic suitability for a species (e.g. Bradley et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Jim\u0026eacute;nez-Valverde et al. \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Tingley et al. \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Despite known limitations (see for example Hui \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), they have shown potential for accurate prediction of alien species spread (Barbet-Massin et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) and response of biological invasions to climate change (Sheppard et al. \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) if properly calibrated.\u003c/p\u003e \u003cp\u003eIn this study, we use Australian Acacia species alien to South Africa as a case study to assess i) the current distribution of the species and their potential impacts, ii) the potential distribution under current climate, and iii) the projected future distribution under climate change scenarios. By overlaying the outputs of the models with the species potential impacts, we produce risk maps for current and projected future climatic conditions and assess potential impacts for protected areas of the Western Cape Province in South Africa.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003eWe detail below how we selected Australian Acacia species and fitted SDM models per species that were then used to project the habitat suitability across South Africa per time horizon, accounting for various climatic scenarios. Then, we explain how we aggregated an impact score from the species-wise SDM projections and the documented impacts per species.\u003c/p\u003e\n\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\n \u003ch2\u003eSpecies selection\u003c/h2\u003e\n \u003cp\u003eWe selected, as a case study, Australian Acacia species (\u003cem\u003eAcacia\u003c/em\u003e Mill. \u003cem\u003esensu stricto\u003c/em\u003e \u0026ndash; synonym: \u003cem\u003eAcacia\u003c/em\u003e subg. \u003cem\u003ePhyllodineae\u003c/em\u003e (DC.) Seringe) which are alien in South Africa. This includes 33 species (Magona et al. \u003cspan class=\"CitationRef\"\u003e2018\u003c/span\u003e; Jansen and Kumschick \u003cspan class=\"CitationRef\"\u003e2022\u003c/span\u003e; Kumschick and Jansen \u003cspan class=\"CitationRef\"\u003e2023\u003c/span\u003e; Supplementary Material Appendix S4.1), of which three species (\u003cem\u003eAcacia crassiuscula\u003c/em\u003e, \u003cem\u003eA. acuminata\u003c/em\u003e and \u003cem\u003eA. koa\u003c/em\u003e) had to be excluded from all analysis because they had less than 10 presence records after the filtering steps detailed below. Acacias are among the most highly impacting alien taxa in South Africa with a diverse range of impacts on native biodiversity (le Ma\u0026icirc;tre et al. 2011; Kumschick \u0026amp; Jansen \u003cspan class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec4\" class=\"Section2\"\u003e\n \u003ch2\u003eOccurrence data\u003c/h2\u003e\n \u003cp\u003eWe gathered distribution data from the native and introduced ranges of the 30 Acacia species, and used it along with bioclimatic variables to build individual ensemble models which were then projected based on current and future climate in South Africa. For each species, worldwide occurrence data were obtained from the Global Biodiversity Information Facility (GBIF, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.gbif.org/\u003c/span\u003e\u003c/span\u003e) using the rgbif package (Chamberlain et al. \u003cspan class=\"CitationRef\"\u003e2022\u003c/span\u003e). For this extraction, we kept all records which had complete coordinates, which were tagged as \u0026ldquo;presence\u0026rdquo;, and which were either human or machine observations or coming from the literature. No time limit was specified. Then we filtered occurrences to reduce the uncertainty arising from questionable records and only kept records which came from a list of 65 trusted GBIF datasets. This list was previously developed for the global distribution of Acacia species (Botella et al. \u003cspan class=\"CitationRef\"\u003e2023\u003c/span\u003e) and is composed of datasets which include at least one documented step of taxonomic verification in the recording of the occurrences. The documented native and introduced countries of each species were also obtained from that study. Besides, we extracted the global land cover raster (approx. 300 m resolution) from the European Space Agency GlobCover product. We excluded records located outside the native and introduced countries of the species (obtained from Botella et al. \u003cspan class=\"CitationRef\"\u003e2023\u003c/span\u003e) or whose 4 km x 4 km cell was located in urban areas based on the GlobCover 2009 global land cover raster (ESA and UCLouvain GlobCover 2009 Project, last accessed 29/01/2024). The latter step was implemented to exclude planted trees. More precisely, we excluded cells classified as \u0026ldquo;Artificial surfaces and associated areas (Urban areas\u0026thinsp;\u0026gt;\u0026thinsp;50%)\u0026rdquo; by GlobCover. We then resampled our remaining records to keep only up to one record per 4 km \u0026times; 4 km cell for each species to reduce the oversampling of some areas (Boria et al. \u003cspan class=\"CitationRef\"\u003e2014\u003c/span\u003e). We obtained a total of 37,940 filtered records (Figure \u003cspan class=\"InternalRef\"\u003eS4\u003c/span\u003e.1 \u0026amp; S4.2). The number of filtered records were imbalanced across the 30 considered species, with 27 for the least sampled (\u003cem\u003eA. adunca\u003c/em\u003e) and 3,738 for the most sampled (\u003cem\u003eA. dealbata\u003c/em\u003e).\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec5\" class=\"Section2\"\u003e\n \u003ch2\u003eClimate variables\u003c/h2\u003e\n \u003cp\u003eTo estimate the climatic suitability for each Acacia species from their occurrence data, we extracted four bioclimatic variables from the Climatologies at High Resolution for Earth\u0026rsquo;s Land Surface Areas database (CHELSA, Karger et al. \u003cspan class=\"CitationRef\"\u003e2017\u003c/span\u003e). We chose the mean annual air temperature (bio1), mean diurnal air temperature range (bio2), annual precipitation (bio12) and precipitation seasonality (bio15) as environmental variables as they are likely to affect plant productivity and survival across the landscape (Mod et al., \u003cspan class=\"CitationRef\"\u003e2016\u003c/span\u003e), and the correlation between each of these variables was moderate (Pearson\u0026rsquo;s |\u003cem\u003er\u003c/em\u003e| \u0026lt; 0.55). These data are available globally at a 4 km \u0026times; 4 km resolution for the period 1981\u0026ndash;2010.\u003c/p\u003e\n \u003cp\u003eThen, we predicted the suitable climatic range of each species per time period in South Africa. The \u0026ldquo;current\u0026rdquo; South African climate was derived from the 1981\u0026ndash;2010 period and the \u0026ldquo;future\u0026rdquo; climate was decomposed in two horizons: mid-century (2041\u0026ndash;2070) and end of the century (2071\u0026ndash;2100). For both time horizons, we extracted the four climatic variables based five Global Climate Models from CHELSA (\u0026ldquo;gfdl-esm4\u0026rdquo;, \u0026ldquo;ipsl-cm6a\u0026rdquo;, \u0026ldquo;mpi-esm1\u0026rdquo;, \u0026ldquo;mri-esm2\u0026rdquo;, \u0026ldquo;ukesm1\u0026rdquo;) for three shared socio-economic pathway scenarios (\u0026ldquo;ssp1.26\u0026rdquo;, \u0026ldquo;ssp3.70\u0026rdquo; and \u0026ldquo;ssp5.85\u0026rdquo;) for the South African region: ssp1.26: scenario with low greenhouse gas (GHG) emission; ssp3.70: scenario with high GHG emission; ssp5.85: scenario with very high GHG emission (Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e and Supplementary Material, Appendix S1).\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec6\" class=\"Section2\"\u003e\n \u003ch2\u003eSpecies distribution models\u003c/h2\u003e\n \u003cp\u003ePseudo-absences selection\u003c/p\u003e\n \u003cp\u003eFor each species we drew ten datasets of pseudo-absences composed of the same number of pseudo-absences as presences. For species with less than 100 presences, we increased the number of pseudo-absences to 100 in each dataset, to reduce the variance (over-fitting) of each model. These pseudo-absences were drawn in the countries being part of the native and introduced ranges of the species assuming that the species likely had an opportunity to establish at pseudo-absence locations. Following recommendations of Barbet-Massin et al. (\u003cspan class=\"CitationRef\"\u003e2012\u003c/span\u003e), pseudo-absences were selected based on the surface range envelope method (Busby \u003cspan class=\"CitationRef\"\u003e1991\u003c/span\u003e), by randomly selecting pseudo-absences outside the usual climatic range of the species occurrences (Thuiller et al. \u003cspan class=\"CitationRef\"\u003e2009\u003c/span\u003e). Given that pseudo-absences were drawn in the native and introduced ranges but outside of the climatic envelope where the species occurrences were recorded (PA.sre.quant\u0026thinsp;=\u0026thinsp;0), we can assume false absences to be rare. Therefore, the final model predictions will approximate the probability of presence conditional to the species presence in the area, i.e. an establishment probability.\u003c/p\u003e\n \u003cp\u003eModel implementation\u003c/p\u003e\n \u003cp\u003eFor each species we built a species distribution model using functions and workflow of the biomod2 package (Thuiller et al. \u003cspan class=\"CitationRef\"\u003e2009\u003c/span\u003e). We applied a random forest classifier (Liaw and Wiener \u003cspan class=\"CitationRef\"\u003e2002\u003c/span\u003e) on the presence/pseudo-absence datasets and the four climatic variables to decipher which climatic conditions are suitable for each species. Each forest was composed of 500 trees, where each branch tests two randomly picked variables (mtry\u0026thinsp;=\u0026thinsp;2), and we impose a minimum of five sampled presence/pseudo-absences per terminal node. We implemented a step of cross-validation for which the presence/pseudo-absence dataset was divided into four blocks along the longitude axis with presences equally balanced between each block (Wenger and Olden \u003cspan class=\"CitationRef\"\u003e2012\u003c/span\u003e). Then, three blocks were used for model training and one for testing. This approach allows for spatial and environmental independence between training and testing sets. Kappa and TSS metrics were computed each time on the testing dataset. As we had ten pseudo-absence datasets and four blocks, we computed 40 random forest models for each species.\u003c/p\u003e\n \u003cp\u003eFor each species we kept only the random forest models with a Kappa and a TSS over 0.7, indicating a robust predictive capacity, and which we refer to as sub-models. We then derived for each species an ensemble of its sub-models that was a \u0026ldquo;committee averaging\u0026rdquo; of the binarized responses given by the sub-models (Thuiller et al. \u003cspan class=\"CitationRef\"\u003e2009\u003c/span\u003e). We used a probability threshold of 0.5 to binarize each submodel response, corresponding to a majority vote among trees of the corresponding random forest. Committee averaged values are the proportion of sub-models predicting the species as present. Here, we considered that committee averaged values can be regarded as a proxy for the probability of a species to establish. Thus, we obtained one projection of establishment probability for South Africa under the current climate and 15 scenarios (five global climate models x three socio-economic pathways) for each time horizon (2041\u0026ndash;2070 and 2071\u0026ndash;2100) as represented in Fig.\u0026nbsp;\u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e. So, for each time horizon and each spatial cell, we averaged the 15 establishment probabilities of the 15 scenarios to obtain the final establishment probability of the species in that cell. This final establishment probability gives an equal weight to the 15 plausible scenarios. Finally, for each time horizon, we produced a potential Acacia species richness map summing the 30 species establishment probability maps, and an associated uncertainty map. The uncertainty map shows standard deviation of the potential richness estimate based on the variability standard deviation of the predicted establishment per species and its variability across the 15 scenarios.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec7\" class=\"Section2\"\u003e\n \u003ch2\u003eRisk maps\u003c/h2\u003e\n \u003cp\u003eAs we aimed to map the accumulated risk of Acacia species potentially causing impacts in different regions in South Africa, we firstly derived a potential impact score per species based on its documented impacts. Data on the impacts of the Acacia species were taken from Jansen and Kumschick (\u003cspan class=\"CitationRef\"\u003e2022\u003c/span\u003e). They assigned categories of Minimal Concern to Massive to each species, based on the impacts reported in the literature and following the International Union for the Conservation of Nature (IUCN) standard for impact classification, the Environmental Impact Classification for Alien Taxa, EICAT (Blackburn et al. \u003cspan class=\"CitationRef\"\u003e2014\u003c/span\u003e; Hawkins et al. \u003cspan class=\"CitationRef\"\u003e2015\u003c/span\u003e; IUCN \u003cspan class=\"CitationRef\"\u003e2020a\u003c/span\u003e, \u003cspan class=\"CitationRef\"\u003eb\u003c/span\u003e; Volery et al. \u003cspan class=\"CitationRef\"\u003e2020\u003c/span\u003e). We used the maximum impact category per species as a precautionary measure, which is the recommended approach used by the IUCN EICAT Standard (IUCN \u003cspan class=\"CitationRef\"\u003e2020a\u003c/span\u003e; Kumschick et al. \u003cspan class=\"CitationRef\"\u003e2024\u003c/span\u003e), noting that other calculation methods are possible depending on the aim of the study (see also Kumschick et al. \u003cspan class=\"CitationRef\"\u003e2024\u003c/span\u003e, Boulesnane-Genguant et al. in prep). We transformed the impact categories into numerical values as follows: Minimal Concern (MC)\u0026thinsp;=\u0026thinsp;0, Minor (MN)\u0026thinsp;=\u0026thinsp;1, Moderate (MO)\u0026thinsp;=\u0026thinsp;2, Major (MR)\u0026thinsp;=\u0026thinsp;3, Massive (MV)\u0026thinsp;=\u0026thinsp;4. 19 species classified as Data Deficient (DD), i.e. no data on impact was found in a comprehensive literature search, and they were assigned a score of 0 (Supplementary Material Appendix S4.1).\u003c/p\u003e\n \u003cp\u003eAs a prediction of local impact risk per species, we multiplied the predicted establishment probability derived from the SDMs in each grid cell by the species impact score. Finally, we summed up the local impact risks across species in each grid cell to obtain the impact risk map (see also Nentwig et al. \u003cspan class=\"CitationRef\"\u003e2010\u003c/span\u003e). Hence, each species occurring adds to the risk of the potential impacts actually manifesting at that site, again using a precautionary approach. Our maps therefore do not represent sites where impacts are definitely going to occur, as the potential establishment of each species alone does not necessarily happen and, even if it does, it does not necessarily result in the worst impact documented for that species. Other factors such as the recipient community and the abundance of the alien species play important roles in the manifestation of impacts. Note also that even if data deficient species are present, they do not affect the impact risk value, as for species whose maximum impact is MC. This might underestimate the impact of some species which actually cause higher impacts but which have not been studied.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\n \u003ch2\u003eChanges in impact risk in protected areas\u003c/h2\u003e\n \u003cp\u003eTo assess the specific change in risk across time on protected areas, we assessed the potential changes of impact risk for protected areas in the Western Cape Province of South Africa. For each protected area, we calculated the mean impact risk score and compared it under current and future (2071\u0026ndash;2100) projected climate. To see whether differences could be found based on the protected area status, we also calculated the mean impact risk score per protected area type. In the Western Cape, protected areas are classified as follows (with increasing level of protection): private nature reserves, provincial nature reserves and national parks. We obtained protected areas data from the South African National Biodiversity Institute (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ewww.bgis.sanbi.org\u003c/span\u003e\u003c/span\u003e).\u003c/p\u003e\n\u003c/div\u003e"},{"header":"Results","content":"\u003cp\u003eFirstly, we present a method to map the potential cumulative impacts of alien species under climate change projections (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, Supplementary Material, Appendix S1). It combines the recorded impacts alien species have had with the projected species richness. Applying this method to Acacia species in South Africa, we find that across all species currently present in the country, the suitable area declines over time and only the south coast in the Eastern Cape province and parts of the Western Cape province still remain suitable for a considerable number of species (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Large parts of the north eastern (NE) part of the country are predicted to become unsuitable for most wattles in future. Indeed, while the current projected impact score lies between 15 and 25 in a relatively large part of the NE area, it is inferior to 5 in nearly all of the area for the end of the century (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Yet, this decrease should be interpreted with caution, as the standard deviation of the end-of-century impact score, due to the variability across socio-economic pathways and climate models, is around 2 in the NE area that currently shows a high risk (Figure \u003cspan refid=\"MOESM3\" class=\"InternalRef\"\u003eS3\u003c/span\u003e). Comparing the current impact risk map to a risk map based on actual records for Acacia species (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e), the general pattern looks similar to potential impacts under current climate, with a less continuous distribution along the east coast, and a marked region of high risk in the Gauteng province. The latter is probably due to the fact that in the map showing impacts based on recorded occurrences, urban areas were not excluded, and much of the Gauteng province would fall under this category (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Furthermore, a positive sampling bias in this highly populated area might be responsible for higher detection rates of the species present and explain the relatively high impact risk when based on actual occurrences.\u003c/p\u003e \u003cp\u003e[insert Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e here]\u003c/p\u003e \u003cp\u003e[insert Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e here]\u003c/p\u003e \u003cp\u003eAlthough species richness of Acacias is not projected to be as high in the southwestern part of the country as along the south coast, the projected risks are similar in both areas (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). That shows that the areas around Cape Town and Hermanus, although not suitable for as many species, might remain suitable for the high impacting species. The impact score in this area remains most likely high compared to other areas even when considering its standard deviation (Figure \u003cspan refid=\"MOESM3\" class=\"InternalRef\"\u003eS3\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eWe found that potential impacts would decrease under climate change within protected areas of the Western Cape. Our results showed a lower impact risk with climate change for protected areas along the coast line, almost no change for northern protected areas and a higher impact risk in few isolated protected areas (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). However, the impact risk differed according to the protection status. The impact risk remains high in areas of high protection status (such as national parks) and low for private protected areas (with lower protection status) (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e).\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eThis study presents, to the best of our knowledge, a first attempt at combining impacts of alien taxa with their projected distribution under climate change to produce risk maps for future climatic scenarios. Our results indicate that the projected suitable areas for alien Acacia species in South Africa are reduced under climate change (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). This is similar to a study looking at \u003cem\u003eA. mangium\u003c/em\u003e and \u003cem\u003eA. auriculiformis\u003c/em\u003e in Brazil, where they found a shift in potential distribution and a reduction in suitable area (Heringer et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). More generally, this is in line with Bellard et al. (\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) who showed in a review of 71 papers covering 423 alien species that climate change is more frequently projected to contribute to a decrease in range size than an increase. However, to assume that the problem of Acacia invasions is going to sort itself out over time might be a bit too optimistic. For once, the sites at highest risk are also the regions where the highly vulnerable fynbos biome is located. This biome is already under pressure from the effects of climate change, including increased temperature and drought, and therefore more vulnerable to other pressures (Slingsby et al. \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Skowno et al. \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2021\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFurthermore, Acacia species which are not currently alien in the country (not yet introduced) were not modelled in this study. That means that other species not in our dataset might behave differently from what we present here. Therefore, this does not exclude the possibility of other Acacia species, and other alien species in general, becoming more problematic in future (see also Sheppard et al. \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Furthermore, Acacias and trees in general are long lived species which can result in a long lag between reduced climate suitability and decline in populations or observed reduction in impacts (Kowarik \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e1995\u003c/span\u003e; Robeck et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). Besides, perennial and predominantly allogame plants like Acacia species often exhibit lag phases of several decades before the beginning of an invasion (Robeck et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). Hence, several Acacia species could still be in a lag phase preceding an invasion in South Africa, and this invasion debt was not accounted for in the trends presented here (e.g., Rouget et al. \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2016\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eAs with any model, there are uncertainties linked to the projections. SDMs are based on several strong assumptions (Guisan et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Hui \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). First, occurrence records should reflect the true performance of the species. However, sampling biases are present in most if not all record databases (Beck et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Moreover, low sample size can highly impact the performance of SDMs and data available may not be sufficient to fully inform the models (Wisz et al. \u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Stockwell et al. 2002). To address these shortcomings, we designed a conservative methodology to keep only the most trustworthy records and limit oversampling in some areas. We also kept records from both the native and the alien range (Broennimann and Guisan \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2008\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eSecond, the species\u0026rsquo; performance should respond directly to the variation of the selected predictors. In this study we chose four climate variables for their known link with plant species survival and development (Mod et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2016\u003c/span\u003e) and their availability in both fine spatial and temporal scales. Moreover, Sheppard et al. (\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), showed that predictions of similar SDMs of the response of three invasive plant species to climate change was highly correlated with field experiments. However, SDMs are correlative and may yield incorrect estimates of habitat suitability if climatic variables are correlated to other unknown variables in the training area (Guisan et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Jarnevich et al. \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Future projections from our fitted models of climatic suitability could be biased by potential confounding factors of climate, as we do not take into account other factors which determine the success of alien species in a new region, such as dispersal capabilities and biotic interactions. Our model fitting procedure implicitly assumed that a species had the opportunity to colonise a large part of any country where it was recorded. Regarding interactions, Australian \u003cem\u003eAcacia\u003c/em\u003e species often have competitive advantages over native plant species and tend to become dominant among plant communities, especially after disturbance (Morris et al. \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2011\u003c/span\u003e), suggesting that the spatial extent of the realised niche would not be strongly restricted by competitive interactions compared to the potential niche. Yet, other interactions, such as mutualistic and trophic interactions with soil fungi (Birnbaum et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) may contribute to constrain the actual range of Acacias and act as confounding factors of climate in fitted SDMs, inducing bias in future projections. Furthermore, factors such as topography could be taken into account to improve the models (Bradley and Mustard \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2006\u003c/span\u003e). Thus, our model predictions could be improved by considering non-climatic drivers (such as soil composition) if these variables were available at a fine resolution.\u003c/p\u003e \u003cp\u003eThird, the species\u0026rsquo; distribution, represented by recorded occurrences, should be stable and fill any available niche in the study environment. Several studies have shown violations of the niche conservatism hypothesis during invasion with niche shifts between the native niche and the introduced niche (e.g. Parravicini et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Broennimann et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Guisan et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Moreover, predictions based on extrapolations on a new territory and with future climatic conditions may not be robust because the data used for model parameterization cannot represent all conditions in the extrapolated region (Elith and Leathwick \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2009\u003c/span\u003e; Barbet-Massin et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Sinclair et al. \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Thus, it must be kept in mind that our suitability and risk maps may be underestimated and should not be taken as a prediction of true future species richness and impacts.\u003c/p\u003e \u003cp\u003eFourth, some future climatic conditions may not have an analog amongst the historical climates of the study area (Williams \u0026amp; Jackson \u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). Yet, SDMs cannot predict exactly how species will respond to conditions that were not used for calibration (Pouteau et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). This could partly explain why most alien species (Bellard et al. \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) including Acacias alien to South Africa, are projected to experience a decrease in the size of their potential range according to our current knowledge. Future work should consider the identification of novel climates so as to avoid putting too much confidence in climates with no current analogues.\u003c/p\u003e \u003cp\u003eCaution is also advised when interpreting the risk maps including the sum of potential impacts of Acacias. Firstly, there are many ways to aggregate impacts, both, within species (calculating one impact value taking into account all impact records for the species) and across species (calculating an impact score for a site where several alien species are present) (Boulesnane-Genguant et al. in prep). Some of the most prominent methods to get one impact value per species have been to sum scores (e.g., Nentwig et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), to calculate a mean value (e.g., Rumlerova et al. 2016), and to take a maximum value (e.g., Blackburn et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) (see also Kumschick et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). They each come with underlying assumptions which can affect the results. The maximum score per species, as used here, was chosen for range shifting species in the US to anticipate if any high impacting invaders are likely to arrive under climate change (Rockwell-Postel et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). However, the applications of scores aggregated across species are scarce (for an example, see Nentwig et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2010\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFurthermore, climate change can not only affect the potential distribution of species, but also modify their invasion behaviour and impacts (e.g., Le Maitre et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Changes to fire regimes could have large impacts on alien species and native ecosystems, especially in the fynbos biome of South Africa where fires have been increasing in frequency and intensity due to alien invasions and climate change (Le Maitre et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Slingsby et al. \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Furthermore, CO\u003csub\u003e2\u003c/sub\u003e concentration could lead to woody plant densification, which is already shown for native woody plants in some southern African habitats (Skowno et al. \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). In the fynbos biome, increased CO\u003csub\u003e2\u003c/sub\u003e could also favour the alien wattles as they are nitrogen fixers as opposed to the native flora which are adapted to low nitrogen conditions (Richardson et al. 2014). In other examples, the synchronisation of the flowering period of native and alien plants may favour the latter, through increased interactions with pollinators to the detriment of the pollination of native species. This is the case in New Zealand with \u003cem\u003eCalluna vulgaris\u003c/em\u003e whose greater phenological plasticity compared to the native species \u003cem\u003eDracophyllum subulatum\u003c/em\u003e means that it can reproduce more easily in areas with a high floral density (Giejsztowt et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Morphological responses to climate change can also increase the competitiveness of alien species. For example, milder winter temperatures in China facilitate the survival of the water hyacinth \u003cem\u003ePontederia crassipes\u003c/em\u003e, and also allow it to develop a greater biomass, forming denser foliage that excludes submerged native plants the following season (You et al. \u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). In terrestrial environments, rising temperatures can also encourage the development of alien plants, reducing the availability of water for native plants, which are then at a disadvantage when it comes to coping with dry spells. This is the case of \u003cem\u003eTamarix\u003c/em\u003e spp. introduced in the United States, which develops greater capacity to capture and use water resources than native riparian species under the effect of drought (Hellmann et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2008\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eGiven all these potential interactions between climate change and impacts of alien species, what we are showing here is not an accurate representation of sites with future impacts of alien Acacia species in South Africa. However, risk maps like these of sites where impacting species could occur under climate change projections can be valuable in helping us prioritise sites for future protection from invasion impacts. Despite the potential shortcomings of the models presented here, our study contributes to our understanding of the impacts of climate change on the risk of alien Acacia invasions, including their impacts. Furthermore, it can aid the prioritisation of clearing actions for alien Acacia species in South Africa and feed into strategies for protected area management and expansion (e.g., Department of Environmental Affairs \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Protected areas are set up to safeguard biodiversity and ecosystem services for the future. We show that for the Western Cape province, most protected areas should experience less impact risk in future. However, the areas with increased risk are of lowest protection, which makes them more vulnerable to negative changes and biodiversity loss. Elevating the protection of these sites and controlling harmful Acacias can help us reduce the loss of high biodiversity areas, as stipulated in the GBF Target 1. Understanding where some of the most damaging invaders might occur under climate change is important so our management is tackling not only current, but also potential future problems.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e \u003ch2\u003eCompeting Interests\u003c/h2\u003e \u003cp\u003eThe authors have no relevant financial or non-financial interests to disclose.\u003c/p\u003e \u003c/p\u003e\u003ch2\u003eFunding\u003c/h2\u003e \u003cp\u003eSK acknowledges the support of the Centre for Invasion Biology (CIB) at Stellenbosch University and the South African Department of Forestry, Fisheries and the Environment (DFFE). SK also acknowledges funding from the B3 project. B3 (Biodiversity Building Blocks for Policy) receives funding from the European Union\u0026rsquo;s Horizon Europe Research and Innovation Programme (through grant no .584 101059592). LJ was supported by the OUSNEXT project funded by the Structure f\u0026eacute;d\u0026eacute;rative de recherche BioST of the Universit\u0026eacute; de La R\u0026eacute;union (PI: RP). OBG was supported by the Centre de Coop\u0026eacute;ration Internationale en Recherche Agronomique pour le D\u0026eacute;veloppement (CIRAD), Fonds Vert and by the Universit\u0026eacute; de La R\u0026eacute;union. This work was co-funded by the European Union: Agricultural Fund for Rural Development (EAFRD), by the Conseil D\u0026eacute;partemental de La R\u0026eacute;union, and by the Centre de Coop\u0026eacute;ration Internationale en Recherche Agronomique pour le D\u0026eacute;veloppement (CIRAD).\u003c/p\u003e\u003ch2\u003eAuthor contributions\u003c/h2\u003e \u003cp\u003eSK and MR conceptualised and designed the study. CB provided distribution data. CB, LJ, OBJ and MR performed the data analyses. The first draft of the manuscript was written by SK and all authors contributed to the writing of the submitted draft and read and approved the final manuscript.\u003c/p\u003e\u003ch2\u003eAcknowledgements\u003c/h2\u003e \u003cp\u003eWe thank Cally Jansen for data collection on impacts of alien Acacia species, and the participants of the workshop on \u0026ldquo;Managing invasions in protected areas\u0026rdquo; held in Stellenbosch from 3\u0026ndash;5 April 2024 for valuable comments received on an earlier draft of the manuscript. SK acknowledges the support of the Centre for Invasion Biology (CIB) at Stellenbosch University and the South African Department of Forestry, Fisheries and the Environment (DFFE), noting that this manuscript does not necessarily reflect the opinions of the DFFE or its employees. SK also acknowledges funding from the B3 project. B3 (Biodiversity Building Blocks for Policy) receives funding from the European Union\u0026rsquo;s Horizon Europe Research and Innovation Programme (through grant no .584 101059592). The views and opinions expressed in the present article are, however, those of the authors only and do not necessarily reflect those of the European Union or the European Commission. Neither the EU nor the EC can be held responsible for them. LJ was supported by the OUSNEXT project funded by the Structure f\u0026eacute;d\u0026eacute;rative de recherche BioST of the Universit\u0026eacute; de La R\u0026eacute;union (PI: RP). OBG was supported by the Centre de Coop\u0026eacute;ration Internationale en Recherche Agronomique pour le D\u0026eacute;veloppement (CIRAD), Fonds Vert and by the Universit\u0026eacute; de La R\u0026eacute;union. This work was co-funded by the European Union: Agricultural Fund for Rural Development (EAFRD), by the Conseil D\u0026eacute;partemental de La R\u0026eacute;union, and by the Centre de Coop\u0026eacute;ration Internationale en Recherche Agronomique pour le D\u0026eacute;veloppement (CIRAD).\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eBarbet-Massin M, Rome Q, Villemant C, Courchamp F (2018) Can species distribution models really predict the expansion of invasive species? PLoS ONE 13(3):e0193085. https://doi.org/10.1371/journal.pone.0193085\u003c/li\u003e\n\u003cli\u003eBarbet-Massin M, Thuiller W, Jiguet F (2010) How much do we overestimate future local extinction rates when restricting the range of occurrence data in climate suitability models? Ecography 33:878\u0026ndash;886\u003c/li\u003e\n\u003cli\u003eBarbet‐Massin M, Jiguet F, Albert CH, Thuiller W (2012) Selecting pseudo‐absences for species distribution models: How, where and how many?. Methods in Ecology and Evolution 3(2):327-338\u003c/li\u003e\n\u003cli\u003eBeck J, B\u0026ouml;ller M, Erhardt A, Schwanghart W (2014) Spatial bias in the GBIF database and its effect on modeling species\u0026rsquo; geographic distributions. Ecological Informatics 19:10\u0026ndash;15. https://doi.org/10.1016/j.ecoinf.2013.11.002\u003c/li\u003e\n\u003cli\u003eBellard C, Thuiller W, Leroy B, Genovesi P, Bakkenes M, Courchamp F (2013) Will climate change promote future invasions?. Global Change Biology 19(12):3740-3748\u003c/li\u003e\n\u003cli\u003eBellard C, Jeschke JM, Leroy B, Mace GM (2018) Insights from modeling studies on how climate change affects invasive alien species geography. Ecology and Evolution 8:5688\u0026ndash;5700. https://doi.org/10.1002/ece3.4098\u003c/li\u003e\n\u003cli\u003eBirnbaum C, Morald TK, Tibbett M, Bennett RG, Standish RJ (2018) Effect of plant root symbionts on performance of native woody species in competition with an invasive grass in multispecies microcosms. Ecology and Evolution 8(17):8652-8664.\u003c/li\u003e\n\u003cli\u003eBlackburn TM, Essl F, Evans T, Hulme PE, Jeschke JM, K\u0026uuml;hn I, Kumschick S, Mrugała A, Markov\u0026aacute; Z, Nentwig W, Pergl J, Py\u0026scaron;ek P, Rabitsch W, Ricciardi A, Richardson DM, Sendek A, Vil\u0026agrave; M, Wilson JRU, Winter M, Genovesi P, Bacher S (2014) A unified classification of alien species based on the magnitude of their environmental impacts. PLoS Biology 12(5): e1001850. DOI: 10.1371/journal.pbio.1001850\u003c/li\u003e\n\u003cli\u003eBoria 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\u003c/li\u003e\n\u003cli\u003eBotella C, Marchante H, Celesti-Grapow L, Brundu G, Geerts S, Ramirez-Albores JE, Gonz\u0026aacute;lez-Moreno P, Ritter M, Richardson DM (2023) The global distribution of Acacia. Wattles: Australian Acacia Species Around the World, 131. https://doi.org/10.1079/9781800622197.0009\u003c/li\u003e\n\u003cli\u003eBradley BA, Wilcove DS, Oppenheimer M (2010) Climate change increases risk of plant invasion in the Eastern United States. Biological Invasions 12:1855\u0026ndash;1872\u003c/li\u003e\n\u003cli\u003eBradley BA, Mustard JF (2006) Characterizing the landscape dynamics of an invasive plant and risk of invasion using remote sensing. Ecological Applications 16(3):1132-1147\u003c/li\u003e\n\u003cli\u003eBroennimann O, Guisan A (2008) Predicting current and future biological invasions: both native and invaded ranges matter. Biol Lett 4:585\u0026ndash;589\u003c/li\u003e\n\u003cli\u003eBroennimann O, Treier UA, M\u0026uuml;ller-Sch\u0026auml;rer H, Thuiller W, Peterson AT, Guisan A (2007) Evidence of climatic niche shift during biological invasion. Ecol Lett 10:701\u0026ndash;709\u003c/li\u003e\n\u003cli\u003eBusby JR (1991) BIOCLIM-a bioclimate analysis and prediction system. Plant Protection Quarterly 6(1):8\u0026ndash;9\u003c/li\u003e\n\u003cli\u003eChamberlain S, Barve V, Mcglinn D, Oldoni D, Desmet P, Geffert L, Ram K (2022) rgbif: Interface to the Global Biodiversity Information Facility API. R package version 3.7.2. https://CRAN.R-project.org/package=rgbif\u003c/li\u003e\n\u003cli\u003eDepartment of Environmental Affairs (2016) National protected areas expansion strategy for South Africa 2016. Department of Environmental Affairs, Pretoria, South Africa. \u003c/li\u003e\n\u003cli\u003eDoak D, Morris W (2010) Demographic compensation and tipping points in climate-induced range shifts. Nature 467:959\u0026ndash;962. doi:10.1038/nature09439\u003c/li\u003e\n\u003cli\u003eElith J, Leathwick J (2009) Species distribution models: Ecological explanation and prediction across space and time. Evolution and Systematics 40:677\u0026ndash;697\u003c/li\u003e\n\u003cli\u003eGiejsztowt J, Classen AT, Deslippe JR (2020) Climate change and invasion may synergistically affect native plant reproduction. Ecology 101(1):e02913\u003c/li\u003e\n\u003cli\u003eGuisan A, Petitpierre B, Broennimann, et al (2014) Unifying niche shift studies: insights from biological invasions. Trends in Ecology and Evolution 29:260\u0026ndash;269. https://doi.org/10.1016/j.tree.2014.02.009\u003c/li\u003e\n\u003cli\u003eGuisan A, Thuiller W, Zimmermann NE (2017) Habitat suitability and distribution models: with applications in R. Cambridge University Press, Cambridge\u003c/li\u003e\n\u003cli\u003eGuisan A, Zimmermann NE (2000) Predictive habitat distribution models in ecology. Ecological Modelling 135(2):147\u0026ndash;186. https://doi.org/10.1016/S0304-3800(00)00354-9\u003c/li\u003e\n\u003cli\u003eHawkins CL, Bacher S, Essl F, Hulme PE, Jeschke JM, K\u0026uuml;hn I, Kumschick S, Nentwig W, Pergl J, Py\u0026scaron;ek P, Rabitsch W, Richardson DM, Vil\u0026agrave; M, Wilson JRU, Genovesi P, Blackburn TM (2015) Framework and guidelines for implementing the proposed IUCN Environmental Impact Classification for Alien Taxa (EICAT). Diversity and Distributions 21(11):1360-1363\u003c/li\u003e\n\u003cli\u003eHellmann JJ, Byers JE, Bierwagen BG, Dukes JS (2008) Five potential consequences of climate change for invasive species. Conservation Biology 22(3):534-543\u003c/li\u003e\n\u003cli\u003eHeringer G, Bueno ML, Meira-Neto JA, Matos FA, Neri AV (2019) Can \u003cem\u003eAcacia mangium\u003c/em\u003e and \u003cem\u003eAcacia auriculiformis\u003c/em\u003e hinder restoration efforts in the Brazilian Atlantic Forest under current and future climate conditions?. Biological Invasions 21(9):2949-2962\u003c/li\u003e\n\u003cli\u003eHui C (2023) The dos and don\u0026rsquo;ts for predicting invasion dynamics with species distribution models. Biological Invasions 25:947\u0026ndash;953. https://doi.org/10.1007/s10530-022-02976-3\u003c/li\u003e\n\u003cli\u003eIPBES (2023) Summary for policymakers of the thematic assessment report on invasive alien species and their control of the Intergovernmental science-policy Platform on Biodiversity and Ecosystem Services. Roy HE, Pauchard A, Stoett P, Renard Truong T, Bacher S, Galil BS, Hulme PE, Ikeda T, Sankaran KV, McGeoch MA, Meyerson LA, Nu\u0026ntilde;ez MA, Ordonez A, Rahlao SJ, Schwindt E, Seebens H, Sheppard AW, Vandvik V (eds.). IPBES secretariat, Bonn, Germany. https://doi.org/10.5281/zenodo.7430692\u003c/li\u003e\n\u003cli\u003eIUCN (2020a) IUCN EICAT Categories and Criteria. The Environmental Impact Classification for Alien Taxa, 1st edn. IUCN, Gland, Switzerland and Cambridge, UK. https://doi.org/10.2305/IUCN.CH.2020.05.en\u003c/li\u003e\n\u003cli\u003eIUCN (2020b). Guidelines for using the IUCN Environmental Impact Classification for Alien Taxa (EICAT) Categories and Criteria: First edition. Version 1.1. Gland, Switzerland and Cambridge, UK. IUCN.\u003c/li\u003e\n\u003cli\u003eJansen C, Kumschick S (2022) A global impact assessment of Acacia species introduced to South Africa. Biological Invasions 24(1):175-187. https://doi.org/10.1007/s10530-021-02642-0\u003c/li\u003e\n\u003cli\u003eJarnevich CS, Stohlgren TJ, Kumar S, Morisette JT, Holcombe TR (2015) Caveats for correlative species distribution modeling. Ecol Inform 29:6\u0026ndash;15\u003c/li\u003e\n\u003cli\u003eJim\u0026eacute;nez-Valverde A, Peterson AT, Sober\u0026oacute;n J, Overton JM, Arag\u0026oacute;n P, Lobo JM (2011) Use of niche models in invasive species risk assessments. Biol Invasions 13:2785\u0026ndash;2797\u003c/li\u003e\n\u003cli\u003eKaplan H, van Niekerk A, Le Roux JJ, Richardson DM, Wilson JR (2014) Incorporating risk mapping at multiple spatial scales into eradication management plans. Biological Invasions 16:691-703\u003c/li\u003e\n\u003cli\u003eKarger DN, Conrad O, B\u0026ouml;hner J, Kawohl T, Kreft H, Soria-Auza RW, Zimmermann NE, Linder HP, Kessler M (2017) Climatologies at high resolution for the earth\u0026rsquo;s land surface areas. Scientific Data 4(1):1-20\u003c/li\u003e\n\u003cli\u003eKowarik I (1995) Time lags in biological invasions with regard to the success and failure of alien species. Plant invasions: general aspects and special problems, 15-38.\u003c/li\u003e\n\u003cli\u003eKumschick S, Jansen C (2023) Evidence-based impact assessment for naturalized and invasive Australian Acacia species. Wattles: Australian Acacia species around the world (eds. D.M. Richardson, J.J. Le Roux \u0026amp; E.M. Marchante). CABI, Wallingford, pp. 359-381.\u003c/li\u003e\n\u003cli\u003eKumschick S, Bertolino S, Blackburn TM, Brundu G, Costello KE, de Groot M, Evans T, Gallardo B, Genovesi P, Govender T, Jeschke JM, Lapin K, Measey GJ, Novoa A, Nunse AL, Probert AF, Py\u0026scaron;ek P, Preda C, Rabitsch W, Roy HE, Smith KG, Tricarico E, Vil\u0026agrave; M, Vimercati G, Bacher S (2024) Using the IUCN Environmental Impact Classification for Alien Taxa (EICAT) to inform decision-making. Conservation Biology 38. DOI:10.1111/cobi.14214\u003c/li\u003e\n\u003cli\u003eLe Maitre DC, Blignaut JN, Clulow A, Dzikiti S, Everson CS, G\u0026ouml;rgens AH, Gush MB (2020) Impacts of plant invasions on terrestrial water flows in South Africa. In Biological Invasions in South Africa (pp. 431-457). Cham: Springer International Publishing.\u003c/li\u003e\n\u003cli\u003eLe Maitre DC, Gaertner M, Marchante E et al. (2011) Impacts of invasive Australian acacias: implications for management and restoration. Diversity and Distributions 17:1015\u0026ndash;1029\u003c/li\u003e\n\u003cli\u003eLenoir J et al. (2008) A significant upward shift in plant species optimum elevation during the 20th century. Science 320:1768\u0026ndash;1771\u003c/li\u003e\n\u003cli\u003eLiaw A, Wiener M (2002) Classification and regression by randomForest. R news 2(3):18-22\u003c/li\u003e\n\u003cli\u003eMagona N, Richardson DM, le Roux JJ, Kritzinger-Klopper S, Wilson JRU (2018) Even well-studied groups of alien species might be poorly inventoried: Australian Acacia species in South Africa as a case study. NeoBiota 39:1-29. https://doi.org/10.3897/neobiota.39.23135\u003c/li\u003e\n\u003cli\u003eManes S et al. (2021) Endemism increases species\u0026rsquo; climate change risk in areas of global biodiversity importance. Biological Conservation 257. https://doi.org/10.1016/j.biocon.2021.109070.\u003c/li\u003e\n\u003cli\u003eMarchioro CA, Krechemer FS (2021) Prevention is better than cure: Integrating habitat suitability and invasion threat to assess global biological invasion risk by insect pests under climate change. Pest Management Science 77(10):4510-4520\u003c/li\u003e\n\u003cli\u003eMod HK, Scherrer D, Luoto M, Guisan A (2016) What we use is not what we know: environmental predictors in plant distribution models. Journal of Vegetation Science 27(6):1308-1322\u003c/li\u003e\n\u003cli\u003eMorris TL, Esler KJ, Barger NN, Jacobs SM, Cramer MD (2011) Ecophysiological traits associated with the competitive ability of invasive Australian acacias. Diversity and Distributions 17(5):898-910\u003c/li\u003e\n\u003cli\u003eNentwig W, Bacher S, Py\u0026scaron;ek P, Vil\u0026agrave; M, Kumschick S (2016) The Generic Impact Scoring System (GISS): a standardized tool to quantify the impacts of alien species. Environmental Monitoring and Assessment 188:315. doi:10.1371/journal.pbio.1001850.\u003c/li\u003e\n\u003cli\u003eNentwig W, K\u0026uuml;hnel E, Bacher S (2010) A generic impact‐scoring system applied to alien mammals in Europe. Conservation Biology 24(1):302-311\u003c/li\u003e\n\u003cli\u003eParravicini V, Azzurro E, Kulbicki M, Belmaker J (2015) Niche shift can impair the ability to predict invasion risk in the marine realm: an illustration using Mediterranean fish invaders. Ecol Lett 18:246\u0026ndash;253\u003c/li\u003e\n\u003cli\u003ePouteau R, Biurrun I, Brunel C, Chytr\u0026yacute; M, Dawson W, Essl F, Fristoe T, Haveman R, Hobohm C, Jansen F, Kreft H, Lenoir J, Lenzner B, Meyer C, Moeslund JE, Pergl J, Py\u0026scaron;ek P, Svenning J-C, Thuiller W, Weigelt P, Wohlgemuth T, Yang Q, van Kleunen M (2021) Potential alien ranges of European plants will shrink in the future, but less so for already naturalized than for not yet naturalized species. Diversity and Distributions 27:2063\u0026ndash;2076. https://doi.org/10.1111/ddi.13378\u003c/li\u003e\n\u003cli\u003eRichardson DM, Allsopp N, D\u0026rsquo;Antonio CM et al. (1994) Plant invasions \u0026ndash; the role of mutualisms. Biol Rev 75:65\u0026ndash;93. https://doi.org/10.1111/j.1469-185X.1999.tb00041.x\u003c/li\u003e\n\u003cli\u003eRobeck P, Essl F, van Kleunen M, Py\u0026scaron;ek P, Pergl J, Weigelt P, Mesgaran MB (2024) Invading plants remain undetected in a lag phase while they explore suitable climates. Nature Ecology \u0026amp; Evolution, 1-12.\u003c/li\u003e\n\u003cli\u003eRockwell-Postel M, Laginhas BB, Bradley BA (2020) Supporting proactive management in the context of climate change: prioritizing range-shifting invasive plants based on impact. Biological Invasions 22:2371-2383\u003c/li\u003e\n\u003cli\u003eRouget M, Robertson MP, Wilson JRU, Hui C, Essl F, Renteria JL, Richardson DM (2016) Invasion debt - quantifying future biological invasions. Diversity and Distributions 22(4):445-456\u003c/li\u003e\n\u003cli\u003eRumlerov\u0026aacute; Z, Vil\u0026agrave; M, Pergl J, Nentwig W, Py\u0026scaron;ek P (2016) Scoring environmental and socioeconomic impacts of alien plants invasive in Europe. Biological Invasions 18:3697-3711\u003c/li\u003e\n\u003cli\u003eSheppard CS, Burns BR, Stanley MC (2014) Predicting plant invasions under climate change: are species distribution models validated by field trials? Glob Change Biol 20:2800\u0026ndash;2814\u003c/li\u003e\n\u003cli\u003eSheppard CS, Burns BR, Stanley MC (2016) Future-proofing weed management for the effects of climate change: is New Zealand underestimating the risk of increased plant invasions?. New Zealand Journal of Ecology 40(3):398-405\u003c/li\u003e\n\u003cli\u003eSinclair SJ, White MD, Newell GR (2010) How useful are species distribution models for managing biodiversity under future climates? Ecology and Society 15(1):8.\u003c/li\u003e\n\u003cli\u003eSkowno AL, Thompson MW, Hiestermann J et al (2017) Woodland expansion in South African grassy biomes based on satellite observations (1990-2013): general patterns and potential drivers. Glob Chang Biol 23:2358\u0026ndash;2369. https://doi.org/10.1111/gcb.13529\u003c/li\u003e\n\u003cli\u003eSkowno AL, Jewitt D, Slingsby JA (2021) Rates and patterns of habitat loss across South Africa\u0026apos;s vegetation biomes. South African Journal of Science 117(1-2):1-5\u003c/li\u003e\n\u003cli\u003eSlingsby JA, Merow C, Aiello-Lammens M, Allsopp N, Hall S, Kilroy Mollmann H, Turner R, Wilson AM, Silander Jr JA (2017) Intensifying postfire weather and biological invasion drive species loss in a Mediterranean-type biodiversity hotspot. Proceedings of the National Academy of Sciences 114(18):4697-4702\u003c/li\u003e\n\u003cli\u003eSmith K (2020) The IUCN Red List and invasive alien species: An analysis of impacts on threatened species and extinctions. International Union for the Conservation of Nature.\u003c/li\u003e\n\u003cli\u003eStockwell DRB, Peterson AT (2002) Effects of sample size on accuracy of species distribution models. Ecol Modell 148:1\u0026ndash;13.\u003c/li\u003e\n\u003cli\u003eThomas CD, Cameron A, Green RE, Bakkenes M, Beaumont LJ, Collingham YC, Erasmus BF, De Siqueira MF, Grainger A, Hannah L, Hughes L, et al. (2004) Extinction risk from climate change. Nature 427(6970):145-148\u003c/li\u003e\n\u003cli\u003eThuiller W, Lafourcade B, Engler R, Ara\u0026uacute;jo MB (2009) BIOMOD\u0026ndash;a platform for ensemble forecasting of species distributions. Ecography 32(3):369-373\u003c/li\u003e\n\u003cli\u003eTingley R, Garc\u0026iacute;a-D\u0026iacute;az P, Arantes CRR, Cassey P (2017) Integrating transport pressure data and species distribution models to estimate invasion risk for alien stowaways. Ecography \u003c/li\u003e\n\u003cli\u003eVil\u0026agrave; M, Espinar JL, Hejda M, Hulme PE, Jaro\u0026scaron;\u0026iacute;k V, Maron JL, Pergl J, Schaffner U, Sun Y, Py\u0026scaron;ek P (2011) Ecological impacts of invasive alien plants: a meta‐analysis of their effects on species, communities and ecosystems. Ecology Letters 14(7):702-708\u003c/li\u003e\n\u003cli\u003eVil\u0026agrave; M, Trillo A, Castro-D\u0026iacute;ez P, Gallardo B, Bacher S (2024)  Field studies of the ecological impacts of invasive plants in Europe. NeoBiota 90:139-159\u003c/li\u003e\n\u003cli\u003eVolery L, Blackburn TM, Bertolino S, Evans T, Genovesi P, Kumschick S, Roy HE, Smith KG, Bacher S (2020) Improving the Environmental Impact Classification for Alien Taxa (EICAT): a summary of revisions to the framework and guidelines. NeoBiota 62:547\u0026ndash;567\u003c/li\u003e\n\u003cli\u003eWenger SJ, Olden JD (2012) Assessing transferability of ecological models: an underappreciated aspect of statistical validation. Methods in Ecology and Evolution 3(2):260-267\u003c/li\u003e\n\u003cli\u003eWilliams JW, Jackson ST (2007) Novel climates, no-analog communities, and ecological surprises. Frontiers in Ecology and the Environment 5(9):475\u0026ndash;482. https://doi.org/10.1890/070037\u003c/li\u003e\n\u003cli\u003eWisz MS, Hijmans RJ, Li J, Peterson AT, Graham CH, Guisan A (2008) Predicting Species Distributions Working Group. Effects of sample size on the performance of species distribution models. Diversity and Distributions 14:763\u0026ndash;773\u003c/li\u003e\n\u003cli\u003eYou W, Yu D, Xie D, Yu L (2013) Overwintering survival and regrowth of the invasive plant \u003cem\u003eEichhornia crassipes\u003c/em\u003e are enhanced by experimental warming in winter. Aquatic Biology 19:45-53\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"biological-invasions","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"binv","sideBox":"Learn more about [Biological Invasions](https://www.springer.com/journal/10530)","snPcode":"10530","submissionUrl":"https://submission.nature.com/new-submission/10530/3","title":"Biological Invasions","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"EICAT, biological invasions, Australian Acacia, SDM, protected area, risk map","lastPublishedDoi":"10.21203/rs.3.rs-4437291/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4437291/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eRisk maps are a useful tool to prioritise sites for management and allocate resources where they are most needed as they can show us where impacts of biological invasions are most likely to happen or expected to be largest. Given the pace of global changes, we need to understand not only the risks under current conditions, but future risks taking into account aspects such as climate change. In this study, we use Australian acacias alien to South Africa as a case study to model their potential distribution under climate change scenarios, and map their potential impacts and the uncertainty related to variability across socio-economic pathways and climatic models, to help us set better priorities for the future. We found that although climatic suitability and therefore the potential high risk area is predicted to decrease under climate change, the highly vulnerable fynbos biome remains an area with high projected impacts. We found that potential impacts would decrease under climate change within protected areas of the Western Cape. However, in this province, national parks and other areas with high levels of protection would experience greater impacts than private protected areas of lower protection status. This can help us prioritise management actions and aid the development of suitable plans to protect biodiversity not only under current, but also future climate conditions.\u003c/p\u003e","manuscriptTitle":"Mapping impacts of alien species on biodiversity in the face of climate change","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-06-11 20:12:18","doi":"10.21203/rs.3.rs-4437291/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"reviewerAgreed","content":"","date":"2024-06-13T06:17:26+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-05-26T14:34:06+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-05-21T12:32:06+00:00","index":"","fulltext":""},{"type":"submitted","content":"Biological Invasions","date":"2024-05-17T10:17:11+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"biological-invasions","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"binv","sideBox":"Learn more about [Biological Invasions](https://www.springer.com/journal/10530)","snPcode":"10530","submissionUrl":"https://submission.nature.com/new-submission/10530/3","title":"Biological Invasions","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"7f8509d2-7236-4a99-a2c5-b95d1a16af1c","owner":[],"postedDate":"June 11th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-12-23T16:03:00+00:00","versionOfRecord":{"articleIdentity":"rs-4437291","link":"https://doi.org/10.1007/s10530-024-03490-4","journal":{"identity":"biological-invasions","isVorOnly":false,"title":"Biological Invasions"},"publishedOn":"2024-12-18 15:57:42","publishedOnDateReadable":"December 18th, 2024"},"versionCreatedAt":"2024-06-11 20:12:18","video":"","vorDoi":"10.1007/s10530-024-03490-4","vorDoiUrl":"https://doi.org/10.1007/s10530-024-03490-4","workflowStages":[]},"version":"v1","identity":"rs-4437291","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4437291","identity":"rs-4437291","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}