Research Article |
Corresponding author: Emily J. McCulloch-Jones ( emilyjoy.conservation@gmail.com ) Academic editor: Joana Vicente
© 2023 Emily J. McCulloch-Jones, Tineke Kraaij, Neil Crouch, Katelyn T. Faulkner.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
McCulloch-Jones EJ, Kraaij T, Crouch N, Faulkner KT (2023) Assessing the invasion risk of traded alien ferns using species distribution models. NeoBiota 87: 161-189. https://doi.org/10.3897/neobiota.87.101104
|
Risk analysis plays a crucial role in regulating and managing alien and invasive species but can be time-consuming and costly. Alternatively, combining invasion and impact history with species distribution models offers a cost-effective and time-efficient approach to assess invasion risk and identify species for which a comprehensive risk analysis should take precedence. We conducted such an assessment for six traded alien fern species, determining their invasion risk in countries where they are traded. Four of the species (Dicksonia antarctica, Dryopteris erythrosora, Lygodium japonicum, and Phlebodium aureum) showed limited global distributions, while Adiantum raddianum and Sphaeropteris cooperi had broader distributions. A. raddianum, however, was the only species found to pose a high invasion risk in two known trade countries – the USA and Australia – and requires a complete risk analysis to determine the appropriate regulatory responses. Dicksonia antarctica, Phlebodium aureum (for New Zealand), and Dryopteris erythrosora (for the USA) posed a medium risk of invasion due to the lack of evidence of impacts, and a complete risk analysis is thus deemed less crucial for these species in these countries. For other species, suitable environments were not predicted in the countries where they are traded, thus the risk of invasion is low, and a complete risk analysis is not required. For species in countries where suitable environments are predicted but no trade information or presence data are available, risk assessments are recommended to better determine the risk posed. Despite the relatively limited potential global distribution of the studied ferns relative to other major plant invaders (e.g., Pinus spp. and Acacia spp.), their history of invasion, documented impacts in pristine environments, and high propagule pressure from trade warrants concern, possibly necessitating legislative and regulatory measures in environmentally suitable regions.
early detection, environmental suitability, horticultural trade, invasion risk, MaxEnt, risk analysis
The intentional or unintentional dispersal of species into areas outside of their native range is facilitated by various pathways of introduction which may be natural or human-mediated (
The most cost-effective means of curbing invasions is through preventative action prior to introduction (
Risk assessments comprise the initial steps of risk analysis and generally consider the likelihood of invasion alongside consequence (negative environmental or socio-economic impacts) (
The consideration of the climatic or environmental suitability of the receiving region for the taxon in question can greatly enhance the predictive capacity of risk assessments (Beaumont et al. 2014;
Species distribution modelling (SDM) is increasingly used to predict the potential distributions of alien species and identify sites that are climatically or environmentally suitable for them. These models have been extensively applied across various taxa from marine life to insects and terrestrial plants, and at various scales from local to global (
Alien ferns are generally understudied and thus are poorly represented in official plant species inventories and in invasive alien plant regulatory lists. Therefore, to detect alien fern species that have been introduced through trade,
We selected six of these alien fern species and used SDMs to determine their potential global distribution. We subsequently considered the results of these models alongside information on the species’ i) invasion status in the countries in which they are traded, ii) invasion history elsewhere, and iii) environmental or socio-economic impacts in their invaded range, to categorise each species, per trading country, in terms of the level of risk posed, and suggest the necessary response in terms of the need for risk analysis. We also indicate additional countries across the globe where the species are not yet known to occur and where risk assessment is necessary.
The six study species all have a history of invasion in numerous countries, are traded in several of the study countries (Canada, the United States of America, the United Kingdom and the Republic of Ireland, South Africa, Australia, and New Zealand), and have been introduced but are not yet naturalised or invasive in the countries where they are traded (
Maximum entropy modelling (MaxEnt) was used for the SDMs in this study as various analyses have proven MaxEnt to be a reliable predictive approach that often outperforms other methods in terms of the accuracy of the predictions, particularly for those related to biological invasions (
When developing SDMs for alien species it is recommended to include occurrence records from both the native and introduced ranges (
The quality of the occurrence records for each species was assessed using the packages Biogeo (
Although ferns are considered habitat specialists, as a group, they are known to have similar broad environmental preferences and generally select for wet habitats with moderate temperatures (i.e., avoiding temperature extremes) (
Candidate predictor variables selected for modelling the potential distribution of the six alien fern species considered in this study.
Variable | Description | Data type | Motivation for selection with regards to ferns |
---|---|---|---|
Landscape variable | |||
Land Cover | Land cover map including 30 classes describing habitat and percentage canopy cover, e.g., tree cover, broadleaved, deciduous, closed (> 40 %); or tree cover, needle-leaved, evergreen, open (15–40 %) (full details, Suppl. material |
Categorical | Accounts for the habitat and light requirements of ferns. Although many species of fern can withstand full sun and occur in bare areas, ferns are most commonly associated with shaded habitats, often in woodlands and forests ( |
WorldClim bioclimatic variables | |||
Bio 1 | Annual Mean Temperature (°C) | Continuous | An important predictor of fern richness and phylogenetic diversity at continental scales ( |
Bio 4 | Temperature Seasonality | Continuous | As for annual mean temperature. SOA. |
Bio 10 | Mean Temperature of the Warmest Quarter (°C) | Continuous | Representative of climatic extremes. Furthermore, although the global distribution of ferns is generally associated with warmer areas in the tropics ( |
Bio 11 | Mean Temperature of the Coldest Quarter (°C) | Continuous | Representative of climatic extremes. A limited number of fern species are adapted to survive sub-zero temperatures, and the majority of species do not show frond freezing tolerance ( |
Bio 12 | Annual precipitation (mm) | Continuous | The variable has been identified as one of the most important predictors of fern richness and phylogenetic diversity at continental scales ( |
Bio 15 | Precipitation seasonality (mm) | Continuous | As for annual precipitation. Additionally, fern occurrence and diversity are strongly associated with areas with many days of rain per year ( |
Bio 16 | Precipitation of the Wettest Quarter (mm) | Continuous | As for precipitation seasonality. Representative of climatic extremes. SOA. |
Bio 17 | Precipitation of the Driest Quarter (mm) | Continuous | Representative of climatic extremes. The different life stages of ferns show varying levels of desiccation tolerance ( |
The six invasive alien fern species selected for the study, their native and invaded ranges, descriptions of the climates and habitats in which they occur, trading countries for which species distribution modelling is required according to
Species | Native range | Invaded range | Climate | Habitat | Trading countries in which species distribution modelling is required | Impact in invaded range |
---|---|---|---|---|---|---|
Adiantum raddianum C.Presl | Mexico to South America | Invaded Hawaii ( |
Tropical and temperate | Herbaceous, terrestrial, lithophytic | CA, USA, AU | Displaces native species |
Dicksonia antarctica Labill. | South-eastern Australia | Naturalised in the United Kingdom ( |
Temperate | Tree fern, terrestrial | CA, USA, NZ | None reported |
Dryopteris erythrosora (D.C.Eaton) Kuntze | Eastern Asia | Naturalised in France and Belgium ( |
Temperate | Herbaceous, terrestrial | CA, USA, UK & RI, AU, NZ | None reported |
Lygodium japonicum (Thunb.) Sw. | Asia | Invasive in south-eastern USA ( |
Tropical and sub-tropical | Herbaceous, climbing | CA | Impacting the economic benefits of pine plantations and smothers indigenous vegetation |
Phlebodium aureum (L.) J.Sm. | South-eastern USA, Caribbean, South America | Invasive in Australia, South Africa ( |
Tropical and sub-tropical | Herbaceous, epiphytic, terrestrial | NZ | None reported |
Sphaeropteris cooperi (F.Muell.) R.M.Tryon | Eastern Australia | Invasive in Hawaii ( |
Temperate, tropical, and sub-tropical | Tree fern, terrestrial | US, UK & RI | Displaces native species and changes soil and plant nutrient dynamics |
Co-linearity can be detrimental to the accuracy of SDM predictions (
It has been recommended that the sites from which background records are selected should be unsuitable for the species, but should be near to the limit of what is suitable (
All models were built using MaxEnt Version 3.4.4 (http://biodiversityinformatics.amnh.org/open_source/MaxEnt/;
All models were run using five-fold cross validation (as all records are used for training and testing in this method) and model performance was evaluated using multiple methods, i) evaluating the Area Under the Curve (AUC) statistic, ii) calculating the Continuous Boyce Index (CBI), iii) assessing the fitted response curves; iv) considering the sensibility of the model in terms of fern ecology (
Information on the invasion status of the study species per trade country and their invasion history (elsewhere) was taken from
Categorisation of alien plant species in the horticultural trade in terms of the risk posed. The categorisation is based on three primary criteria: history of invasion, evidence of impact; and environmental suitability. Each level of risk is aligned with a suggestion in terms of the requirement for risk analysis. See also Suppl. material
Level of risk posed | Criteria | Requirement for risk analysis |
---|---|---|
High | Invasive or potentially invasive (i.e., invasive somewhere in the world) species for which suitable environments are available in the focus country, and for which impacts are known in its invaded range. | Complete risk analysis needed |
Medium | Invasive or potentially invasive (i.e., invasive somewhere in the world) species for which suitable environments are available in the focus country, but for which no impacts have been recorded in its invaded range. | Risk analysis may be necessary, but resources are better focussed on high risk species |
Low | Invasive or potentially invasive (i.e., invasive somewhere in the world) species, but for which no suitable environments exist in the focus country. | No further analysis necessary |
For all the models AUC values ranged between 0.76 and 0.95, indicating moderate to high performance, and the CBI values were > 0.97 indicating that predictions were consistent with the distribution of the occurrence records (Table
Model evaluation statistics for the SDMs for six selected fern species in trade. Results for both the Area Under the Curve (AUC) and Continuous Boyce Index (CBI) are shown.
Species | AUC | CBI |
---|---|---|
Adiantum raddianum | 0.76 | 1 |
Dicksonia antarctica | 0.88 | 1 |
Dryopteris erythrosora | 0.88 | 0.97 |
Lygodium japonicum | 0.85 | 1 |
Phlebodium aureum | 0.91 | 1 |
Sphaeropteris cooperi | 0.95 | 1 |
The Jack-knife test of variable contribution showed that the most important predictor differed among the species, but land cover was the most important predictor for more than one species (D. antarctica and P. aureum) (Table
Average percent contribution of the environmental predictors used in the SDMs for each species. The most important predictor for each species is in bold.
Land Cover | Temperature seasonality | Mean temperature of the warmest quarter | Mean temperature of the coldest quarter | Precipitation seasonality | Precipitation of the wettest quarter | Precipitation of the driest quarter | |
---|---|---|---|---|---|---|---|
Adiantum raddianum | 4.3 | 8.7 | 8 | 39.7 | 1 | 7.5 | 30.8 |
Dicksonia antarctica | 67.4 | 14.5 | 12.7 | 0.2 | 0 | 0.2 | 5 |
Dryopteris erythrosora | 9.4 | 7.1 | 18.4 | 3.7 | 0.1 | 39 | 22.3 |
Lygodium japonicum | 9.4 | 15.1 | 0.6 | 1.3 | 0.2 | 17.9 | 55.5 |
Phlebodium aureum | 35 | 0.3 | 0.1 | 28.7 | 10.6 | 1.6 | 23.7 |
Sphaeropteris cooperi | 17 | 25.1 | 12.6 | 21 | 0.1 | 9.3 | 14.9 |
Average contribution | 23.8 | 11.8 | 8.7 | 15.8 | 2 | 11.7 | 25.3 |
The predicted potential global distributions for most species spanned relatively few continents, with the exception of A. raddianum and S. cooperi for which suitable environments were predicted over several continents and in a greater number of countries across the globe when compared to the other study species (Fig.
The global potential distribution of six invasive alien fern species as predicted by the species distribution models. In some instances, the identified suitable environments are not easily observed or have been superimposed with multiple species occurrence records. Insets have been used in these cases to improve visibility.
Our literature search showed that environmental or socio-economic impacts have been recorded for three of the study species (A. raddianum, L. japonicum and S. cooperi; Table
Requirement for risk analysis of species in the countries in which they are traded based on, i) their invasion history, ii) invasion status in the country (all introduced and not yet naturalised or invasive), iii) whether the species is known to have impacts in its alien range, and iv) the availability of suitable environments according to the species distribution models (see Table
Species | Country | Risk rating and requirement for risk analysis |
---|---|---|
Adiantum raddianum | CA | Low risk – no need for a risk analysis |
USA | High risk – needs a complete risk analysis | |
AUS | High risk – needs a complete risk analysis | |
Dicksonia antarctica | CA | Low risk – no need for a risk analysis |
USA | Low risk – no need for a risk analysis | |
NZ | Medium risk – risk analysis needed, but not immediately | |
Dryopteris erythrosora | CA | Low risk – no need for a risk analysis |
USA | Medium risk – risk analysis needed, but not immediately | |
UK & RI | Low risk – no need for a risk analysis | |
AUS | Low risk – no need for a risk analysis | |
NZ | Low risk – no need for a risk analysis | |
Lygodium japonicum | CA | Low risk – no need for a risk analysis |
Phlebodium aureum | NZ | Medium risk – risk analysis needed, but not immediately |
Sphaeropteris cooperi | USA | Low risk – no need for a risk analysis |
UK & RI | Low risk – no need for a risk analysis |
A broad environmental and climatic tolerance is exhibited in many popular horticultural species (
Despite the satisfactory performance of the models, a few of the response curves were truncated and slightly abnormal, indicating that the occurrence records of some of the assessed species might not characterise their full fundamental niche (i.e., the full set of conditions in which a species can survive in the absence of biotic interactions;
Variables that contributed highly to the models reflected well the documented biological and environmental limitations typical for most fern species. The large contribution of precipitation variables, for example, was unsurprising as ferns generally require moist environments (
Although many ferns can colonise disturbed and altered habitats (
The consideration of environmental suitability in conjunction with the invasion status of a species and the knowledge of their invasion and impact history has enabled us to classify the study species in terms of the invasion risk they pose in the countries in which they are traded, and thus get an indication of whether they require a complete risk analysis. The USA, Australia, and New Zealand show the greatest potential for invasions by these traded alien ferns, with suitable environments available for at least four of the six study species which now constitute medium or high risk species that require risk analysis. The USA is of particular concern, as this country already has extensive documentations of fern invasions and is a prominent trader in alien fern species (
The USA, Australia, and New Zealand are most in need of trade regulations, specifically for the species identified as medium and high risk to prevent the occurrence of, or increase in, invasions. As official regulation is contingent on risk analysis, this study provides crucial information by identifying species that require risk analysis. Efforts to detect and manage escaped populations are necessary as suitable environments are available and invasive populations may remain undetected in these countries. Furthermore, countries that possess suitable environments, but lack documented occurrences of these species should conduct risk assessments. Based on the assessment outcomes, they can determine whether it is necessary to perform risk analysis and implement trade regulations to prevent future invasions. It is important to note that all species considered in this study are highly popular in trade and introduction via this pathway is thus likely (
This study is the first to model the potential global distribution of multiple invasive alien fern species, identifying countries susceptible to invasion and informing the need for risk analyses. The models reveal relatively limited potential global distributions for these ferns compared to other major invaders, but their association with undisturbed habitats, such as forests, and their documented impacts in such habitats raises concern. The projections of environmental suitability have allowed us to complete risk assessments for the studied species to inform their requirement for risk analyses in the countries in which they are traded, thus initiating the early stages of management action. The USA, Australia, and New Zealand stand out as potential hotspots for invasion by traded alien ferns, necessitating management interventions and on-ground population detection for high or medium risk species. This risk assessment approach serves as a valuable management tool, highlighting focal species for each trade country and supporting efficient resource allocation in alien species management and regulation.
We thank all data collection assistants for their time and help with the study. Support and funding for this work was provided through the DSI-NRF Centre for Excellence for Invasion Biology (C•I•B) and the South African National Biodiversity Institute (SANBI). Katelyn T. Faulkner and Emily J. McCulloch-Jones thank the South African Department of Forestry, Fisheries and the Environment (DFFE) for funding, noting that this publication does not necessarily represent the views or opinions of DFFE or its employees.
Supplementary information
Data type: docx
Explanation note: (1): Comprehensive methods on the modelling procedure applied in the study. The sub-sections are divided as in the main text to improve readability; (2): The total number of occurrence records for each study species obtained from the Global Biodiversity Information Facility (GBIF; https://www.gbif.org/) and the number of occurrence records available for modelling post-data cleaning; (3): A breakdown of the land cover variable which was acquired from the ESA CCI Land Cover project (http://maps.elie.ucl.ac.be/CCI/viewer/download.php). Left aligned numbers and land cover types are main categories, and right aligned numbers and land cover types are sub-categories; (4): Decision tree used for categorising the priority of species for risk analysis (after