Research Article |
Corresponding author: Hector Zumbado-Ulate ( zumbadohector@gmail.com ) Academic editor: Darren Kriticos
© 2023 Hector Zumbado-Ulate, Tyler E. Schartel , Gregory S. Simmons, Matthew P. Daugherty.
This is an open access article distributed under the terms of the CC0 Public Domain Dedication.
Citation:
Zumbado-Ulate H, Schartel TE, Simmons GS, Daugherty MP (2023) Assessing the risk of invasion by a vineyard moth pest guild. NeoBiota 86: 169-191. https://doi.org/10.3897/neobiota.86.100579
|
Biological invasions are most effectively managed when identified in their early stages, which often hinges on robust surveillance programs. The recent invasion of the European grapevine moth (Lobesia botrana) in California suggests that viticultural areas in the western United States may face severe economic consequences from this and other Tortricid and Pyralid moth species if they were to establish. To gain insights into the risk these grapevine pests pose, we used occurrence records for L. botrana and four other moths native to Europe or the eastern United States and selected environmental variables to predict the extent of climatically suitable areas and potential pest co-occurrence along the West Coast of the United States. A suite of models was generated using MaxEnt with species-specific tuning of model settings. Overall, the results confirmed high suitability for L. botrana to establish across much of the study region, driven largely by high monthly variability in precipitation and low elevation. Two species were predicted to have intermediate suitability to establish over the study region (i.e., grape tortrix moth, Argyrotaenia ljungiana; grape berry moth, Paralobesia viteana), while two others had low suitability (i.e., European grape berry moth, Eupoecilia ambiguella; Christmas berry webworm, Cryptoblabes gnidiella). The highest predicted potential for co-occurrence was between L. botrana and P. viteana, accounting for 19% of the total viticulture area, followed by L. botrana and A. ljungiana for 11% of the study area. These results may help with the optimization of surveillance efforts by indicating which species or areas should be prioritized for the deployment of invasive pest detection programs with pheromone traps. Indeed, given the apparent potential for co-occurrence of multiple moth pests in certain areas, our results may inform where single or multi-lure traps should be deployed as a more cost-efficient monitoring tool.
detection trapping, invasion risk, pest surveillance, species distribution model
Invasive species are a significant threat to global biodiversity (
Biological invasions in the United States cost more than $100 billion annually and are increasing in frequency (
A recent invader of particular importance to California’s wine, raisin, and table grape industry is the European grapevine moth (EGVM), Lobesia botrana (Lepidoptera: Tortricidae) [Denis & Schiffermüller]; one of several Lepidopteran agricultural pests that have proven capable of rapid geographic range expansion (
In the United States, EGVM was first detected in California in late 2009 (
As part of a larger study evaluating the factors that contributed to the successful eradication of EGVM in California,
EGVM is only one of several moth pests of grapevines with the potential to be highly destructive. Other species in the families Tortricidae and Pyralidae have proven to be significant pests in other viticultural regions (
Monitoring using pheromone-baited traps is commonly employed for pest management, and for early detection of invasive insects (
We gathered occurrence records from the native and invaded ranges of five high-risk lepidopteran pests of grapevines and selected a number of environmental variables to quantify invasion risk along the West Coast of the United States (
We focused on the invasive potential of five grapevine pests (EGBM, EGVM, GBM, GTM, and CBW) in grape-growing regions along the West Coast of the United States, in portions of California, Oregon, and Washington (5–40°N, 70–118°W; Suppl. material
Risk assessments based on extrapolations of species distribution model (SDM) predictions to regions and periods different from the conditions used to calibrate the model (i.e., model transferability) are effective for pest management and species conservation planning (
We downloaded occurrence data for all five pest species in their native and other invaded ranges from 1960 to the present date from the Global Biodiversity Information Facility (GBIF; www.gbif.org; Suppl. material
All 19 BIOCLIM variables (
Here, we provide an overview of our climatic suitability modeling methodology following the ODMAP (Overview, Data, Model, Assessment, and Prediction) protocol for species distribution models (
Model selection was conducted using the highest average of the area under the curve of the receiver-operating characteristic (‘AUC mean’), the standardized true skill statistic (sTSS), and the average of the 10-percentile training omission rate (‘10.or.pt mean’). For selected models, we estimated the percent contribution of each selected abiotic predictor and generated response curves by comparing the probability of each pest species’ presence relative to each abiotic predictor (
To identify areas that may be susceptible to the establishment of multiple moth species, we used two alternative thresholds to generate binary predictions (raster absence-presence maps) of the potential range of each species across the study region. Binary predictions were transformed into polygons to quantify the extent of climatically suitable areas (ESH;
To identify areas of co-occurrence, binary predictions were combined into a single raster. Because the binary predictions only had values of zero and one, the resulting cumulative raster displayed values between zero (no pest species predicted to occur in a pixel) and five. Then, this cumulative raster was transformed into a polygon to estimate the potential co-occurrence of multiple species, from two to five species. To calculate ESH we transformed the projected coordinate system of binary predictions and AVAs from WGS84 to NAD 1983 Albers contiguous USA (ESRI 102003).
Finally, a principal component analysis (PCA) was generated to visualize the environmental space where multiple pest species are predicted to overlap. For this, we selected the ten predictors with highest contributions to the SDMs (Table
Percent contribution (% C) and permutation importance (% P) of selected abiotic predictors in species-specific climatic suitability models (EGBM = European grape berry moth, Eupoecilia ambiguella; EGVM = European grapevine moth, Lobesia botrana; GBM = Grape berry moth, Paralobesia viteana; GTM = Grape tortrix moth, Argyrotaenia ljungiana; CBW = Christmas berry webworm, Cryptoblabes gnidiella). The two predictors with the highest contributions are in bold.
Predictor | EGBM | EGVM | GBM | GTM | CBW |
---|---|---|---|---|---|
% C (% P) | % C (% P) | % C (% P) | % C (% P) | % C (% P) | |
BIO2; Mean diurnal range (°C) | 2.7 (0.5) | – | 10.5 (25.6) | 32.2 (47) | 69.5 (20.1) |
BIO3; Isothermality (°C) | 7.7 (6.8) | 6.1 (5.3) | – | 35.5 (35.6) | – |
BIO7; Temperature annual range (°C) | – | 11.4 (13.9) | – | – | – |
BIO8; Mean temperature of wettest quarter (°C) | 6 (7.8) | 3.9 (5.3) | – | 0 (0) | – |
BIO13; Precipitation of wettest month (mm) | 7.1 (14.7) | 6.1 (7.2) | – | – | – |
BIO14; Precipitation of driest month (mm) | – | – | 0 (0) | – | 0.7 (7.4) |
BIO15; Precipitation seasonality (mm) | 2.2 (5.9) | 57.6 (42.1) | – | 13.5 (2.4) | 0.5 (3.7) |
BIO18; Precipitation of warmest quarter (mm) | – | 0 (0) | 27.1 (41.2) | 0.2 (0.4) | 0 (0.2) |
BIO19; Precipitation of coldest quarter (mm) | 3.8 (0) | – | – | – | 1.2 (0.1) |
Elevation (m) | 57.6 (57.6) | 12.6 (10) | 0 (0) | 14.4 (1.6) | 7.9 (6.1) |
EPQ; Emberger’s pluviothermic quotient | – | – | – | 0.3 (2.1) | – |
gDD5; growingDegDays5 (°C)1 | 4.4 (5.1) | – | – | – | – |
mTW; Minimum temperature of warmest month (°C) | – | 0.3 (1.3) | 56.5 (23) | 3.8 (10.5) | 12 (0) |
PETDQ; PET of driest quarter (mm)2 | – | 1.4 (14.1) | 3.7 (2.7) | 0.1 (0.4) | 3.4 (42.1) |
PETS; PET seasonality (mm)2 | 8.5 (1.6) | 0.6 (0.9) | 2.2 (7.4) | – | 4 (18.2) |
PETWQ; PET of wettest quarter (mm)2 | – | – | – | – | 0.9 (2.1) |
One preferred model was identified for each species (Suppl. material
Sixteen abiotic predictors were retained among the five models selected (Table
Of the five focal species evaluated, the area estimated to be moderate to highly climatically suitable for establishment was highest for EGVM, especially along the coastline, in western and central areas of Washington and Oregon, and central regions of California (Fig.
Extent of climatically suitable areas (ESH) and corresponding percent of the total area of viticulture regions (% VR) for five moth species (EGBM = European grape berry moth, Eupoecilia ambiguella; EGVM = European grapevine moth, Lobesia botrana; GBM = grape berry moth, Paralobesia viteana; GTM = grape tortrix moth, Argyrotaenia ljungiana; CBW = Christmas berry webworm, Cryptoblabes gnidiella) using two binary thresholds for suitability: 10-percentile lowest omission rate threshold (10.or.pt) and the maximum training sensitivity plus specificity threshold (maxSS).
Species | Threshold | |||
---|---|---|---|---|
10.or.pt | maxSS | |||
ESH (km2) | % VR | ESH (km2) | % VR | |
EGBM | 21640 | 12.1 | 1134.1 | 0.6 |
EGVM | 115605 | 64.6 | 64894.6 | 36.3 |
GBM | 36940 | 20.6 | 64.8 | 0.0 |
GTM | 22776 | 12.7 | 109.6 | 0.1 |
CBW | 1074 | 0.6 | 123.0 | 0.1 |
Number of species1 | ESH (km2) | % VR | ESH (km2) | % VR |
0 | 37917.8 | 21.2 | 113329.8 | 63.3 |
1 | 95758.2 | 53.5 | 65396.1 | 36.6 |
2 | 34154.6 | 19.1 | 115.5 | 0.1 |
3 | 9586.5 | 5.4 | 79.6 | 0.0 |
4 | 1402.6 | 0.8 | 0.0 | 0.0 |
5 | 101.3 | 0.1 | 0.0 | 0.0 |
Climatic suitability map for the European grapevine moth (EGVM), Lobesia botrana, in viticultural regions along the West Coast of the United States A continuous climatic suitability estimates B binary predictions of climatically suitable areas based on the 10-percentile lowest omission rate threshold.
Climatic suitability was relatively moderate throughout the study region for two pest species. For GBM, the most climatically suitable regions occurred in small patches across viticulture areas in Washington and Oregon, and very small portions along the coast of California (Fig.
Climatic suitability estimates and binary predictions based on the 10-percential lowest omission rate threshold for four moth species A, B grape berry moth, Paralobesia viteana C, D grape tortrix moth (GTM), Argyrotaenia ljungiana E, F European grape berry moth (EGBM), Eupoecilia ambiguella G, H Christmas berry webworm (CBW), Cryptoblabes gnidiella.
Finally, the vast majority of the study region was projected to have relatively low climatic suitability for two focal species. For EGBM, the most climatically suitable regions for establishment occurred in central Washington and far northern Oregon (Fig.
Implementing the 10.or.pt binary threshold revealed that approximately 25% of the overall area of viticulture regions was predicted to be climatically suitable for pest co-occurrence (Fig.
Regions predicted to be climatically suitable for multiple species A occurrence or co-occurrence of up to five of the moth species based on 10-percentile lowest omission rate thresholds B regions of potential co-occurrence of the three species with the highest overall suitability in the study region: European green vine moth (EGVM) Lobesia botrana, grape berry moth (GBM), Paralobesia viteana, and grape tortrix moth (GTM), Argyrotaenia ljungiana C principal component analysis depicting climatic envelopes and environmental predictors (Table
Traditionally, researchers have concentrated on understanding biological invasions at large scales (
Once invasive species settle, effective management becomes difficult and costly to carry out, especially in regions affected by multiple pest species (
Scientists have linked successful invasive species to high abundance, wide distribution in their native ranges, and distinct traits that ease establishment and spread (e.g.,
Our results are consistent with those of a prior global analysis of EGVM suitability (
The prior invasion of EGVM in California, and its ultimate eradication, may yield important lessons for future responses to invasive insects (
We found that in addition to EGVM, two other moth species, GBM and GTM, might find moderate expanses of climatically suitable areas. Since the expansion of grape-growing regions in North America, GBM has increased in abundance and distribution, but it has not been observed in the western United States. However, viticulture regions in Washington and Oregon seem to offer suitable conditions, as this species performs well in temperatures in seasonal humid environments (
An additional benefit of this work was to model multiple species’ pest risk establishment probabilities which can also aid efforts to develop tools, such as multi-lure pheromone traps, that can simultaneously monitor multiple species while reducing the costs and time-intensive nature of monitoring efforts (
Although a growing number of studies have modeled the distribution of multiple species or assessed co-occurrence of multiple species through joint SDMs or occupancy models (
Detection and surveillance efforts are important components of early pest management strategies but are often costly and time-intensive (
Our findings suggest that most resources should be used to avoid a secondary spread of EGVM in the viticulture regions of the West Coast of the United States. Additionally, given the apparent potential for coexistence of the European species in some areas (
Thanks to M. Cooper, L. Varela and other members of the EGVM Technical Working Group for their input during the early phases of this research, and to R. Broadway for providing access to EGVM trapping records. We also thank K. Bigsby and A. Tripodi at the Animal and Plant Health Inspection Service (APHIS) for constructive comments that helped to improve the manuscript. This work was supported by funds from Animal and Plant Health Inspection Service, cooperative agreements AP20PPQS&T00C192 and AP21PPQS&T00C166 to M. Daugherty.
Assessing the risk of invasion by a vineyard moth pest guild
Data type: figures and tables (word document)
Explanation note: table S1. 34 environmental variables considered during the development of species distribution models for five lepidopteran grapevine pests. table S2. Climatic suitability modeling methods as organized by the ODMAP framework (