Supplementary material 2 from: Ruzzier E, Scaccini D, Tirozzi P, Orioli V, Dondina O, Di Giulio A, Pozzebon A, Bani L (2025) Predicting the global distribution and invasion scenarios of the Spotted Lanternfly, Lycorma delicatula (White, 1845) (Hemiptera, Fulgoridae). NeoBiota 103: 267-298. https://doi.org/10.3897/neobiota.103.154246

Predicting the global distribution and invasion scenarios of Lycorma delicatula (Hemiptera: Fulgoridae) Enrico Ruzzier 1,4, Davide Scaccini 2, Pietro Tirozzi 3, Valerio Orioli 3, Olivia Dondina 3,4, Andrea Di Giulio 1,4, Alberto Pozzebon 2 & Luciano Bani 3,4 1 Department of Science, Università Roma Tre, Rome, Italy 2 Department of Agronomy, Food, Natural resources, Animals, Environment (DAFNAE), University of Padua, Italy 3 Department of Earth and Environmental Sciences, University of Milano-Bicocca, Milan, Italy 4 National Biodiversity Future Center (NBFC), Palermo, Italy Email of the corresponding author: [email protected] SUPPLEMENTARY INFORMATION Fig. S1. Distribution of the 20 sets of PAs + BPs used for the habitat suitability model (HSM) calibration in the narrow native calibration range. In green urban-biased PAs (number of bias PAs = number of occurrences) and in black random BPs (number of random BPs = 4 times the number of occurrences); in red Lycorma delicatula occurrences. Fig. S2. Distribution of the 5 sets of BPs used for the bioclimatic suitability model (BSM) calibration in the wide calibration range. In green random BPs (number of random BPs = 5 times the number of occurrences); in red Lycorma delicatula occurrences. Fig. S3. Response curves of covariates included in the Lycorma delicatula habitat suitability models (HSM). Fig. S4. Response curves of covariates included in the Lycorma delicatula bioclimatic suitability models (BSM). a) b) c) Fig. S5. Habitat suitability map for Lycorma delicatula in the narrow native calibration range obtained projecting the habitat suitability model (HSM); a) lower limit of the 95% confidence interval; b) mean; c) upper limit of the 95% confidence interval. a) b) c) Fig. S6. Bioclimatic suitability map for Lycorma delicatula in the wide calibration range obtained projecting the bioclimatic suitability model (BSM); a) lower limit of the 95% confidence interval; b) mean; c) upper limit of the 95% confidence interval. Fig. S7. Overall suitability map for Lycorma delicatula in the narrow native calibration range. Fig. S8. Boyce of the overall suitability map for Lycorma delicatula in the narrow native calibration range. Fig.S15. Eastern and southeast Asia (ESA) projections of Lycorma delicatula suitability maps; a) habitat suitability map; b) bioclimatic suitability map. Occurrence points: blue: native range; red: invaded range. Fig.S16. USA projections of Lycorma delicatula suitability maps; a) habitat suitability map; b) bioclimatic suitability map. Red points: L. delicatula occurrences. Fig.S17. Details of USA (East Coast) projections of Lycorma delicatula suitability maps; a) habitat suitability map; b) bioclimatic suitability map. Red points: L. delicatula occurrences. Fig.S18. European projections of Lycorma delicatula suitability maps; a) habitat suitability map; b) bioclimatic suitability map. Red point: L. delicatula putative occurrence.