Phenological Mapping of Invasive Insects: Comparison
View Latest Version
Please note this is a comparison between Version 2 by Lindsay Dong and Version 1 by Brittany Barker.

Phenological maps can depict the development and seasonal activities (phenology) of invasive insects at area-wide scales, such as counties, states, or entire nations. When regularly updated using real-time and forecast climate data, these maps may improve the timeliness of early detection and control tactics that target specific life stages. Rapid responses to invasive insects may increase the likelihood that populations are eradicated or controlled before they can spread or increase in size. Phenological maps may also be used to assess pest establishment risk, investigate pest–host interactions, and measure climate-driven changes in pest phenology.

  • degree-day model
  • forecast
  • pest
  • climate

1. Introduction

The protection of agricultural and natural resources depends on the precise timing of surveillance, monitoring, and management of invasive insect populations [1,2,3][1][2][3]. Scheduling pest control tactics according to calendar dates and expectations of the “normal” time in which seasonal activities (phenology) of pests have occurred in previous years is often ineffective because rates of insect development often vary annually due to variations in weather [4,5,6,7][4][5][6][7]. Conversely, modeling the phenology of an insect species using information on their life cycle requirements and climate data for the current year may increase the precision of estimates of dates when important seasonal events occur, such as the first adult emergence and egg hatch [4,7,8][4][7][8]. This information may improve the effectiveness and cost efficiency of early detection and management programs because these programs often target life stages that are more observable (e.g., larvae vs. adults of wood-boring beetles) or more vulnerable to control tactics, such as pesticide treatments [3,9,10,11,12,13][3][9][10][11][12][13].
Degree-day models are widely used in decision support systems that predict the phenology of agricultural insect pests because of their simplicity and ability to accommodate multiple species with varying life histories [8,9,10,14,15,16,17,18,19,20][8][9][10][14][15][16][17][18][19][20]. The development and phenology of an organism in a degree-day model is driven by heat accumulation above a lower temperature threshold (and oftentimes below an upper temperature threshold) over a daily or weekly time step (Figure 1) [4,5,7,21,22,23,24][4][5][7][21][22][23][24]. The lower developmental temperature is often referred to as the base temperature. Typical examples of degree-day models that have been used for many years include those that predict first egg hatch of the codling moth [Cydia pomonella (L.)] in tree fruits [13,25,26[13][25][26][27][28],27,28], first emergence of the western cherry fruit fly (Rhagoletis indifferens Curran) [29], and adult flight, egg hatch, and larval development of the spongy (formerly “gypsy”) moth [Lymantria dispar (L.)] [30,31,32,33,34][30][31][32][33][34].
Insects 15 00006 g001
Figure 1. Degree-day (DD) accumulations for 1 January to 1 September predicted by site-based and spatial phenology models. (A) Plot of date vs. DD accumulations depicts dates of adult flight activity at a single site. (B) An elevation map with predictions for multiple sites is shown with a key of DD requirements for adult flight activity. (C) Phenological map of the same area as (B) shows continuous predictions. The X and Y axis in (B,C) indicate longitude and latitude, respectively.
Most degree-day models for invasive insect species use climate data for a specific site, such as a weather station. Site-based model predictions, such as degree-day accumulations and dates of phenological events (e.g., first adult flight), are typically displayed in tabular and/or graphical formats (Figure 1A). This information is useful for decision support for pest surveillance or management in small areas, such as a forest parcel, fruit orchard, or vineyard [3,9,10,13,35][3][9][10][13][35]. However, site-based models are less applicable for decision support at area-wide scales, such as large counties or states, because phenology may vary spatially due to geographic factors (e.g., latitude, elevation, and continental effects), anthropogenic disturbances, and biological factors, such as spatial variations in population development and host plant availability.
Model predictions for multiple sites (e.g., in the form of number values, color-coded pins, or other expressions) can be overlayed on a base map to visualize geographic variation in insect phenology [8]. For example, degree accumulations for each site can be shown as text labels on an elevation map, which can be referenced to a key of degree-day requirements for important phenological events (Figure 1B). Nonetheless, maps with discontinuous model predictions may be of limited use for decision-makers in areas that lack model predictions, particularly for topographically and climatically complex areas where phenology can vary over short distances.
The continuous mapping of phenology (hereafter phenological mapping) depicts model predictions over an entire geographic area (Figure 1C). Phenological mapping of pests became more common beginning in the 1990s with advancements in computers and a growing number of digitized geographic datasets and geographic information systems (GIS) software options [32,39–42]. Digital spatial climate datasets developed over the past ca. 30 years have increased the robustness and timeliness of phenological maps. 

2. Data Requirements for Degree-Day Models

Detailed reviews on data requirements and methods for degree-day modeling of insects already exist [17,23,24,48,49][17][23][24][36][37] and are, therefore, only summarized here. Degree-day models are often developed using experimentally collected data on temperature–development relationships to estimate parameters, such as developmental rates, developmental temperature thresholds, duration of life stages, and stage-specific events [17,50,51,52][17][38][39][40]. These data are combined with daily temperature data, typically minimum and maximum temperatures (Tmin and Tmax, respectively), to estimate degree-days using various calculation methods [4,7,14,24,50,53][4][7][14][24][38][41]. A start date or biological event, such as the first flight, is needed to synchronize the insect phenology model to field populations [24,54,55][24][42][43].

3. Types of Phenological (Degree-Day) Maps

3.1. Generic Degree-Day Map

