
==== Front
Ecol Evol
Ecol Evol
10.1002/(ISSN)2045-7758
ECE3
Ecology and Evolution
2045-7758
John Wiley and Sons Inc. Hoboken

10.1002/ece3.70252
ECE370252
ECE-2023-10-01896.R2
Biogeography
Research Article
Research Article
Globally suitable areas for Lycorma delicatula based on an optimized Maxent model
Zhao et al.
Zhao Zhengxue https://orcid.org/0000-0003-1764-8690
1 2 3 4
Yang Lin 1 2 3
Chen Xiangsheng 1 2 3 chenxs3218@163.com

1 Institute of Entomology Guizhou University Guiyang China
2 Provincial Special Key Laboratory for Development and Utilization of Insect Resources of Guizhou Guizhou University Guiyang China
3 Guizhou Key Laboratory for Agricultural Pest Management of Mountainous Region Guizhou University Guiyang China
4 College of Agriculture Anshun University Anshun China
* Correspondence
Xiangsheng Chen, Institute of Entomology, Guizhou University, Guiyang, China.
Email: chenxs3218@163.com

20 9 2024
9 2024
14 9 10.1002/ece3.v14.9 e7025207 8 2024
31 10 2023
16 8 2024
© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

Lycorma delicatula, a globally invasive pest, has caused considerable economic losses in many countries. Determining the potential distribution range of L. delicatula is crucial for its effective management and control; however, our understanding of this species remains limited. In this study, Maxent model with occurrence records and environmental variables were fit first and then optimized by selecting the best combination of feature classes and regularization multipliers using the lowest score of corrected Akaike information criterion. Subsequently, we predicted global suitable areas for L. delicatula both currently and in the future (2041–2060, 2061–2080, and 2081–2100). The results indicated that the mean temperature of the driest quarter is the most important environmental variable limiting L. delicatula distribution. Currently, the suitable areas are concentrated in East Asia (mainly in China, South Korea, and Japan), central and eastern United States, and southern Europe. Compared with current environmental conditions, in all future climate scenarios, the number of suitable areas for L. delicatula increased. In addition, we revealed that suitable areas are likely to expand northward in the future. Our study results suggest that policymakers and governments should prioritize the development of pest management measures in suitable areas for L. delicatula, especially in high suitable areas, to control this invasive pest and minimize global economic losses.

In future climate scenarios, the suitable areas for L. delicatula were larger than the current ones. Suitable areas are likely to expand northward in the future.

climate change
invasive pest
Lycorma delicatula
Maxent model
pest management
suitable areas
Science and Technology Support Program of Guizhou Province20201Y129 National Natural Science Foundation of China 10.13039/501100001809 31860209 source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:20.09.2024
Zhao, Z. , Yang, L. , & Chen, X. (2024). Globally suitable areas for Lycorma delicatula based on an optimized Maxent model. Ecology and Evolution, 14 , e70252. 10.1002/ece3.70252
==== Body
pmc1 INTRODUCTION

Over the past two centuries, thousands of insect species have been transported worldwide (mostly unintentionally) and have colonized outside their native ranges (Liebhold et al., 2018). Invasive insects have significantly negatively affected agricultural productivity (Ziska et al., 2011), forest resources (Boyd et al., 2013), and human health (Juliano & Philip Lounibos, 2005), incurring considerable costs. It has been estimated that invasive insects cause global annual losses of at least $70 billion, with associated health costs exceeding $6.9 billion (Bradshaw et al., 2016). Various strategies have been implemented to control these invasive insects, showing varying degrees of effectiveness (McLaughlin & Dearden, 2019). The most effective strategy is to prevent their introduction (Gallien et al., 2012). Species distribution models are helpful to achieve this goal by determining their potential distribution under current and future climate conditions.

Lycorma delicatula (Hemiptera: Fulgoridae) is a planthopper species initially described in Shanxi, Shandong, and Hebei provinces in China (Liu, 1939). This pest has over 70 host plant species (Jung et al., 2022), such as black locusts (Robinia pseudoacacia), cottonwoods, willows (Salix spp.), grapes (Vitis spp.), and apples (Malus spp.) (Xiao, 1992; Zhang, 1993). This broad host range allows L. delicatula to infest various habitats, including agricultural, urban, suburban, and natural forested areas (Urban & Leach, 2023). Currently, L. delicatula has rapidly invaded many countries, including the United States (such as Maryland, New Jersey, Delaware, Virginia, and Ohio), Japan, and South Korea (Lee et al., 2019), and has become a highly notorious invasive pest. The nymphs and adults of L. delicatula frequently aggregate and harm plants directly by feeding on plant sap (Dara et al., 2015). Recent studies have shown that feeding by L. delicatula is most damaging to grapevines (Harner et al., 2022), but is less damaging than previously presumed to hardwood tree species (Hoover et al., 2023) and tree fruit (Nixon et al., 2023). For example, the damage to agriculture in the eastern United States occurred primarily in the vineyard system (Harner et al., 2022), resulting in strongly reduced yield (up to 90%), fruit quality, and sometime caused vine decline over multiple years (Urban, 2020). Moreover, a recent study predicted that L. delicatula have a high establish probability in the grape‐producing counties of California by 2033 in the western United States (Jones et al., 2022).

Numerous studies on L. delicatula have been conducted from different perspectives, including its biology (Keena & Nielsen, 2021; Kreitman et al., 2021), ethology (Keller et al., 2020; Leach & Leach, 2020), and control measures (Leach et al., 2019). These studies aimed to provide crucial information for controlling this pest. Notably, the basic condition for pest management is determining their potential distribution, which is currently not well understood for L. delicatula.

Species distribution models have become effective tools for projecting the distribution of invasive insects (Elith, 2017; Huang et al., 2019; Liu & Shi, 2020; Ramos et al., 2018). To the best of our knowledge, three studies have used these models to predict the global distribution of L. delicatula. Jung et al. (2017) identified potential high‐risk areas in the United States, Brazil, Mexico, Congo, China, Japan, and South Korea based on the CLIMEX model, but this study relied solely on distribution data from South Korea. Wakie et al. (2020) observed highly suitable areas in Asia, Oceania, South America, North America, Africa, and Europe using the Maxent model. Nevertheless, they did not optimize two key parameters (i.e., feature classes and regularization multipliers), which affected its accuracy (Phillips & Dudík, 2008) and the determination of the distribution range. Huron et al. (2022) also established the potential risk for the 50 U.S. states and 223 countries by the Maxent model. Notably, a common limitation from three studies aforementioned is that they do not predicted the global distribution range of L. delicatula under future climate change. The lack of knowledge about future distribution is not conducive to the management and control of pests, because the distribution of pests often shifts in the future relative to the current time, which makes that existing pest management strategies may become inadequate and require modification. Thus, designing an excellent species distribution model for forecasting the future global distribution of L. delicatula is urgent.

The Maxent model is based on a machine learning algorithm that estimates the distribution (geographic range) of a target species by determining the probability distribution of maximum entropy while considering constraints derived from environmental variables at several occurrences (Phillips et al., 2017). This model is the most commonly used species distribution model owing to its superior performance over other models (Merow et al., 2013; Nair & Peterson, 2023; Phillips et al., 2006; Zhou et al., 2023). Thus, in this study, the Maxent model was employed to predict the current and future potential suitable areas of L. delicatula globally. Our study mainly addresses two key questions: (1) Where are the current and future potential suitable areas located? (2) How will the potential suitable areas change under future climate conditions?

2 MATERIALS AND METHODS

2.1 Species distribution data

The 15,599 raw occurrence records of L. delicatula were collected from the literature and Global Biodiversity Information Facility (GBIF, 2023). The GBIF occurrence records with common errors in biological collections (e.g., sea coordinates and coordinate–country mismatches) and high coordinate uncertainty (>20 km) were cleaned using the CoordinateCleaner Package via R 4.2.1 (Zizka et al., 2019). Furthermore, to reduce the influence of sampling bias on our prediction results, the occurrence records were spatially thinned at a distance of 20 km, based on the spThin package via R 4.2.1 (Aiello‐Lammens et al., 2015). This process retains the maximum number of occurrence records for a given 20 km using a randomization approach. Consequently, 464 occurrence records were obtained globally, and they are distributed in China (e.g., Guizhou and Yunan provinces), South Korea, southern Japan, and the United States (e.g., North Carolina, Tennessee, and Virginia states) (Figure 1).

FIGURE 1 Global occurrence records of L. delicatula.

2.2 Environmental variables