Generic degree-day maps show the current degree-day accumulations based on one or more standard lower temperature thresholds. For example, degree-day maps at Michigan State University’s Enviroweather [https://www.enviroweather.msu.edu (accessed on 19 December 2023)] use a standard lower temperature threshold of 50 °F (10 °C), and those at USPest.org [https://uspest.org/wea/index.html#DDMAPS/ (accessed on 19 December 2023)] use thresholds of 32, 41, or 50 °F (0, 5, or 10 °C, respectively; Figure 2). Degree-day maps at the University of Wisconsin’s AgWeather Vegetable Disease & Insect Forecasting Network [https://agweather.cals.wisc.edu/vdifn?p=insect (accessed on 19 December 2023)] use several different species-specific thresholds.
Insects 15 00006 g002
Figure 2. Example cumulative degree-day (DD) mapping products for Washington state generated online at USPest.org [https://uspest.org (accessed on 19 December 2023)]. All maps depict DD accumulation [base 41 °F (5 °C)] between 1 January and 31 August of 2020 (calculated using the single triangle method). (A) Map generated using an older version of the program [https://uspest.org/cgi-bin/usmapmaker.pl (accessed on 19 December 2023)] Awith a spatial resolution of 800 m. Black dots indicate the locations of weather stations used in the correction of monthly PRISM-based cumulative DDs. (B) Map produced by the newer program [https://uspest.org/dd/mapper (accessed on 19 December 2023)] that uses daily PRISM data with a coarser (4 km) spatial resolution (no weather station correction required). (C) Downscaled output produced using the newer program (800 m resolution).
Generic degree-day maps provide a good, general reference to show how the season is progressing, especially when compared to averages, such as 30-year normal degree-day maps. However, predictions of current degree-day accumulations are not matched up with insect life stages, so they must be used with experience and care for guiding surveying and management activities. 
Enviroweather’s degree-day models use 1 March as the start date and includes only Michigan and Wisconsin maps, whereas models at USPest.org at https://uspest.org/ (accessed on 20 December 2023) uses 1 January as the start date and covers all of the contiguous United States (CONUS) with separate mapping of major subregions (e.g., Midwest, Northwest, Southeast) and of states in the Pacific Northwest. Both websites are updated daily to produce maps based on climate data for the current year and historical averages (30-year normals). An additional map depicts the difference between historical and current year degree-day accumulations. 
The degree-day mapping program at USPest.org [https://uspest.org/cgi-bin/usmapmaker.pl (accessed on 19 December 2023)] first appeared online in 1998 as a decision support tool for pest management in Oregon. The mapping region was expanded to include the entire Pacific Northwest region by 2002, and then to include CONUS at a higher spatial resolution (800 m) by 2005 [56,57][44][45]. Degree-days based on 30-year normals (centered on 1995) are calculated using gridded monthly temperature data from the PRISM database, while daily real-time degree-days are calculated using data from thousands of weather stations in the USPest.org collection of public networks. For the map creation, the degree-day mapping program uses a ‘climatologically aided interpolation’ method (sometimes more generally referred to as ‘delta correction’) that uses a gridded climate dataset, such as PRISM, to improve the interpolation of a site-based dataset, such as recent station observations (Figure 2A) [45,58][46][47]. 
Constructed in 2017, a second custom online phenology mapping program at USPest.org [https://uspest.org/dd/mapper (accessed on 19 December 2023)] was developed to offer an alternative and simpler workflow that uses real-time daily PRISM temperature grids and does not require data from multiple weather stations for correction (Figure 2B). However, the input data has a lower resolution (4 km), which, although adequate for most state-level maps, would be insufficient for small states or single growing regions that are topographically complex, such as Hood River County, Oregon. To address this issue, this version of the degree-day mapping program includes an option to downscale resulting degree-day maps to 800 m (Figure 2C) using a custom distance–elevation weighted regression algorithm.
Insects 15 00006 g003
Figure 3. Degree-day lookup table map for the brown marmorated stink bug [Halyomorpha halys (Stål)] produced as part of the Spatial Analytic Framework for Advanced Risk Information Systems (SAFARIS) Field Operations Weekly map series [59][48]. The map uses a degree-day lookup table to associate cumulative degree-days with predicted life stages present across the contiguous United States on 22 May 2023. Thus, it provides a “snapshot in time” of phenology model predictions for a specific date. Reproduced with permission from SAFARIS, Brown Marmorated Stink Bug (Halymorpha halys) Phenological Stages. Published by SAFARIS, U.S. Department of Agriculture (USDA), and North Carolina State University, 2023. Available online: https://safaris.cipm.info (accessed on 19 December 2023).

3.2. Degree-Day Lookup Table Map

Degree-day lookup table maps show the current life stages or phenological events of an organism that correspond to specified values or ranges of accumulated degree-days for a specified date (Figure 3). Thus, degree-day accumulations, which are depicted in generic degree-day maps, are matched to specific points or events during the life cycle. For insects, life cycle points (and events) could typically include the egg stage present, egg hatch, larval stage present, pupal stage present, adult emergence and presence, and egg laying. The simplicity of the approach and its applicability to multiple organisms has sustained its use for several years.

Specific features that make the degree-day lookup table approach so common include:
  • A relatively straightforward workflow. The workflow of generating a degree-day lookup table map involves using gridded daily Tmin and Tmax data to calculate degree-day accumulations between a start date (usually 1 January, although some models use other start dates, such as 1 March) and a specified end date. Degree-day lookup tables are then used to associate degree-day accumulations with life stages, and output maps depict the results with color tables and legends.
  • The use of common base thresholds for multiple species. As degree-day lookup table maps are relatively simple and generic, there is the potential to use the same lower temperature threshold base maps for multiple species. This contrasts with more complex models that would require a separate base map for each case because of the application of different parameter values, including lower and upper thresholds, calculation methods, start dates, or diapause. The use of upper thresholds is relatively rare, at least for degree-day lookup table maps. For example, the SAFARIS FO Weekly maps produced by models for the old world bollworm [Helicoverpa armigera (Hübner)] and brown marmorated stink bug [Halyomorpha halys (Stål)] are constructed using the same 54 °F base maps, with no upper threshold. 
  • An ability to provide a “snapshot in time” for a single date. This allows, for example, regular updates that provide a gradually changing view of the current or near-future status of insect phenology. For example, degree-day lookup table maps produced by the Degree-Day, establishment Risk, and Phenological event maps (DDRP) platform [16] every 2−3 days depict the life stage and generation of insects on the map issue date. The USA National Phenology Network’s Pheno Forecast maps take advantage of 7-day National Digital Forecast Database (NDFD) forecasts to provide a 1-week “look ahead” prediction for CONUS [47][49]. SAFARIS PestCAST maps include a 1-month forecast using a 7-day NDFD forecast followed by three weeks of recent 20-year average PRISM data [15].
  • Relatively simple design requirements. Degree-day lookup table maps can be designed as very simple visualization tools, such as by designing legend items to display only the stage or activity of interest (e.g., adult flight or egg hatch). Other stages and activities can then be represented as merged entries. The practice of focusing end users on a single target event represents a clear trade-off in reducing complexity (of multiple life stages) for users who may need clear directions in implementing surveillance or management actions.

3.3. Phenological Event Map

In contrast to a degree-day lookup table map, a phenological event map depicts the dates on which accumulating degree-days reach a value (target degree-day total) that corresponds with a selected phenological event [16,32,40,61][16][32][50][51]. Phenological event maps may offer the following advantages over degree-day lookup table maps to support monitoring and surveillance programs for invasive insects:
  • Standardization. Mapping dates of phenological events allows for the standardization of legends and color tables across multiple species and events. For example, the colors assigned to each range of dates in the legend (e.g., 1–8 January = dark blue, 9–16 January = medium blue, etc.) can be applied to several events within a species, as well the same or different events in other species. As an example, Figure 4 shows a phenological event map produced by DDRP that depicts the average date of egg laying by the overwintering generation of the light brown apple moth [Epiphyas postvittana (Walker)] for 2023.
  • Operationally ready. Phenological event maps could be considered a more operational (tactical) product than degree-day lookup tables because they predict dates of events for a particular life stage, potentially up to weeks or months into the future.
  • Simpler comparisons and expression of error rates. Phenological event maps allow for more direct comparisons of year-to-year variations of events than generic degree-day and degree-day lookup table maps. It is a relatively simple recordkeeping and reporting exercise to express differences in dates in the form of days difference [12,27,28][12][27][28]. This approach can also be used to express errors between predicted and observed events as discussed below (“8. Model validation”).
Insects 15 00006 g006
Figure 4. Phenological event map depicting the earliest date of egg laying by the overwintered (OW) generation of the light brown apple moth (Euphyas postvittana) for 2023 produced by DDRP (map issue date: 13 December 2023).

4. Applications of Phenological Maps

Phenological maps produced using real-time and forecast climate data have the potential to support the early detection of invasive pests. For example, regularly updated phenological maps on the SAFARIS platform provide decision support for the Cooperative Agricultural Pest Survey (CAPS) program [15[15][16],16], which conducts national and statewide surveys for exotic plant pests in the United States deemed to be of regulatory significance to the United States Department of Agriculture (USDA) Animal and Plant Health Inspection Service’s (APHIS) Plant Protection and Quarantine (PPQ) program [64][52].

Phenological maps can support pest managers in timing treatments or other control tactics that target certain life stages. For example, phenological maps of egg hatch and larval development for the spongy moth were developed to support the timing of insecticidal sprays conducted for “stop the spread” programs in the eastern United States [30,31,32][30][31][32].  

Maps that depict a pest’s potential number of generations per year (i.e., voltinism) may help identify areas at risk of establishment because persistence in a new area requires a life cycle completion [15,65,66,67][15][53][54][55]. Additionally, maps of voltinism provide insight into expected levels of pest growth and subsequent damages to host plants for a given year [65,68,69][53][56][57]. Phenological maps may be used to assess the potential impacts of climate change on invasive insect phenology. For example, some studies have modeled the phenology of invasive insects under future climate change scenarios to estimate the impacts of global warming and altered precipitation patterns on important phenological events, potential voltinism, and population growth [42,68–72]. Mapping the extent to which critical time periods of insect life cycles coincide with phenological windows of host plant suitability may provide insight into pest establishment risk and outbreaks dynamics [34,73–76]. 

5. Gridded Climate Data

Maps used for within-season decision support of invasive insects depend on having access to real-time daily Tmin and Tmax data with spatial resolutions that are appropriate for the needs of decision-makers. For example, phenological maps at a 4 km resolution are generally sufficient to support pest surveillance programs for the entire CONUS [15], but are probably not appropriate for smaller scales, such as a county or city. Real-time PRISM data with a spatial resolution of 4 km are freely available, and higher resolution (800 m) data can be purchased from the PRISM group. Real-time DDRP forecasts at USPest.org are produced using PRISM data (4 km resolution) as climatic inputs, whereas monthly updated North America Multi-Model Ensemble (NMME) 7-month forecasts or recent 10-year average PRISM data (calculated on a bimonthly basis) are used to predict pest phenology up to the end of the year [16].

Phenological mapping for within-season decision support in areas outside of the United States is typically hindered by a lack of real-time gridded daily Tmin and Tmax data. However, historical datasets may be used for model development and validation, such as those for Europe [83][58], continental North America and Hawaii [84[59][60],85], Brazil [79[61][62],86], China [80,81][63][64], India [82][65], and Bangladesh, Nepal, and Pakistan [78][66]. 

Some phenological mapping studies overcame an absence of readily available gridded daily climate data by interpolating weather station data over a landscape of interest using custom software [31,40,41,87,88,89,90][31][50][67][68][69][70][71]. For example, the GEO-BUG platform offered four automated interpolation methods to map the date at which a pest insect species reached a specified life stage in the United Kingdom [41,88][67][69]. Interpolation methods commonly applied to Tmin and Tmax estimates include those based on distance analyses (e.g., inverse distance weighted and spline interpolation) or geostatistics (e.g., kriging and multiple regression) [33,58,88,89,90,91][33][47][69][70][71][72].

6. Potential Sources of Error and Uncertainty

Common sources of error in insect phenology models include natural population variability, microclimatic factors, anthropogenic disturbances (e.g., land use patterns), and biotic factors (e.g., migration, host quality, competition, predation, and disease) [17,23,40,48,52,87,88][17][23][36][40][50][68][69].
Phenology models should somehow, but seldom do, include estimates of monitoring or sampling errors [93][73]. The best that can often be accomplished is to rely on observations collected by researchers or that have been verified by multiple people, such as “Research Grade” observations in the iNaturalist database [https://www.inaturalist.org (accessed on 19 December 2023)]. However, there is no guarantee that sampling errors can be minimized even when using published or verified observations [94][74].
It should be noted, however, that model and sampling errors are often lower than errors resulting from intrinsic population variability. Many species populations, for reasons including selection pressures due to climate variability in time and space, are spread widely phenologically [23,95][23][75]. Sometimes, this is referred to as seasonal plasticity [96][76].

7. Increasing Model Realism While Maintaining Simplicity

A common question in phenology modeling discussions is whether to use a simple model that is adaptable to multiple species versus a more complex, single-species model that can potentially deliver greater realism and accuracy [17,49,50,52,55][17][37][38][40][43]. Using simple linear degree-day accumulations to display either developmental stages (via lookup table) or dates of phenological events are often adequate in their predictive accuracy, while still simple enough to be readily adapted for dozens, or even hundreds, of species. For example, validation analyses of certain models used by the DDRP, Pheno Forecast, and SAFARIS platforms have revealed evidence of overall good predictive performances [15,16,46,47][15][16][49][77]. A population-based approach to phenological mapping may improve predictive accuracy because it helps account for developmental variation that naturally occurs within insect populations [98,99,100,101,102,103][78][79][80][81][82][83]. DDRP is a population modeling platform [16] that is most similar to the still relevant and still in use “grandfather” of phenology modeling platforms, Predictive Extension Timing Estimator (PETE) [8,27][8][27]. Both platforms use a cohort approach to population modeling that involves tracking the development of population cohorts through all life stages over the year using a daily time step. Cohorts may start development at different times, which produces a distribution of times in which they transition into a new life stage (e.g., egg to larva) or undergo a particular phenological event [99,103][79][83]. For many, if not most, insects in temperate climates, the use of single-factor degree-day models (i.e., temperature) have proven adequate as evidenced by their successful implementation for many agricultural pest species [9,10,15,20][9][10][15][20]. However, multifactor models that include additional driving variables, such as chilling, day length, and water availability, may be appropriate for certain species [48,61,70,73,92,104,105,106,107,108,109][36][51][84][85][86][87][88][89][90][91][92]. As an example, the brown marmorated stink bug, an important pest that has recently invaded and spread across the United States, has a reproductive diapause that terminates in mid-spring in response to photoperiod [114][93]. Modeling of this pest would likely be improved by using a two-phases approach in which photoperiod-cued diapause termination is followed by a degree-day response phase for spring and summer development [106,107][89][90]. Models that account for geographic variation in life history traits may be necessary to accurately predict phenology of insects in new environments [23,52,105,119][23][40][88][94]. For example, Grevstad et al. [105][88] used different values for the critical day length parameter in spatialized phenology models for the Japanese knotweed psyllid owing to genetically based differences in this species’ critical photoperiod in the native range. However, defining optimal parameter values would create the need for laboratory studies of insects sampled from across the species’ distribution, which may be infeasible. Even if these data were available for the entire known distribution, they may not be useful for newly introduced species that have unknown geographic origins [61][51]. Therefore, it is most prudent, with this lack of full understanding of the life cycle, to model what is understood and simply state the model’s shortcomings so that they can be considered by managers.

8. Model Validation

Errors in predictions of phenology models should be estimated when validation data are available [20,55,90][20][43][71]. Model overprediction, in which events or life stages are predicted later than observed dates, is typically more problematic than underprediction because decision-makers may miss the best opportunity to detect or manage populations. Communicating potential model errors to end users of phenological maps may allow them to adjust the timing of their surveillance and management activities accordingly [55][43]. Phenology model validation requires a set of observations not used in model development. Observations could be resampled using jackknifing or cross-validation methods, split into separate sets (e.g., 75% for development, 25% for validation), or originate from different years or areas [31,87,90,120,121][31][68][71][95][96]. Potential sources of data for model validation include peer-reviewed literature, reports, graduate theses, unpublished monitoring studies, and online databases such as iNaturalist, Nature’s Notebook [https://www.usanpn.org/natures_notebook (accessed on 19 December 2023)], and the National Agricultural Pest Information System [https://cmr.earthdata.nasa.gov/search/concepts/C1214608223-SCIOPS (accessed on 19 December 2023)] [15,16,46,47][15][16][49][77].  Validating a spatial phenology involves extracting gridded model predictions for each georeferenced observation using GIS software. Ideally, the geographic precision of a phenological observation should be equal to or greater than the spatial resolution of climatic datasets used for modeling. For CONUS, this would correspond to spatial resolutions of 800 m to 4 km depending on the climatic dataset. The statistics for phenological model validation, keeping in line with the principle of parsimony, should be done using readily understood terms. Summary statistics may include the mean difference between predicted and observed dates (day of year, DOY), along with variability estimates, such as 95% confidence intervals, to provide an estimate of the overall error including the bias [55,87,90][43][68][71]. The bias may be estimated as the average amount by which predicted DOYs are greater than observed DOYs, whereas the mean absolute error (i.e., the average absolute difference) can estimate the expected number of days a given prediction might be in error. Model performance may also be evaluated using statistical tests, such as confusion matrices or equivalence tests [46,47,90]. 

9. Conclusions

Phenological maps can provide insight into the development and seasonal activities of invasive insects at area-wide scales, such as counties, states, or entire nations. Several web-based platforms offer generic degree-day maps, degree-day lookup table maps, and phenological event maps to support within-season decision-making for the detection and control of invasive insects for CONUS. With the development of real-time gridded climate datasets for more regions of the world, phenological maps could become more commonly used for planning pest surveillance and management activities. Phenological maps may also be used to assess establishment risk and to investigate pest–host interactions and climate-driven changes in pest phenology.

References

  1. Bowers, J.H.; Malayer, J.R.; Martínez-López, B.; LaForest, J.; Bargeron, C.; Neeley, A.D.; Coop, L.; Barker, B.S.; Mastin, A.J.; Parnell, S.; et al. Surveillance for early detection of high-consequence pests and pathogens. In Tactical Sciences for Biosecurity of Animal and Plant Systems; Cardwell, K.F., Bailey, K.L., Eds.; IGI Global: Hershey, PN, USA, 2022; pp. 120–177.
  2. Reaser, J.K.; Veatch, S.D.; Burgiel, S.W.; Kirkey, J.; Brantley, K.A.; Burgos-Rodri, J. The early detection of and rapid response (EDRR) to invasive species: A conceptual framework and federal capacities assessment. Biol. Invasions 2020, 22, 1–19.
  3. Vänninen, I. Advances in Insect Pest and Disease Monitoring and Forecasting in Horticulture; Burleigh Dodds Science Publishing: Sawston, UK, 2022.
  4. Wilson, L.T.; Barnett, W.W. Degree-days: An aid in crop and pest management. Calif. Agric. 1983, 37, 4–7.
  5. Herms, D.A. Using degree-days and plant phenology to predict pest activity. In IPM (Integrated Pest Management) of Midwest Landscapes, Minnesota Agricultural Experiment Station Publication SB-07645; Krischik, V., Davidson, J., Eds.; Minnesota Agricultural Experiment Station: St. Paul, MN, USA, 2004; pp. 49–59.
  6. Ascerno, M.E. Insect phenology and integrated pest management. J. Arboric. 1991, 17, 13–15.
  7. Zalom, F.G.; Goodell, P.B.; Wilson, L.T.; Barnett, W.W.; Bentley, W.J. Degree-Days: The Calculation and Use of Heat Units in Pest Management; University of California, Division of Agriculture and Natural Resources: Berkeley, CA, USA, 1983.
  8. Gage, S.H.; Whalon, M.E.; Miller, D.J. Pest event scheduling system for biological monitoring and pest management. Environ. Entomol. 1982, 11, 1127–1133.
  9. Ferguson, A.W.; Skellern, M.P.; Johnen, A.; Richthofen, V.; Watts, N.P.; Bardsley, E.; Murray, A.; Cook, S.M. The potential of decision support systems to improve risk assessment for pollen beetle management in winter oilseed rape. Pest Manag. Sci. 2016, 72, 609–617.
  10. Jones, V.P.; Brunner, J.F.; Grove, G.G.; Petit, B.; Tangren, G.V.; Jones, W.E. A web-based decision support system to enhance IPM programs in Washington tree fruit. Pest Manag. Sci. 2010, 66, 587–595.
  11. Rossi, V.; Sperandio, G.; Caffi, T.; Simonetto, A.; Gilioli, G. Critical success factors for the adoption of decision tools in IPM. Agronomy 2019, 9, 710.
  12. Cormier, D.; Chouinard, G.; Pelletier, F.; Vanoosthuyse, F.; Joannin, R. An interactive model to predict codling moth development and insecticide application effectiveness. IOBC-WPRS Bull. 2016, 112, 65–70.
  13. Jones, V.P. Using phenology models to estimate insecticide effects on population dynamics: Examples from codling moth and obliquebanded leafroller. Pest Manag. Sci. 2021, 77, 1081–1093.
  14. Roltsch, W.J.; Zalom, A.G.; Strand, J.F.; Pitcairn, M.J.; Frank, J.R.; Ann, G.Z.; Strand, J.F.; Pitcairn, M.J. Evaluation of several degree-day estimation methods in California climates. Int. J. Biometeorol. 1999, 42, 169–176.
  15. Takeuchi, Y.; Tripodi, A.; Montgomery, K. SAFARIS: A spatial analytic framework for pest forecast systems. Front. Insect Sci. 2023, 3, 1198355.
  16. Barker, B.S.; Coop, L.; Wepprich, T.; Grevstad, F.S.; Cook, G. DDRP: Real-time phenology and climatic suitability modeling of invasive insects. PLoS ONE 2020, 15, e0244005.
  17. Coop, L.B.; Barker, B.S. Advances in understanding species ecology: Phenological and life cycle modeling of insect pests. In Integrated Management of Insect Pests: Current and Future Developments; Kogan, M., Heinrichs, E., Eds.; Burleigh Dodds Science Publishing: Sawston, UK, 2020; pp. 43–96.
  18. Welch, S.M.; Croft, B.A.; Brunner, J.F.; Michels, M.; Welch, B.; Croft, J.; Brunner, M.; Michels, S.M. PETE: An extension phenology modeling system for management of multi-species pest complex. Environ. Entomol. 1978, 7, 482–494.
  19. Nietschke, B.S.; Magarey, R.D.; Borchert, D.M.; Calvin, D.D.; Jones, E. A developmental database to support insect phenology models. Crop Prot. 2007, 26, 1444–1448.
  20. Orlandini, S.; Magarey, R.D.; Woo Park, E.; Sporleder, M.; Kroschel, J. Methods of agroclimatology: Modeling approaches for pests and diseases. In Agroclimatology: Linking Agriculture to Climate, Agronomy Monograph 60; Hatfield, J.H., Sivakuma, M.V.K., Prueger, J.H., Eds.; Wiley: Hoboken, NJ, USA, 2018; pp. 453–488. Available online: https://acsess.onlinelibrary.wiley.com/doi/book/10.2134/agronmonogr60 (accessed on 18 December 2023).
  21. Arnold, C.Y.Y. Maximum-minimum temperatures as a basis for computing heat units. Proc. Soc. Hortic. Sci. 1960, 76, 682–692.
  22. Wang, J.Y. A critique of the heat unit approach to plant response studies. Ecology 1960, 41, 785–790.
  23. Chuine, I.; Régnière, J. Process-based models of phenology for plants and animals. Annu. Rev. Ecol. Evol. Syst. 2017, 48, 159–182.
  24. Mirhosseini, M.A.; Fathipour, Y.; Reddy, G.V.P.R. Arthropod development’s response to temperature: A review and new software for modeling. Ann. Entomol. Soc. Am. 2017, 110, 507–520.
  25. Knight, A.L. Adjusting the phenology model of codling moth (Lepidoptera: Tortricidae) in Washington State apple orchards. Environ. Entomol. 2007, 36, 1485–1493.
  26. Brunner, J.F.; Hoyt, S.C.; Wright, M.A. Codling Moth Control—A New Tool for Timing Sprays; Cooperative Extension Bulletin, 1072; Washington State University: Pullman, WA, USA, 1987.
  27. Barros-Parada, W.; Knight, A.L.; Fuentes-Contreras, E. Modeling codling moth (Lepidoptera: Tortricidae) phenology and predicting egg hatch in apple orchards of the Maule Region, Chile. Chil. J. Agric. Res. 2015, 75, 57–62.
  28. Jorgensen, C.D.; Martinsen, M.E.; Westover, L.J. Validating Michigan State University’s codling moth model (MOTHMDL) in an arid environment. Gt. Lakes Entomol. 1979, 12, 203–212.
  29. Song, Y.H.; Coop, L.B.; Omeg, M.; Rield, H. Development of a phenology model for predicting Western cherry fruit fly, Rhagoletis indifferens Curran (Diptera: Tephritidae), emergence in the mid Columbia area of the western United States. J. Asia. Pac. Entomol. 2003, 6, 187–192.
  30. Régnière, J.; Sharov, A. Phenology of Lymantria dispar (Lepidoptera: Lymantriidae), male flight and the effect of moth dispersal in heterogeneous landscapes. Int. J. Biometeorol. 1998, 41, 161–168.
  31. Régnière, J.; Sharov, A. Simulating temperature-dependent ecological processes at the sub-continental scale: Male gypsy moth flight phenology as an example. Int. J. Biometeorol. 1999, 42, 146–152.
  32. Schaub, L.P.; Ravlin, F.W.; Gray, D.R.; Logan, J.A. Landscape framework to predict phenological events for gypsy moth (Lepidoptera: Lymantriidae) management programs. Environ. Entomol. 1995, 24, 10–18.
  33. Russo, J.M.; Liebhold, A.M.; Kelley, J.G.W. Mesoscale weather data as input to a gypsy moth (Lepidoptera: Lymantriidae) phenology model. J. Econ. Entomol. 1993, 86, 838–844.
  34. Foster, J.R.; Townsend, P.A.; Mladenoff, D.J. Mapping asynchrony between gypsy moth egg-hatch and forest leaf-out: Putting the phenological window hypothesis in a spatial context. For. Ecol. Manag. 2013, 287, 67–76.
  35. Samietz, J.; Graf, B.; Höhn, H.; Schaub, L.P.; Höpli, H.U. SOPRA: Forecasting tool for fruit tree pest insects. Rev. Suisse Vitic. Arboric. Hortic. 2007, 39, 187–194.
  36. Delahaut, K. Insects. In Phenology: An Integrative Environmental Science; Schwartz, M., Ed.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2003; pp. 405–419.
  37. Damos, P.T.; Savopoulou-Soultani, M. Temperature-driven models for insect development and vital thermal requirements. Psyche 2012, 2012, 123405.
  38. Pruess, K.P. Degree-day methods for pest management. Environ. Entomol. 1983, 12, 613–619.
  39. Campbell, A.; Frazer, B.D.; Gilbert, N.; Gutierrez, A.P.; Mackauer, M. Temperature requirements of some aphids and their parasites. J. Appl. Ecol. 1974, 11, 431–438.
  40. Rebaudo, F.; Rabhi, V.B. Modeling temperature-dependent development rate and phenology in insects: Review of major developments, challenges, and future directions. Entomol. Exp. Appl. 2018, 166, 607–617.
  41. Baskerville, G.L.; Emin, P. Rapid estimation of heat accumulation from maximum and minimum temperatures. Ecology 1969, 50, 514–517.
  42. Riedl, H.; Croft, B.A.; Howitt, A.J. Forecasting codling moth phenology based on pheromone trap catches and physiological-time models. Can. Entomol. 1976, 108, 449–460.
  43. Welch, S.M.; Croft, B.A.; Michels, M.F. Validation of pest management models. Environ. Entomol. 1981, 10, 425–432.
  44. Coop, L. What’s New—Online IPM Weather Data, Degree-Day, and Plant Disease Risk Models 2000–2003. Available online: https://uspest.org/wea/weanew03.html (accessed on 22 November 2023).
  45. Coop, L. What’s New—IPM Pest and Plant Disease Models and Forecasting—For Agricultural, Pest Management, and Plant Biosecurity Decision Support in the US. Available online: https://uspest.org/wea/weanew0409.html (accessed on 22 November 2023).
  46. Daly, C. Guidelines for assessing the suitability of spatial climate data sets. Int. J. Climatol. 2006, 26, 707–721.
  47. Willmott, C.J.; Robeson, S.M. Climatologically aided interpolation (CAI) of terrestrial air temperature. Int. J. Climatol. 1995, 15, 221–229.
  48. SAFARIS Brown Marmorated Stink Bug (Halymorpha halys) Phenological Stages. In Spatial Analytic Framework for Advanced Risk Information Systems (SAFARIS); U.S. Dept. of Agriculture (USDA): Washington, DC, USA; North Carolina State University: Raleigh, NC, USA, 2022; Available online: https://safaris.cipm.info (accessed on 23 May 2022).
  49. Crimmins, T.M.; Gerst, K.L.; Huerta, D.G.; Marsh, R.L.; Posthumus, E.E.; Rosemartin, A.H.; Switzer, J.; Weltzin, J.F.; Coop, L.B.; Dietschler, N.; et al. Short-term forecasts of insect phenology inform pest management. Ann. Entomol. Soc. Am. 2020, 113, 139–148.
  50. Régnière, J. Generalized approach to landscape-wide seasonal forecasting with temperature-driven simulation models. Environ. Entomol. 1996, 25, 869–881.
  51. Grevstad, F.S.; Coop, L.B. The consequences of photoperiodism for organisms in new climates. Ecol. Appl. 2015, 25, 1506–1517.
  52. CAPS Resource and Collaboration Site. Animal and Plant Health Inspection Service (APHIS), Plant Protection and Quarantine (PPQ). Available online: https://caps.ceris.purdue.edu (accessed on 18 December 2023).
  53. Jarvis, C.H.; Baker, R.H.A. Risk assessment for nonindigenous pests: 2. Accounting for interyear climate variability. Divers. Distrib. 2001, 7, 237–248.
  54. Jarvis, C.H.; Baker, R.H.A. Risk assessment for nonindigenous pests: I. Mapping the outputs of phenology models to assess the likelihood of establishment. Divers. Distrib. 2001, 7, 223–235.
  55. Sporleder, M.; Juarez, H.; Simon, R.; Kroschel, J. ILCYM-Insect life cycle modeling: Software for developing temperature-based insect phenology models with applications for regional and global pest risk assessments and mapping. In Proceedings of the 15th Triennial ISTRC Symposium of the International Society for Tropical Root Crops (ISTRC), Lima, Peru, 2–6 November 2009; pp. 216–223.
  56. Stoeckli, S.C.; Hirschi, M.; Spirig, C.; Calanca, P.; Rotach, M.W.; Samietz, J. Impact of climate change on voltinism and prospective diapause induction of a global pest insect—Cydia pomonella (L.). PLoS ONE 2012, 7, e35723.
  57. Jakoby, O.; Lischke, H.; Wermelinger, B. Climate change alters elevational phenology patterns of the European spruce bark beetle (Ips typographus). Glob. Chang. Biol. 2019, 25, 4048–4063.
  58. Cornes, R.C.; van der Schrier, G.; van den Besselaar, E.; Jones, P.D. An ensemble version of the E-OBS temperature and precipitation data sets. J. Geophys. Res. 2018, 123, 9391–9409.
  59. Thornton, M.M.; Strestha, R.; Wei, Y.; Thornton, P.E.; Kao, S.-C.; Wilson, B.E. Daymet: Daily Surface Weather Data on a 1-Km Grid for North America, Version 4; Oak Ridge National Laboratory Distributed Active Archive Center: Oak Ridge, TN, USA, 2020.
  60. Thornton, P.E.; Shrestha, R.; Thornton, M.; Kao, S.C.; Wei, Y.; Wilson, B.E. Gridded daily weather data for North America with comprehensive uncertainty quantification. Sci. Data 2021, 8, 190.
  61. Xavier, A.C.; Scanlon, B.R.; King, C.W.; Alves, A.I. New and improved Brazilian daily weather gridded data (1961–2020). Int. J. Climatol. 2022, 42, 8390–8404.
  62. Xavier, A.C.; King, C.W.; Scanlon, B.R. Daily gridded meteorological variables in Brazil (1980–2013). Int. J. Climatol. 2016, 36, 2644–2659.
  63. Fang, S.; Mao, K.; Xia, X.; Wang, P.; Shi, J.; Bateni, S.M.; Xu, T.; Cao, M.; Heggy, E.; Qin, Z. Dataset of daily near-surface air temperature in China from 1979 to 2018. Earth Syst. Sci. Data 2022, 14, 1413–1432.
  64. Qin, R.; Zhao, Z.; Xu, J.; Ye, J.S.; Li, F.M.; Zhang, F. HRLT: A high-resolution (1 d, 1 km) and long-term (1961–2019) gridded dataset for surface temperature and precipitation across China. Earth Syst. Sci. Data 2022, 14, 4793–4810.
  65. Nengzouzam, G.; Hodam, S.; Bandyopadhyay, A.; Bhadra, A. Spatial and temporal trends in high resolution gridded temperature data over India. Asia-Pac. J. Atmos. Sci. 2019, 55, 761–772.
  66. Ali, S.; Bhutta, Z.A.; Reboita, M.S.; Goheer, M.A.; Ebrahimi, S.; Rozante, J.R.; Kiani, R.S.; Muhammad, S.; Khan, F.; Rahman, M.M.; et al. A 5-km gridded product development of daily temperature and precipitation for Bangladesh, Nepal, and Pakistan from 1981 to 2016. Geosci. Data J. 2023; in press.
  67. Jarvis, C.H. GEO_BUG: A geographical modelling environment for assessing the likelihood of pest development. Environ. Model. Softw. 2001, 16, 753–765.
  68. Jarvis, C.H.; Stuart, N. Accounting for error when modelling with time series data: Estimating the development of crop pests throughout the year. Trans. GIS 2001, 5, 327–343.
  69. Jarvis, C.H.; Baker, R.H.A.; Morgan, D. The impact of interpolated daily temperature data on landscape-wide predictions of invertebrate pest phenology. Agric. Ecosyst. Environ. 2003, 94, 169–181.
  70. Jarvis, C.H.; Collier, R.H. Evaluating an interpolation approach for modelling spatial variability in pest development. Bull. Entomol. Res. 2002, 92, 219–231.
  71. Racca, P.; Zeuner, T.; Jung, J.; Kleinhenz, B. Model validation and use of geographic information systems in crop protection warning service. In Precision Crop Protection: The Challenge and Use of Heterogeneity; Oerke, E.C., Gerhards, R., Menz, G., Sikora, R., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 259–276.
  72. Jarvis, C.H.; Stuart, N. A comparison among strategies for interpolating maximum and minimum daily air temperatures. Part I: The selection of “guiding” topographic and land cover variables. J. Appl. Meteorol. 2001, 40, 1060–1074.
  73. Primack, R.B.; Gallinat, A.S.; Ellwood, E.R.; Crimmins, T.M.; Schwartz, M.D.; Staudinger, M.D.; Miller-Rushing, A.J. Ten best practices for effective phenological research. Int. J. Biometeorol. 2023, 67, 1509–1522.
  74. Probert, A.F.; Wegmann, D.; Volery, L.; Adriaens, T.; Bakiu, R.; Bertolino, S.; Essl, F.; Gervasini, E.; Groom, Q.; Latombe, G.; et al. Identifying, reducing, and communicating uncertainty in community science: A focus on alien species. Biol. Invasions 2022, 24, 3395–3421.
  75. Nijhout, H.F. Development and evolution of adaptive polyphenisms. Evol. Dev. 2003, 5, 9–18.
  76. Neslon, R.J.; Denlinger, D.L.; Somers, D.E. Photoperiodism. The Biological Calendar; Oxford University Press: New York, NY, USA, 2010.
  77. Barker, B.S.; Coop, L.; Duan, J.J.; Petrice, T.R. An integrative phenology and climatic suitability model for emerald ash borer. Front. Insect Sci. 2023, 3, 1239173.
  78. Shaffer, P.L. Prediction of variation in development period of insects and mites reared at constant temperature. Environ. Entomol. 1983, 12, 1012–1019.
  79. Wagner, T.L.; Wu, H.-I.; Feldman, R.M.; Sharpe, P.J.H.; Coulson, R.N. Multiple-cohort approach for simulating development of insect populations under variable temperatures. Ann. Entomol. Soc. Am. 1985, 78, 691–704.
  80. Sharpe, P.J.H.; Curry, G.L.; DeMichele, D.W.; Cole, C.L. Distribution model of organism development times. J. Theor. Biol. 1977, 66, 21–38.
  81. Yonow, T.; Zalucki, M.P.; Sutherst, R.W.; Dominiak, B.C.; Maywald, G.F.; Maelzer, D.A.; Kriticos, D.J. Modelling the population dynamics of the Queensland fruit fly, Bactrocera (Dacus) tryoni: A cohort-based approach incorporating the effects of weather. Ecol. Modell. 2004, 173, 9–30.
  82. Howe, R.W. Temperature effects on embryonic development in insects. Annu. Rev. Entomol. 1967, 10, 15–42.
  83. Curry, G.U.Y.L.; Feldman, R.M.; Smith, C. A stochastic model of a temperature-dependent population. Theor. Popul. Biol. 1978, 13, 197–213.
  84. Jönsson, A.M.; Pulatov, B.; Linderson, M.L.; Hall, K. Modelling as a tool for analysing the temperature-dependent future of the Colorado potato beetle in Europe. Glob. Chang. Biol. 2013, 19, 1043–1055.
  85. Ward, S.F.; Moon, R.D.; Herms, D.A.; Aukema, B.H. Determinants and consequences of plant–insect phenological synchrony for a non-native herbivore on a deciduous conifer: Implications for invasion success. Oecologia 2019, 190, 867–878.
  86. Suppo, C.; Bras, A.; Robinet, C. A temperature- and photoperiod-driven model reveals complex temporal population dynamics of the invasive box tree moth in Europe. Ecol. Modell. 2020, 432, 109229.
  87. Rowley, C.; Cherrill, A.; Leather, S.R.; Pope, T.W. Degree-day based phenological forecasting model of saddle gall midge (Haplodiplosis marginata) (Diptera: Cecidomyiidae) emergence. Crop Prot. 2017, 102, 154–160.
  88. Grevstad, F.S.; Wepprich, T.; Barker, B.S.; Coop, L.B.; Shaw, R.; Bourchier, R.S. Combining photoperiod and thermal responses to predict phenological mismatch for introduced insects. Ecol. Appl. 2022, 32, e2557.
  89. Ogburn, E.C.; Ohmen, T.M.; Huseth, A.S.; Reisig, D.D.; Kennedy, G.G.; Walgenbach, J.F. Temperature-driven differences in phenology and habitat suitability for brown marmorated stink bug, Halyomorpha halys, in two ecoregions of North Carolina. J. Pest Sci. 2023, 96, 373–387.
  90. Nielsen, A.L.; Chen, S.; Fleischer, S.J. Coupling developmental physiology, photoperiod, and temperature to model phenology and dynamics of an invasive heteropteran, Halyomorpha halys. Front. Physiol. 2016, 7, 165.
  91. Quesada-Moraga, E.; Valverde-García, P.; Garrido-Jurado, I. The effect of temperature and soil moisture on the development of the preimaginal Mediterranean fruit fly (Diptera: Tephritidae). Environ. Entomol. 2012, 41, 966–970.
  92. Ma, G.; Tian, B.-L.; Zhao, F.; Wei, G.-S.; Hoffmann, A.A.; Ma, C.-S. Soil moisture conditions determine phenology and success of larval escape in the peach fruit moth, Carposina sasakii (Lepidoptera, Carposinidae): Implications for predicting drought effects on a diapausing insect. Appl. Soil Ecol. 2017, 110, 65–72.
  93. McDougall, R.N.; Ogburn, E.C.; Walgenbach, J.F.; Nielsen, A.L. Diapause termination in invasive populations of the brown marmorated stink bug (Hemiptera: Pentatomidae) in response to photoperiod. Environ. Entomol. 2021, 50, 1400–1406.
  94. Bean, D.W.; Dalin, P.; Dudley, T.L. Evolution of critical day length for diapause induction enables range expansion of Diorhabda carinulata, a biological control agent against tamarisk (Tamarix spp.). Evol. Appl. 2012, 5, 511–523.
  95. Efron, B.; Gong, G. A leisurely look at the bootstrap, the jackknife, and cross-validation. Am. Stat. 1983, 37, 36–48.
  96. Hevesi, J.A.; Istok, J.D.; Flint, A.L. Precipitation estimation in mountainous terrain using multivariate geostatistics. Part I: Structural analysis. J. Appl. Meteorol. Climatol. 1992, 31, 661–676.
More
© Text is available under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/)
Video Production Service
Feedback