The 19 bioclimatic variables (Bio1–Bio19) and altitude data were widely used to predict the suitable areas of insects in the species distribution models (e.g., Wei et al., 2018, 2019; Zhao et al., 2024a). Bioclimatic variables are related to temperature and precipitation, and represent annual trends (e.g., annual precipitation), seasonality (e.g., isothermality), and limiting factors (e.g., mean temperature of the wettest quarter and driest quarter) (WorldClim website, https://www.worldclim.org). In this study, 19 bioclimatic variables recorded from 1970 to 2000 were selected. To reduce collinearity among 20 environmental variables, we followed a specific process. First, environmental variable values from the occurrence records were extracted using ArcGIS 10.7. Second, pairwise Pearson correlation coefficients (r) were calculated among the environmental variables using SPSS 25. Only one variable was retained when the correlation coefficient of a pair of variables |r| was ≥0.85. Third, variance inflation factor (VIF) value were computed for the retained variables using the same software, and the variable with the highest VIF value was omitted until it was <5. Finally, five environmental variables were selected for subsequent analysis: isothermality (Bio3), mean temperature of the wettest quarter (Bio8), mean temperature of the driest quarter (Bio9), annual precipitation (Bio12), and altitude.

Future bioclimatic variables were derived from three timeframes (2041–2060, 2061–2080, and 2081–2100) based on data from the Coupled Model Intercomparison Project Phase 6 with two shared socioeconomic pathways (SSP126 and SSP585). To control uncertainty in future climate models, we averaged climate data from the CanESM5, IPSL‐CM6A‐LR, and MIROC6 models (Zhao et al., 2024a). Moreover, we used current altitude data for future predictions, which are expected to undergo minimal changes over the next few decades.

All environmental variables were downloaded from the WorldClim website (https://www.worldclim.org) and had a spatial resolution of 5 arc min.

2.3 Optimized Maxent model

The feature classes (FC) and regularization multipliers (RM) in the Maxent model have affected its performance (Phillips & Dudík, 2008; Radosavljevic et al., 2014), indicating that these two parameters should be optimized (Muscarella et al., 2014). A FC is a simple mathematical transformation of an environmental variable and determine the kinds of constraints allowed in a model (Velasco & González‐Salazar, 2019). More FC enables more flexible and complex fits to the observed data, but they may require more data (Merow et al., 2013; Phillips & Dudík, 2008). RM is a value that seeks to balance model fit and complexity by warranting that predicted values do not fit too exactly the empirical constraints of environmental variables (Phillips & Dudík, 2008; Velasco & González‐Salazar, 2019). Here, using the ENMeval package in R 4.2.1, we employed RM values ranging from 0.5 to 4 (increments of 0.5) and eight different FC combinations (L, LQ, LQH, LQHP, LQHPT, QHP, QHPT, and HPT; L = linear, Q = quadratic, H = hinge, P = product, and T = threshold) to select the model with the lowest score of corrected Akaike information criterion (Muscarella et al., 2014). Consequently, LQHPT for FC and 1 for RM values were obtained. Five replications with cross‐validation and randomly 10,000 background points located across the globe were selected to run the model, and the cloglog output format was used, as it is the most appropriate for estimating the probability of presence (Phillips et al., 2017).

Model performance was evaluated using the area under the receiver operating characteristic curve (AUC) and the true skill statistic (TSS), both of which are considered excellent when their values exceed 0.9 and 0.8, respectively (Ben Rais Lasram et al., 2010; Bogawski et al., 2019). The relative importance of the environmental variables was assessed using the jackknife test.

2.4 Changes in suitable areas

We applied the threshold of maximum training sensitivity plus specificity (0.17) to divide suitable/unsuitable areas. Suitable areas were further divided into three levels: low suitable area (0.17–0.4), moderate suitable area (0.4–0.6), and high suitable area (>0.6).

To describe the changes in suitable areas under different climatic conditions, we identified patterns of contraction, expansion, and no change under future conditions relative to those under the current conditions. This analysis was conducted using ArcGIS 10.7 with SDMtoolbox (Brown, 2014).

3 RESULTS

3.1 Model evaluation

All AUC and TSS values for each run of the Maxent model exceeded 0.9 (Table 1). Moreover, the mean AUC and TSS values of the five replicates were 0.9716 and 0.9188, respectively, indicating the excellent performance of the model established in this study (Table 1).

TABLE 1 AUC and TSS values for each run and the mean from the Maxent model.

Replications	AUC	TSS	
1	0.9711	0.9076	
2	0.9715	0.9379	
3	0.9755	0.9219	
4	0.9723	0.9117	
5	0.968	0.9152	
Mean	0.9716	0.9188	

3.2 Importance of environmental variables

The results of the jackknife test revealed relative importance of the five environmental variables in determining L. delicatula distribution (Figure 2). The mean temperature of the driest quarter (Bio9) emerged as the most important environmental variable, followed by annual precipitation (Bio12), isothermality (Bio3), and mean temperature of the wettest quarter (Bio8) (Figure 2). Evidently, altitude showed the lowest importance (Figure 2).

FIGURE 2 Environmental variable importance using the jackknife test.

Based on the response curves (Figure 3), the suitability of the environmental variables for L. delicatula was as follows: 0.48%–37.72% for isothermality (Bio3), 1.57–47.03°C for the mean temperature of the wettest quarter (Bio8), −8.38 to 28.49°C for the mean temperature of the driest quarter (Bio9), 426.15–7417.3 mm for annual precipitation (Bio12), and −651.7 to 6616.7 m for altitude.

FIGURE 3 Response curves of the relationship between the presence probability of L. delicatula and five environmental variables.

3.3 Patterns and changes in suitable areas

Under current environmental conditions, suitable areas for L. delicatula were mainly distributed in East Asia, central and eastern United States, and southern Europe (Figure 4), totaling an area of 7.89 × 106 km2 (Table 2). Among them, the high, moderate, and low suitable regions covered areas of 1.89 × 106, 1.62 × 106, and 4.38 × 106 km2, respectively (Table 2). In East Asia, suitable areas were mainly concentrated in China, South Korea, and Japan (Figure 4). The central and eastern United States were predicted to be highly suitable areas for L. delicatula (Figure 4). Conversely, in southern Europe, low suitable areas were dominant, with high and moderate suitable areas being minimal (Figure 4).

FIGURE 4 Currently suitable areas for L. delicatula on a global scale.

TABLE 2 Areas with suitability under current and future climatic conditions (km2).

Levels	Current	2041–2060	2061–2080	2081–2100	
SSP126	SSP585	SSP126	SSP585	SSP126	SSP585	
Total	7.89 × 106	1.17 × 107	1.24 × 107	1.15 × 107	1.4 × 107	1.16 × 107	1.64 × 107	
Low	4.38 × 106	7.39 × 106	8.19 × 106	7.34 × 106	9.7 × 106	7.35 × 106	1.19 × 107	
Moderate	1.62 × 106	2.12 × 106	2.20 × 106	2.1 × 106	2.6 × 106	2.18 × 106	3.21 × 106	
High	1.89 × 106	2.10 × 106	2.09 × 106	2.13 × 106	1.73 × 106	2.09 × 106	1.38 × 106	

Under the 2041–2060 SSP126 scenario, suitable areas for L. delicatula were projected to cover an area of 1.17 × 107 km2 (Table 2), indicating an increasing trend relative to current climatic conditions (Figure 5a). The high, low, and moderate suitable areas also increased consistently (Table 2). Under the SSP126 scenario of 2041–2060, the range expansion region reached an area of 4.51 × 106 km2 (Table 3) and was concentrated in the north of suitable areas, mainly distributed in China, Europe, the United States, and Canada (Figure 6a). The region with no change in suitability covered the largest area at 7.1 × 106 km2 (Table 3). Additionally, the range contraction region covered a limited area of 4.85 × 105 km2 (Table 3). According to the 2041–2060 SSP585 scenario, suitable areas for L. delicatula covered an area of 1.24 × 107 km2 (Table 2), which was greater than that under current conditions (Figure 5b). Compared with the SSP126 scenario, the areas of regions with range expansion and contraction were greater (Figure 6b), whereas the area of region with no change decreased (Table 3).

FIGURE 5 Suitable areas for L. delicatula during 2041–2060 on a global scale. SSP126 and SSP585 represent future climate scenarios for shared socioeconomic pathways 126 and 585, respectively.

TABLE 3 Changes in suitable areas under future conditions compared with the current conditions (km2).

Types	2041–2060	2061–2080	2081–2100	
SSP126	SSP585	SSP126	SSP585	SSP126	SSP585	
Range expansion	4.51 × 106	5.54 × 106	4.52 × 106	7.62 × 106	4.5 × 106	1.06 × 107	
Range contraction	4.85 × 105	6.51 × 105	5.4 × 105	1.17 × 106	4.7 × 105	1.73 × 106	
No change	7.1 × 106	6.94 × 106	7.05 × 106	6.42 × 106	7.12 × 106	5.86 × 106	

FIGURE 6 Suitable area changes for L. delicatula during 2041–2060 compared with the current environmental conditions. SSP126 and SSP585 represent future climate scenarios for shared socioeconomic pathways 126 and 585, respectively.

Under the 2061–2080 SSP126 scenario, the distribution pattern of suitable areas and the three different suitability levels were expected to increase compared to the current conditions (Figure 7a, Table 2). Similar to the two 2041–2060 scenarios, suitable areas expanded primarily toward the north (Figure 8a). Notably, in the SSP585 scenario, suitable area was predicted to occupy a larger distribution range relative to SSP126 scenario (Figure 7b), and the expansion of suitable areas into the north increased compared to the SSP585 scenario for 2041–2060 (Figure 8b). In this timeframe, the expansion, contraction, and no change areas covered 7.62 × 106 km2, 1.17 × 106 km2, and 6.42 × 106 km2, respectively (Table 3).

FIGURE 7 Suitable areas for L. delicatula during 2061–2080 on a global scale. SSP126 and SSP585 represent future climate scenarios for shared socioeconomic pathways 126 and 585, respectively.

FIGURE 8 Suitable area changes for L. delicatula during 2061–2080 compared with the current environmental conditions. SSP126 and SSP585 represent future climate scenarios for shared socioeconomic pathways 126 and 585, respectively.

Under the 2081–2100 SSP126 scenario, suitable areas for L. delicatula were projected to increase compared with the current environmental conditions (Figure 9a) and reach 1.16 × 107 km2 (Table 2), with high suitable areas predicted to cover 2.09 × 106 km2 (Table 2). The increasing suitable areas were predicted to span 4.5 × 106 km2 (Table 3), primarily expanding to the north (Figure 10a). Areas with no change and contraction measured 7.12 × 106 km2 and 4.7 × 105 km2, respectively (Table 3). Under the SSP585 scenario, the projected suitable areas reached a maximum value (1.64 × 107 km2) in all future climate scenarios, whereas high suitable areas become the least (Figure 9b, Table 2). Notably, the suitable areas expanding northward exhibited the widest geographic range (Figure 10b), substantially surpassing those of other future climate scenarios.

FIGURE 9 Suitable areas for L. delicatula during 2081–2100 on a global scale. SSP126 and SSP585 represent future climate scenarios for shared socioeconomic pathways 126 and 585, respectively.

FIGURE 10 Suitable area changes for L. delicatula during 2081–2100 compared with the current environmental conditions. SSP126 and SSP585 represent future climate scenarios for shared socioeconomic pathways 126 and 585, respectively.

4 DISCUSSION

Based on occurrence records and environmental variables, we developed an optimized Maxent model to predict current and future global suitable areas for L. delicatula and further explore changes in these suitable areas. Two evaluation metrics (AUC and TSS) consistently demonstrated the exceptional performance of the model.

The model results showed that the mean temperature of the driest quarter (Bio9) was the most important environmental variable influencing L. delicatula distribution. The importance of this variable was demonstrated in the distribution of this pest in the United States using the Maxent model (Wakie et al., 2020). Consistent with previous studies on distribution modeling of other insects (Abou‐Shaara et al., 2021; Amaro et al., 2022; Wei et al., 2018), temperature‐related variables were more important than precipitation‐related variables. This could be attributed to insects being poikilotherms, with temperature strongly affecting their physiological activities (Zhang et al., 2023) and thus predominantly determining their geographical range. In addition, it is crucial to acknowledge that altitude variables undeniably influence insect distribution, as shown in previous studies (Schuldt & Assmann, 2009, 2011). However, in this study, the effect of altitude on L. delicatula distribution is lowest in all environmental variables.

The most important economic crops damaged by L. delicatula is grapevines. California grows more grapes than any other state in the United States and produces 82% of the national grape crop. Several previous studies have predicted that California is suitable for L. delicatula to survive in the current and future (Jung et al., 2017; Jones et al., 2022; Wakie et al., 2020), whereas our predictions results showed that California is not suitable. Thus, this suggests that concerns about L. delicatula being able to damage California grapevines may be overdone. In addition, previous research results also show that the top five grape‐growing Europe countries including Spain, France, Italy, Portugal, and Romania are also environmentally suitable for L. delicatula (Jung et al., 2017; Wakie et al., 2020). In our study, suitable areas were not found in Portugal, and the remaining four countries had smaller range of suitable areas. However, the prediction results further indicated that these suitable areas will exist in the future, suggesting that the potential negative economic impact of L. delicatula on grape production is permanent. It is worth noting that Wakie's study (Wakie et al., 2020) also used the Maxent model, but did not strictly optimize the combinations of FC and RM that seriously affect the model performance and prediction results; however, we did. In this case, suitable areas predicted by our study are more accurate.

Numerous studies have demonstrated that pest distributions can either expand or contract under the influence of climate change (Lee et al., 2021; Ramasamy et al., 2021; Wang et al., 2021; Wang et al., 2020; Xu et al., 2020; Zhao et al., 2024b). In our study, we observed that the suitable habitat for L. delicatula is expected to expand significantly (mainly northward) in the future, indicating increased invasive potential. One typical region affected by this expansion is Europe. Under current climate conditions, suitable areas appear only in southern Europe. However, under all future climate scenarios, especially the SSP585 scenarios  for 2061–2080 and 2081–2100, suitable areas are projected to have expanded into northern Europe. This expansion may be due to the region becoming warmer in the future, according to our findings that temperature is the main factor affecting L. delicatula distribution (Figure 2) and the occurrence rate in high‐temperature conditions is expected to increase (Figure 3). Additionally, we observed that the high suitable areas did not consistently increase in the future compared to the present, which aligns with findings from previous studies (e.g., Gao et al., 2022; Wei et al., 2019). Enhanced management strategies are essential for regions with high suitable habitats, such as the central and eastern United States. Furthermore, we noted that the suitable areas for L. delicatula under the SSP585 scenario consistently exceeded those under the SSP126 scenario. Since SSP585 represents a high‐forcing scenario and SPP126 represents a low‐forcing scenario, the former predicts higher temperatures in the future. Consequently, our results are consistent with the understanding that warmer temperatures increase opportunities for L. delicatula population establishment (Figure 3).

The suitable areas determined in this study are valuable for developing pest management strategies to effectively control L. delicatula. First, policymakers and governments should prioritize substantial efforts to manage L. delicatula in its current and future suitable areas, especially in highly suitable areas. Second, a range of appropriate and sustained proactiveness management strategies, such as monitoring and strict quarantine measures, should be established in infested areas. Control measures must be promptly developed, and this pest should be eliminated upon discovery.

A common assumption in species distribution models is that species distributions are in equilibrium with environments (i.e., species occur in all climatically suitable areas). For invasive species, the invasion process may not be complete; therefore, such species may not have occupied all suitable environments (Wilson et al., 2007). Consequently, the modeling for invasive species may be biased (Broennimann & Guisan, 2008). To solve this problem, a classical approach is to train the models on native distributional areas (i.e., only using native occurrence records) because these areas were considered to have the advantage of achieving a higher probability of distributional equilibrium (Jiménez‐Valverde et al., 2011). Nevertheless, although the inclusion of invasion records in species distribution models is controversial, it is widely argued that these records can better approximate future invasion potential (Morey & Venette, 2020). Most importantly, recent studies have shown that models using occurrence records from both native and invasive areas accurately predict potential distribution of invasive species (e.g., Etges et al., 2023; Magory, 2019). Thus, the global potential distribution of L. delicatula predicted by this study model is still valuable.

Species distribution is influenced by a combination of biotic and abiotic factors, and incorporating both factors into species distribution models can improve the accuracy of the model (de Araújo et al., 2014; Giannini et al., 2013; Lippitt et al., 2008). Unfortunately, in this study, we did not consider a small range of biotic factors that influence L. delicatula distribution when constructing the Maxent model. Host plants influence the geographical range of L. delicatula, as demonstrated by a recent study by Jung et al. (2022). This study reported that tree height and diameter at the root collar of Ailanthus altissima trees, one of the main host plants for L. delicatula, significantly increased the number of egg masses in this pest, and thus Ailanthus altissima trees favor population growth and increases invasion risk. Additionally, studies have shown that human activity increases the risk of L. delicatula invasion (Cook et al., 2021; Ladin et al., 2023). Therefore, for more accurate predictions in the future, host plants and human activity variables such as the human influence index should be included as predictor variables in species distribution models.

5 CONCLUSIONS

We developed an optimized Maxent model to predict the current and future suitable areas for the invasive pest L. delicatula worldwide. Our model results predicted that East Asia, central and eastern United States, and southern Europe are the main suitable areas under current environmental conditions. However, in future climate scenarios, the suitable areas for L. delicatula were larger than the current ones, indicating a significant increase in the invasive potential of this pest. Furthermore, our results indicated a northward expansion trend of suitable areas. This study provides valuable insights into the risk of introducing L. delicatula under both current and future conditions and can be used to develop appropriate pest management strategies.

AUTHOR CONTRIBUTIONS

Zhengxue Zhao: Conceptualization (lead); formal analysis (lead); visualization (lead); writing – original draft (lead). Lin Yang: Conceptualization (supporting); data curation (lead); formal analysis (supporting). Xiangsheng Chen: Methodology (lead); visualization (supporting); writing – original draft (supporting).

FUNDING INFORMATION

This research was funded by the National Natural Science Foundation of China (No. 31860209) and the Science and Technology Support Program of Guizhou Province (No. 20201Y129).

CONFLICT OF INTEREST STATEMENT

We declare no conflict of interest in connection with the work submitted.

ACKNOWLEDGEMENTS

We thank the Aiello‐Lammens Matthew E for providing help in spThin package. We also like to express our sincere thanks to editor and the two anonymous reviewers.

DATA AVAILABILITY STATEMENT

Raw distribution data, environment variables, and R code are available in figshare: https://doi.org/10.6084/m9.figshare.24460603.
==== Refs
REFERENCES

Abou‐Shaara, H. , Alashaal, S. A. , Hosni, E. M. , Nasser, M. G. , Ansari, M. J. , & Alharbi, S. A. (2021). Modeling the invasion of the large hive beetle, Oplostomus fuligineus, into North Africa and South Europe under a changing climate. Insects, 12 (4 ), 275.33804941
Aiello‐Lammens, M. E. , Boria, R. A. , Radosavljevic, A. , Vilela, B. , & Anderson, R. B. (2015). spThin: An R package for spatial thinning of species occurrence records for use in ecological niche models. Ecography, 38 (5 ), 541–545.
Amaro, G. , Fidelis, E. G. , de Medeiros, C. M. , & da Silva, R. S. (2022). Risk analysis of the spread of the quarantine pest mite Schizotetranychus hindustanicus in Brazil. Experimental and Applied Acarology, 88 (3–4 ), 263–275.36344861
Ben Rais Lasram, F. , Guilhaumon, F. , Albouy, C. , Somot, S. , Thuiller, W. , & Mouillot, D. (2010). The Mediterranean Sea as a ‘cul‐de‐sac’ for endemic fishes facing climate change. Global Change Biology, 16 (12 ), 3233–3245.
Bogawski, P. , Damen, T. , Nowak, M. M. , Pedziwiatr, K. , Wilkin, P. , Mwachala, G. , Pierzchalska, J. , & Wiland‐Szymanska, J. (2019). Current and future potential distributions of three dracaena Vand. Ex L. species under two contrasting climate change scenarios in Africa. Ecology and Evolution, 9 (12 ), 6833–6848.31380019
Boyd, I. L. , Freer‐Smith, P. H. , Gilligan, C. A. , & Godfray, H. C. J. (2013). The consequence of tree pests and diseases for ecosystem services. Science, 342 (6160 ), 1235773.24233727
Bradshaw, C. J. , Leroy, B. , Bellard, C. , Roiz, D. , Albert, C. , Fournier, A. , Barbet‐Massin, M. , Salles, J. M. , Frédéric, S. , & Courchamp, F. (2016). Massive yet grossly underestimated global costs of invasive insects. Nature Communications, 7 (1 ), 12986.
Broennimann, O. , & Guisan, A. (2008). Predicting current and future biological invasions: Both native and invaded ranges matter. Biology Letters, 4 (5 ), 585–589.18664415
Brown, J. L. (2014). SDMtoolbox: A python‐based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods in Ecology and Evolution, 5 (7 ), 694–700.
Cook, R. T. , Ward, S. F. , Liebhold, A. M. , & Fei, S. (2021). Spatial dynamics of spotted lanternfly, Lycorma delicatula, invasion of the northeastern United States. NeoBiota, 70 , 23–42.
Dara, S. K. , Barringer, L. , & Arthurs, S. P. (2015). Lycorma delicatula (Hemiptera: Fulgoridae): A new invasive pest in the United States. Journal of Integrated Pest Management, 6 (1 ), 20.
de Araújo, C. B. , Marcondes‐Machado, L. O. , Costa, G. C. , & Silman, M. (2014). The importance of biotic interactions in species distribution models: A test of the Eltonian noise hypothesis using parrots. Journal of Biogeography, 41 (3 ), 513–523.
Elith, J. (2017). Predicting distributions of invasive species. In A. P. Robinson , T. Walshe , M. A. Burgman , & M. Nunn (Eds.), Invasive species: Risk assessment and management (pp. 93–129). Cambridge University Press.
Etges, M. F. , Martínez‐Lanfranco, J. A. , & Guadagnin, D. L. (2023). Spread risk assessment of invasive axis deer using bioclimatic niche models. Biological Invasions, 25 (2 ), 485–498.
Gallien, L. , Douzet, R. , Pratte, S. , Zimmermann, N. E. , & Thuiller, W. (2012). Invasive species distribution models‐how violating the equilibrium assumption can create new insights. Global Ecology and Biogeography, 21 (11 ), 1126–1136.
Gao, X. , Zhao, Q. , Wei, J. , & Zhang, H. (2022). Study on the potential distribution of Leptinotarsa decemlineata and its natural enemy Picromerus bidens under climate change. Frontiers in Ecology and Evolution, 9 , 786436.
GBIF . (2023). Global Biodiversity Information Facility. https://www.gbif.org/
Giannini, T. C. , Chapman, D. S. , Saraiva, A. M. , Alves‐dos‐Santos, I. , & Biesmeijer, J. C. (2013). Improving species distribution models using biotic interactions: A case study of parasites, pollinators and plants. Ecography, 36 (6 ), 649–656.
Harner, A. D. , Leach, H. L. , Briggs, L. , & Centinari, M. (2022). Prolonged phloem feeding by the spotted lanternfly, an invasive planthopper, alters resource allocation and inhibits gas exchange in grapevines. Plant Direct, 6 (10 ), e452.36226305
Hoover, K. , Iavorivska, L. , Lavely, E. K. , Uyi, O. , Walsh, B. , Swackhamer, E. , Johnson, A. , & Eissenstat, D. M. (2023). Effects of long‐term feeding by spotted lanternfly (Hemiptera: Fulgoridae) on ecophysiology of common hardwood host trees. Environmental Entomology, 52 (5 ), 888–899.37643187
Huang, M. , Ge, X. , Shi, H. , Tong, Y. , & Shi, J. (2019). Prediction of current and future potential distributions of the eucalyptus pest Leptocybe invasa (Hymenoptera: Eulophidae) in China using the CLIMEX model. Pest Management Science, 75 (11 ), 2958–2968.30868710
Huron, N. A. , Behm, J. E. , & Helmus, M. R. (2022). Paninvasion severity assessment of a U.S. grape pest to disrupt the global wine market. Communications Biology, 5 (1 ), 655.35788172
Jiménez‐Valverde, A. , Peterson, A. T. , Soberón, J. O. J. M. , Overton, J. M. , Aragón, P. , & Lobo, J. M. (2011). Use of niche models in invasive species risk assessments. Biological Invasions, 13 , 2785–2797.
Jones, C. , Skrip, M. M. , Seliger, B. J. , Jones, S. , Wakie, T. , Takeuchi, Y. , Petras, V. , Petrasova, A. , & Meentemeyer, R. K. (2022). Spotted lanternfly predicted to establish in California by 2033 without preventative management. Communications Biology, 5 (1 ), 558.35676315
Juliano, S. A. , & Philip Lounibos, L. (2005). Ecology of invasive mosquitoes: Effects on resident species and on human health. Ecology Letters, 8 (5 ), 558–574.17637849
Jung, J. M. , Jung, S. , Byeon, D. H. , & Lee, W. H. (2017). Model‐based prediction of potential distribution of the invasive insect pest, spotted lanternfly L. delicatula (Hemiptera: Fulgoridae), by using CLIMEX. Journal of Asia‐Pacific Biodiversity, 10 (4 ), 532–538.
Jung, M. , Kho, J. W. , Gook, D. H. , Lee, Y. S. , & Lee, D. H. (2022). Dispersal and oviposition patterns of L. delicatula (Hemiptera: Fulgoridae) during the oviposition period in Ailanthus altissima (Simaroubaceae). Scientific Reports, 12 (1 ), 9972.35705683
Keena, M. A. , & Nielsen, A. L. (2021). Comparison of the hatch of newly laid Lycorma delicatula (Hemiptera: Fulgoridae) eggs from the United States after exposure to different temperatures and durations of low temperature. Environmental Entomology, 50 (2 ), 410–417.33399184
Keller, J. A. , Johnson, A. E. , Uyi, O. , Wurzbacher, S. , Long, D. , & Hoover, K. (2020). Dispersal of L. delicatula (Hemiptera: Fulgoridae) nymphs through contiguous, deciduous forest. Environmental Entomology, 49 (5 ), 1012–1018.32814958
Kreitman, D. , Keena, M. A. , Nielsen, A. L. , Hamilton, G. , & Brent, C. (2021). Effects of temperature on development and survival of nymphal Lycorma delicatula (Hemiptera: Fulgoridae). Environmental Entomology, 50 (1 ), 183–191.33269378
Ladin, Z. S. , Eggen, D. A. , Trammell, T. L. E. , & D'Amico, V. (2023). Human‐mediated dispersal drives the spread of the spotted lanternfly (Lycorma delicatula). Scientific Reports, 13 (1 ), 1098.36658159
Leach, H. , & Leach, A. (2020). Seasonal phenology and activity of spotted lanternfly (Lycorma delicatula) in eastern US vineyards. Journal of Pest Science, 93 (4 ), 1215–1224.
Leach, H. , Biddinger, D. J. , Krawczyk, G. , Smyers, E. , & Urban, J. M. (2019). Evaluation of insecticides for control of the spotted lanternfly, Lycorma delicatula, (Hemiptera: Fulgoridae), a new pest of fruit in the northeastern U.S. Crop Protection, 124 , 104833.
Lee, C. M. , Lee, D. S. , Kwon, T. S. , Athar, M. , & Park, Y. S. (2021). Predicting the global distribution of Solenopsis geminata (Hymenoptera: Formicidae) under climate change using the MaxEnt model. Insects, 12 (3 ), 229.33800231
Lee, D. H. , Park, Y. L. , & Leskey, T. C. (2019). A review of biology and management of L. delicatula (Hemiptera: Fulgoridae), an emerging global invasive species. Journal of Asia‐Pacific Entomology, 22 (2 ), 589–596.
Liebhold, A. M. , Yamanaka, T. , Roques, A. , Augustin, S. , Chown, S. L. , Brockerhoff, E. G. , & Pyšek, P. (2018). Plant diversity drives global patterns of insect invasions. Scientific Reports, 8 (1 ), 12095.30108295
Lippitt, C. D. , Rogan, J. , Toledano, J. , Sangermano, F. , Eastman, J. R. , Mastro, V. , & Sawyer, A. (2008). Incorporating anthropogenic variables into a species distribution model to map gypsy moth risk. Ecological Modelling, 210 (3 ), 339–350.
Liu, G. (1939). Some extracts from the history of entomology in China. Psyche, 46 , 23–28.
Liu, Y. , & Shi, J. (2020). Predicting the potential global geographical distribution of two icerya species under climate change. Forests, 11 (6 ), 684.
Magory Cohen, T. , McKinney, M. , Kark, S. , & Dor, R. (2019). Global invasion in progress: Modeling the past, current and potential global distribution of the common myna. Biological Invasions, 21 , 1295–1309.
McLaughlin, G. M. , & Dearden, P. K. (2019). Invasive insects: Management methods explored. Journal of Insect Science, 19 (5 ), 17.
Merow, C. , Smith, M. J. , & Silander, J. A. (2013). A practical guide to MaxEnt for modeling species' distributions: What it does, and why inputs and settings matter. Ecography, 36 (10 ), 1058–1069.
Morey, A. C. , & Venette, R. C. (2020). Minimizing risk and maximizing spatial transferability: Challenges in constructing a useful model of potential suitability for an invasive insect. Annals of the Entomological Society of America, 113 (2 ), 100–113.
Muscarella, R. , Galante, P. J. , Soley‐Guardia, M. , Boria, R. A. , Kass, J. M. , Uriarte, M. , Anderson, R. P. , & McPherson, J. (2014). ENMeval: An R package for conducting spatially independent evaluations and estimating optimal model complexity for Maxent ecological niche models. Methods in Ecology and Evolution, 5 (11 ), 1198–1205.
Nair, R. R. , & Peterson, A. T. (2023). Mapping the global distribution of invasive pest Drosophila suzukii and parasitoid Leptopilina japonica: Implications for biological control. PeerJ, 11 , e15222.37123003
Nixon, L. J. , Barnes, C. , Deecher, E. , Madalinska, K. , Nielsen, A. , Urban, J. , & Leskey, T. C. (2023). Evaluating deployment strategies for spotted lanternfly (Lycorma delicatula Hemiptera: Fulgoridae) traps. Journal of Economic Entomology, 116 (2 ), 426–434.36881709
Phillips, A. J. , & Dudík, M. (2008). Modeling of species distributions with Maxent new extensions. Ecography, 31 (2 ), 161–175.
Phillips, S. J. , Anderson, R. P. , & Schapire, R. E. (2006). Maximum entropy modeling of species geographic distributions. Ecological Modelling, 190 , 231–259.
Phillips, S. J. , Anderson, R. P. , Dudík, M. , Schapire, R. E. , & Blair, M. E. (2017). Opening the black box: An open‐source release of Maxent. Ecography, 40 (7 ), 887–893.
Radosavljevic, A. , Anderson, R. P. , & Araújo, M. (2014). Making better Maxent models of species distributions: Complexity, overfitting and evaluation. Journal of Biogeography, 41 (4 ), 629–643.
Ramasamy, M. , Das, B. , & Ramesh, R. (2021). Predicting climate change impacts on potential worldwide distribution of fall armyworm based on CMIP6 projections. Journal of Pest Science, 95 (2 ), 841–854.
Ramos, R. S. , Kumar, L. , Shabani, F. , & Picanco, M. C. (2018). Mapping global risk levels of Bemisia tabaci in areas of suitability for open field tomato cultivation under current and future climates. PLoS One, 13 , e0198925.29902221
Schuldt, A. , & Assmann, T. (2009). Environmental and historical effects on richness and endemism patterns of carabid beetles in the western Palaearctic. Ecography, 32 (5 ), 705–714.
Schuldt, A. , & Assmann, T. (2011). Belowground carabid beetle diversity in the western Palaearctic‐effects of history and climate on range‐restricted taxa (Coleoptera, Carabidae). Zookeys, 100 , 461–474.
Urban, J. M. (2020). Perspective: Shedding light on spotted lanternfly impacts in the USA. Pest Management Science, 76 (1 ), 10–17.31525270
Urban, J. M. , & Leach, H. (2023). Biology and management of the spotted lanternfly, Lycorma delicatula (Hemiptera: Fulgoridae), in the United States. Annual Review of Entomology, 68 (1 ), 151–167.
Velasco, J. A. , & González‐Salazar, C. (2019). Akaike information criterion should not be a “test” of geographical prediction accuracy in ecological niche modelling. Ecological Informatics, 51 , 25–32.
Wakie, T. T. , Neven, L. G. , Yee, W. L. , & Lu, Z. (2020). The establishment risk of L. delicatula (Hemiptera: Fulgoridae) in the United States and globally. Journal of Economic Entomology, 113 (1 ), 306–314.31579914
Wang, F. , Wang, D. , Guo, G. , Zhang, M. X. , Lang, J. Y. , & Wei, J. F. (2021). Potential distributions of the invasive barnacle scale Ceroplastes cirripediformis (Hemiptera: Coccidae) under climate change and implications for its management. Journal of Economic Entomology, 114 (1 ), 82–89.33184624
Wang, R. L. , Yang, H. , Wang, M. T. , Zhang, Z. , Huang, T. T. , Wen, G. , & Li, Q. (2020). Predictions of potential geographical distribution of Diaphorina citri (Kuwayama) in China under climate change scenarios. Scientific Reports, 10 (1 ), 9202.32513980
Wei, J. F. , Zhao, Q. , Zhao, W. Q. , & Zhang, H. F. (2018). Predicting the potential distributions of the invasive cycad scale Aulacaspis yasumatsui (Hemiptera: Diaspididae) under different climate change scenarios and the implications for management. Peer J, 6 , e4832.29844981
Wei, J. , Li, X. , Lu, Y. , Zhao, L. , Zhang, H. , & Zhao, Q. (2019). Modeling the potential global distribution of Phenacoccus madeirensis green under various climate change scenarios. Forests, 10 , 773.
Wilson, J. R. , Richardson, D. M. , Rouget, M. , Procheş, Ş. , Amis, M. A. , Henderson, L. , & Thuiller, W. (2007). Residence time and potential range: Crucial considerations in modelling plant invasions. Diversity and Distributions, 13 (1 ), 11–22.
Xiao, G. (1992). Forest insect of China (2nd ed.). Forest Research Institute.
Xu, D. P. , Li, X. Y. , Jin, Y. W. , Zhuo, Z. H. , Yang, H. J. , Hu, J. M. , & Wang, R. L. (2020). Influence of climatic factors on the potential distribution of pest Heortia vitessoides Moore in China. Global Ecology and Conservation, 23 , e01107.
Zhang, Y. B. , Zhao, Z. X. , Wang, Y. J. , & Liu, T. L. (2023). Suitable habitats for Cicadella viridis and Evacanthus interruptus (Hemiptera: Cicadellidae) with global climate change. Journal of Entomological Science, 58 (2 ), 215–229.
Zhang, Z. Z. (1993). Forest insects. China Forestry Publishing House.
Zhao, Z. X. , Feng, X. L. , Wang, Y. J. , Zhou, Z. X. , & Zhang, Y. B. (2024b). Potential suitability areas of Sitobion miscanthi in China based on the MaxEnt model: Implications for management. Crop Protection, 183 , 106755.
Zhao, Z. X. , Yang, L. , Long, J. K. , Chang, Z. M. , & Chen, X. S. (2024a). Predicting suitable areas for Metcalfa pruinosa (Hemiptera: Flatidae) under climate change and implications for management. Journal of Insect Science, 24 (3 ), 7; ieae053.
Zhou, Y. , Tao, J. , Yang, J. , Zong, S. , & Ge, X. (2023). Niche shifts and range expansions after invasions of two major pests: The Asian longhorned beetle and the citrus longhorned beetle. Pest Management Science, 79 (9 ), 3149–3158.37013934
Ziska, L. H. , Blumenthal, D. M. , Runion, G. B. , Hunt, E. R., Jr. , & Diaz‐Soltero, H. (2011). Invasive species and climate change: An agronomic perspective. Climatic Change, 105 (1–2 ), 13–42.
Zizka, A. , Silvestro, D. , Andermann, T. , Azevedo, J. , Duarte Ritter, C. , Edler, D. , Farooq, H. , Herdean, A. , Ariza, M. , Scharn, R. , Svantesson, S. W. , Zizka, V. , & Antonelli, A. (2019). Coordinate cleaner: Standardized cleaning of occurrence records from biological collection databases. Methods in Ecology and Evolution, 10 (5 ), 744–751.
