Next Article in Journal
An Optimized SIFT-OCT Algorithm for Stitching Aerial Images of a Loblolly Pine Plantation
Previous Article in Journal
Comparison of Multiple Machine Learning Models for Estimating the Forest Growing Stock in Large-Scale Forests Using Multi-Source Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Statistics for Analyzing Climatic Conditions Influencing Lymantria dispar Outbreaks

1
Department of Biosystems Machinery Engineering, Chungnam National University, Daejeon 34134, Korea
2
Division of Forest Diseases and Insect Pests, National Institute of Forest Science, Seoul 02455, Korea
3
Department of Applied Biology, Chungnam National University, Daejeon 34134, Korea
4
Department of Smart Agriculture Systems, Chungnam National University, Daejeon 34134, Korea
*
Author to whom correspondence should be addressed.
Forests 2022, 13(9), 1474; https://doi.org/10.3390/f13091474
Submission received: 20 August 2022 / Revised: 6 September 2022 / Accepted: 10 September 2022 / Published: 13 September 2022
(This article belongs to the Section Forest Biodiversity)

Abstract

:
The spongy moth (Lymantria dispar) is a forest pest that damages a variety of trees in North America and Eurasia. A spongy moth outbreak occurred in part of South Korea in 2020 and caused severe damage to domestic forests and human society. Since the occurrence of spongy moths is influenced by climatic factors, this study examines the causes of spongy moth outbreaks by analyzing the temporal and spatial differences in climatic factors, influencing spongy moth occurrence using specimens collected during field surveys. Climatic factors were identified using global occurrence coordinates to compare the weather characteristics of spongy moth occurrence in domestic regions, using the kernel density function. Spatial and temporal comparisons were performed for monthly weather factors obtained from field surveys in 2020 and 2021 in areas with high and low spongy moth larvae densities. Spongy moth outbreaks may result from particular combinations of variable seasonality in temperature and precipitation, including high temperatures during cold periods and low precipitation during developmental periods.

1. Introduction

The spongy moth (Lymantria dispar) is an agricultural and forest pest that has caused damage to various crops and trees, as well as health impacts from mass occurrence [1,2]. It is generally found in temperate areas from northern Africa to Eurasia, reaching the Japanese archipelago, and has now been introduced to France and North America [3,4,5]. In South Korea, spongy moths have been reported throughout the country, except in the southwest regions, causing health concerns in urban residences and forest damage [2,6], thus necessitating the prediction of possible conditions for preventing outbreaks.
Climatic factors have a considerable influence on the distribution and population of insects [7]. Spongy moth habitats are predominantly distributed in temperate regions, suggesting that regional climate determines their distribution range [1]. The effect of climatic conditions on spongy moths is likely a dominant factor influencing their occurrence and population growth [8]. Excessively high or low temperatures which vary by season may increase the mortality of spongy moth eggs, with high precipitation potentially influencing the infection rate for the deadly virus in spongy moths [9,10,11,12,13]. Optimal growth and the survival rate of spongy moth larvae are likely related to temperature, suggesting that spongy moth occurrence is strongly dependent on climatic conditions [14,15].
In 2020, there was a large outbreak of spongy moths in South Korea, causing damage to a wide range of forests and social issues of discomfort due to a large quantity of eggs in residential areas [16,17]. In contrast, spongy moths were not frequently recorded in South Korea in 2019 and 2021, suggesting temporal difference of environments affects the domestic outbreak. In addition, spatial difference of the spongy moth population even in 2020 outbreak indicates the necessity for examining environmental variability determined as the Moran effect [8]. Therefore, this study examines the climatic characteristics of areas where spongy moths are distributed and identifies the weather factors affecting spongy moth outbreaks through spatial and temporal comparisons based on field surveys in areas with large populations of spongy moths. Findings from this research are expected to provide a baseline to facilitate forecasting future spongy moth outbreaks.

2. Materials and Methods

2.1. Distribution of the Spongy Moth to Extract Climate Data

The global occurrence coordinates for spongy moths were obtained from the Global Biodiversity Information Facility (GBIF) and cross-checked with the Center for Agriculture and Bioscience International (CABI) to confirm reliable occurrence points [18,19]. To remove sampling bias and autocorrelation, spatial filtering was performed using the SDM toolbox in ArcMap (version 10.4.1; ESRI, Redlands, CA, USA), resulting in 3934 occurrence coordinates for the spongy moth [20,21] (Figure 1a). The meteorological data were extracted from the final occurrence coordinates.

2.2. Global Meteorological Data from Areas of Spongy Moth Occurrence

Global climate data including the monthly average temperature, monthly minimum temperature, monthly maximum temperature, and monthly precipitation from 1970 to 2000 were obtained from WorldClim with 1 km resolution (www.worldclim.org (accessed on 6 December 2021)) [22]. As the world climate data were provided as an image file, we first converted it to raster format and then extracted the values for four meteorological data points at the coordinates of occurrence using ArcMap.

2.3. Domestic Weather Data

Weather data at the national scale, including the monthly average temperature, monthly minimum temperature, monthly maximum temperature, and monthly total precipitation, were obtained from the Korea Meteorological Administration (KMA) for January to December 2019–2021 (Figure 1). Given that the spongy moth outbreak occurred in 2020 but it was not observed nationwide in 2019 or 2021, weather data were obtained from 2019 to 2021. The weather data were extracted for six areas with high spongy moth larvae density, or where spongy moths were not observed, to compare the spatial and temporal differences in the weather data.

2.4. Domestic Occurrence Data and Field Survey

The data recording areas with high and low density of spongy moth in 2020 was obtained from the National Institute of Forest Science. This data did not include quantitative population density, but indicated the level of density. A field survey was conducted to investigate the density of gypsy moths from March to September 2021 at six sites with high or low density of the spongy moth (Figure 1b). Field surveys were performed once or twice per month at each site, with burlap traps installed on approximately 20 trees. Species of trapped trees, the host plants of spongy moths, were Acer pseudosieboldianum, Ginkgo biloba, Malus floribunda, Cornus officinalis, Juglans regia, Quercus mongolica, Magnolia kobus, Diospyros kaki and Zelkova serrata [23,24,25]. Spongy moths were observed in egg masses in March and September, in the larval state from April to June, and in the pupa and adult state from June to August. The samples collected were stored in absolute ethanol at Chungnam National University, Korea. These field surveys were conducted from May to October 2021, using the same method used for the occurrence regions. A total of 1768 eggs and 1133 larvae were found in March and April, respectively, but all disappeared in June and no adults were found at all. In addition, a total of 20 adults were only found in three low-density areas during the whole field survey.

2.5. Statistical Analysis for Global Climatic and Domestic Weather Factors

To analyze the climatic characteristics for areas where spongy moths have been recorded, the kernel density function was used for four types of global climatic data, as well as data from South Korea. By comparing the similarity of the distribution patterns for the meteorological data, climatic characteristics suitable for spongy moth occurrence were spatially examined. For the domestic meteorological data, ANOVA and post-hoc tests (Tukey tests) were used to compare meteorological data for time and locations with high and low density. We examined whether there were significant differences in weather factors for time and regions according to the larvae density of spongy moths. The statistical analyses were conducted using R version 4.1.2 after grouping the monthly weather factors for each of the three regions classified by density with a significance level of 0.05 [26].

3. Results

3.1. Characterizing Climatic Conditions Suitable for Spongy Moth Occurrence

The distribution of climatic factors at global occurrence coordinates for spongy moths was analyzed using the kernel density function and temporally compared to areas where the domestic spongy moth outbreak was surveyed (Figure 2). Globally, 90% or more areas showed an average monthly temperature, maximum monthly temperature, minimum monthly temperature, and monthly precipitation of −7–21 °C, 2–27 °C, −11–15 °C, and 30–115 mm, respectively. The maximum frequencies were observed at 17.5 °C, 22.8 °C, 12.2 °C, and 63 mm for average monthly temperature, maximum monthly temperature, minimum monthly temperature, and monthly precipitation, respectively. The range of favorable climatic factors was found to be an average monthly temperature of 10–20°C, with an optimal maximum range of 0–25 °C.
For weather factors extracted from the domestic field survey areas from 2020, when the spongy moth outbreak was reported, the temperature ranges were estimated to be from –30 °C with the highest frequency at 6.5 °C, 25.2 °C, −0.3 °C, and 34 mm for the average monthly temperature, maximum monthly temperature, minimum monthly temperature, and monthly precipitation, respectively. In 2021, when there was no spongy moth outbreak in South Korea, there was an average temperature range of −3–25 °C with a maximum frequency of 15.8 °C. The average maximum and minimum temperatures were 3–29 °C, and −8–21 °C, respectively, while the monthly precipitation range was 6–237 mm with maximum frequencies of 41.8 mm.
When comparing the domestic weather characteristics of the area where the spongy moth occurred in 2020 and the global distribution areas, the global distribution was lower than in South Korea in 2020, suggesting that the South Korean temperatures were higher than the conditions to which the species has been adapted. The precipitation distribution encompassed a wider range with high variation, which suggests high precipitation seasonality in 2020, when spongy moth occurrence was high. The global and domestic climate distribution overlap over a wide range and variation patterns are similar, suggesting the potential for the use of climatic preference as a factor in determining the distribution of spongy moths. The distribution of weather factors in 2021, when there was no outbreak, was similar to that of 2020, except for precipitation.

3.2. Spatial Analysis of Domestic Weather Conditions

The kernel density function is useful for identifying favorable weather conditions, but a similar distribution range suggests that it is necessary to analyze outbreak factors using more detailed weather data. In this study, detailed weather factors were recorded monthly during one year for the field survey areas with and without spongy moth outbreaks, which were statistically compared to consider spatial differences. As a result, all four weather factors were found to be significantly different during the winter season, while precipitation and maximum temperature were or were not different depending on the month (Table 1). Statistically significant differences were found for all weather factors in December, and each temperature factor showed a difference of 4–5 °C and precipitation of approximately 6 mm. Temperature factors for November, January, and February differed by location, suggesting that winter temperatures may have affected the spongy moth outbreak. Significant differences were observed for precipitation in February, March, July, September, and December. March precipitation decreased by approximately 40 mm compared with February in high-density areas, which was not reported in low-density areas. There was a significant difference despite the small difference in average precipitation in March compared with other months, suggesting that there was little variation.

3.3. Temporal Analysis of Domestic Weather Conditions

A temporal comparison was conducted for weather factors in 2020, when there were many spongy moths, and in 2019, and 2021, when there were few spongy moths (Table 2). A highly significant difference was observed for precipitation, which differed during the winter and spring seasons, in addition to during August and October. Precipitation was substantially higher in January, February, and August, but was lower from March to May and October. A significant difference can be seen in the precipitation in January and February, which is the winter season. However, since the spongy moth remains as eggs during this season, the effect of precipitation might be less than that of temperature. For temperature, a significant difference was only found during the winter season by approximately 2–6 °C according to month, and the most substantial difference was observed for the minimum temperature in January. In January, all four weather factors were different, with temperatures 2–6 °C lower in 2020 than in 2019 and 2021, with a higher precipitation of approximately 50 mm. As suggested by findings from the spatial comparison, the temperature during the winter season may affect spongy moth outbreaks.

4. Discussion

In previous studies, the times of spongy moth outbreaks have been studied in relation to various environmental factors, of which climate has been shown to be dominant [27]. Consequently, spatiotemporal analyses of field surveys and weather factors between 2019 and 2021 were conducted in this study to consider weather factors preceding the outbreak in 2020 that were different from those preceding the non-outbreak years of 2019 and 2021. Global and domestic climate factors were not found to significantly favor the outbreak in 2020. However, the domestic weather at the occurrence sites had slightly higher temperature ranges and variable precipitation by season. Compared with the weather at the low-density sites, the weather at the occurrence sites varied with more peaks in the kernel density function, suggesting that seasonal variation in weather may be a factor causing spongy moth outbreaks [4]. The difference between the global occurrence data, including both European and Asian spongy moths, and domestic field surveys collecting Asian spongy moths may have resulted in slight differences in climatic data due to their own biological characteristics [28]. In addition, outbreaks typically progress through several phases, including population growth, peak, and declining years, suggesting the lasting effect of climate. For this reason, it is important to analyze climatic factors with consideration to seasonal variation and spatial presence or absence of a species to determine the reason for the outbreak.
A spatial comparative analysis was performed by selecting an area according to the density of larvae at the time of the spongy moth outbreak in 2020. This was done to confirm the role of weather factors by revealing the effect of weather difference by location, because simple analysis on weather over time in the same site could not explain the spatial difference in 2020. Therefore, we tried to confirm the role of weather factors by showing that climate differences due to spatial differences influence the occurrence of spongy moths. For this analysis, it was possible to identify different weather factors by month, suggesting that the winter weather affected spongy moth outbreaks. The eggs of spongy moths begin to hatch when exposed to a temperature of approximately 10 °C after being exposed to low temperatures [29], suggesting the appropriate hatching timing (March) may be slightly advanced because of the relatively high temperature in February. This means that spongy moth eggs in low-density areas hatch earlier than those in high-density areas. As the larvae hatched earlier in February, when the average monthly temperature was 4.6 °C because of early hatching, the minimum growth temperature of 7 °C was not reached, which likely influenced the larval mortality rate [15]. Monthly precipitation appears to affect the spatial differences in spongy moth outbreaks. Precipitation is a factor that substantially affects the spread of the fungal pathogen Entomophaga maimaiga and the LdNPV virus, which is fatal to spongy moths, and high precipitation limits the population of spongy moths [10,30]. Therefore, high winter temperature and precipitation when spongy moths begin to develop after hatching may result in a low density of spongy moths. Warm winter temperature and adequate precipitation at the time of development may be factors leading to gypsy moth outbreaks. In particular, the high winter temperature, enough to cause early hatching due to the high temperature in winter, was achieved in 2020 in contrast to the normal winter temperature in 2021, causing outbreak of the spongy moth. There might be other factors affecting the spongy moth population, such as topology, host plants, and local biodiversity [31], but they are relatively constant compared with weather factors. In addition, all six sites of the field survey were near mountainous areas, suggesting that surrounding environment was similar. Therefore, a spatial comparison of weather factors could find a potential cause of the spongy moth outbreak.
The temporal comparison confirmed the effect of winter weather on spongy moth outbreaks. The daily minimum temperature in January, the coldest month in South Korea, was approximately −3–6 °C in 2020, while it was −8–4 °C in 2019, and −9–3 °C in 2021. This low winter temperature may have increased the mortality rate with the exposure of spongy moth eggs to extreme weather. If the winter temperature reaches −29 °C, the supercooling point of the spongy moth egg, or if there is a period in which the temperature is maintained at approximately −15 °C to −20 °C, spongy moth outbreaks may not occur [13]. When considering the high winter temperatures recorded in areas with low density in 2020, it is likely that excessively high or low winter temperatures can inhibit the occurrence of spongy moths. As in the case of spatial analysis, precipitation was significantly different during spring, suggesting high levels of precipitation around the main period of development for the spongy moth larvae after hatching limited spongy moth outbreaks. Moderate winter temperatures and low spring precipitation may cause outbreaks of spongy moths. Consequently, it may be possible to forecast outbreaks of spongy moths in spring based on the previous winter weather. This is helpful for establishing a strategy for monitoring and controlling the spongy moth, which can subsequently reduce losses in forests and agriculture, including the fruit industry that suffers due to spongy moth outbreaks.
Although weather factors were analyzed as major factors influencing the occurrence of spongy moths in South Korea, the population density of the Lepidoptera family, which has a cyclic population, is known to be affected by parasitoids [32]. Previous studies have suggested that weather factors such as temperature and precipitation were related to the population of parasitoids [33,34], suggesting the role of weather and climate factors on the spongy moth population. In contrast, another study showed less correlation between parasitoids population and weather factors [35]. Since Lepidopteran populations are often cyclic, usually with long time intervals of outbreak [32], long-term observation may be required to accurately identify the outbreak of the gypsy moth. Compared to the species distribution modeling method that directly evaluates potential distribution of a species, this climate analysis method cannot directly predict the potential distribution area. However, it can easily identify the climatic characteristics of the habitat area without a complex modeling process, and can identify noticeable climatic differences through quantitative comparisons. Therefore, a more reliable result can be produced by comparing the climatic characteristics of the potentially habitable area through the prediction by species distribution modeling and data based on distribution information [36,37].

5. Conclusions

In this study, we aimed to determine the cause of spongy moth outbreaks in South Korea by analyzing the climatic conditions. There are some drawbacks in this study, including the absence of quantified population density applicable for investigating the correlation between meteorological factors and population, a limited number of sites, and other environmental factors that affect spongy moth development. However, this study applied simple spatiotemporal statistics which did not require quantitative population data to explore the role of climatic factors, and proposed potential reasons for the outbreak, which were winter weather and spring precipitation, along with weather seasonality. According to data released by the Korea Meteorological Administration, winter temperatures in 2020 were abnormally high throughout the country, with an average of 3.1 °C, the warmest winter since 1973. The increase in mortality due to the cold winter weather did not occur in areas with spongy moth outbreaks, whereas in the southern regions, the early hatching time accelerated due to the excessively hot winter temperatures likely inhibiting spongy moth occurrence due to insufficient food supplies. High precipitation during the spring season increases fungal pathogens, limiting the development of spongy moths. Therefore, outbreaks of spongy moths can be forecasted through the analysis of climate factors in winter and spring, which is helpful for establishing a control strategy.

Author Contributions

Conceptualization, W.-H.L. and S.-H.J., methodology, J.-W.S., S.-H.Y. and W.-H.L.; software, J.-W.S., J.-M.J. and S.-H.Y.; validation, J.-W.S. and J.-M.J.; formal analysis, J.-W.S. and H.-M.B.; investigation, J.-W.S., H.-M.B. and S.-H.J.; resources, Y.-W.N.; data curation, J.-W.S. and J.-M.J.; writing—original draft preparation, J.-W.S. and W.-H.L.; writing—review and editing, J.-W.S., S.-H.J. and W.-H.L.; visualization, J.-W.S.; supervision, S.-H.J. and W.-H.L.; project administration, Y.-W.N. and S.-H.J.; funding acquisition, S.-H.J. and W.-H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the “Investigation of larval occurrence characteristics of gypsy moths (Project no. FE0702-2021-01-2021)” funded by the National Institute of Forest Science in Korea.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to institutional policy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Montgomery, M.E.; Wallner, W.E. The gypsy moth. A westward migrant. In Dynamics of Forest Insect Populations: Patterns, Causes, Implications; Berryman, A.A., Ed.; Plenum Press: New York, NY, USA, 1988; pp. 353–376. ISBN 978-1-4899-0789-9. [Google Scholar]
  2. APHIS, U. Pest Alert: Asian Gypsy Moth. Agriculture USDo, Ed. 2016. Available online: https://www.aphis.usda.gov/publications/plant_health/content/printable_version/fs_phasiangm.pdf (accessed on 9 March 2022).
  3. Leonard, D.E. Recent developments in ecology and control of the gypsy moth. Annu. Rev. Entomol. 1974, 19, 197–229. [Google Scholar] [CrossRef]
  4. Elkinton, J.S.; Liebhold, A.M. Population dynamics of gypsy moth in North America. Ann. Rev. Entomol. 1990, 35, 571–596. [Google Scholar] [CrossRef]
  5. Keena, M.A.; Cote, M.J.; Grinberg, P.S.; Wallner, W.E. World distribution of female flight and genetic variation in Lymantria dispar (Lepidoptera: Lymantriidae). Environ. Entomol. 2014, 37, 636–649. [Google Scholar] [CrossRef]
  6. Tuthill, R.W.; Canada, A.T.; Wilcock, K.; Etkind, P.H.; O’Dell, T.M.; Shama, S.K. An epidemiologic study of gypsy moth rash. Am. J. Public Health 1984, 74, 799–803. [Google Scholar] [CrossRef] [PubMed]
  7. Halsch, C.A.; Shapiro, A.M.; Fordyce, J.A.; Nice, C.C.; Thorne, J.H.; Waetjen, D.P.; Forister, M.L. Insects and Recent Climate Change. Proc. Natl. Acad. Sci. USA 2021, 118, e2002543117. [Google Scholar] [CrossRef]
  8. Williams, D.W.; Liebhold, A.M. Influence of weather on the synchrony of gypsy moth (Lepidoptera: Lymantriidae) outbreaks in New England. Environ. Entomol. 1995, 24, 987–995. [Google Scholar] [CrossRef]
  9. Summers, J.N. Effect of Low Temperature on the Hatching of Gipsy-Moth Eggs; US Department of Agriculture: Washington, DC, USA, 1922. [CrossRef]
  10. D’Amico, V.; Elkinton, J.S. Rainfall Effects on Transmission of Gypsy Moth (Lepidoptera: Lymantriidae) Nuclear Polyhedrosis Virus. Environ. Entomol. 1995, 24, 1144–1149. [Google Scholar] [CrossRef]
  11. Thompson, L.M.; Faske, T.M.; Banahene, N.; Grim, D.; Agosta, S.J.; Parry, D.; Tobin, P.C.; Johnson, D.M.; Grayson, K.L. Variation in growth and developmental responses to supraoptimal temperatures near latitudinal range limits of gypsy moth Lymantria dispar (L.), an expanding invasive species. Physiol. Entomol. 2017, 42, 181–190. [Google Scholar] [CrossRef]
  12. Fält-Nardmann, J.J.J.; Ruohomäki, K.; Tikkanen, O.P.; Neuvonen, S. Cold hardiness of Lymantria monacha and L. dispar (Lepidoptera: Erebidae) eggs to extreme winter temperatures: Implications for predicting climate change impacts. Ecol. Entomol. 2018, 43, 422–430. [Google Scholar] [CrossRef]
  13. Ananko, G.G.; Kolosov, A.V. Asian gypsy moth (Lymantria dispar L.) populations: Tolerance of eggs to extreme winter temperatures. J. Therm. Biol. 2021, 102, 103123. [Google Scholar] [CrossRef]
  14. Logan, J.A.; Casagrande, P.A.; Liebhold, A.M. Modeling environment for simulation of gypsy moth (Lepidoptera: Lymantriidae) larval phenology. Environ. Entomol. 1991, 20, 1516–1525. [Google Scholar] [CrossRef]
  15. Limbu, S.; Keena, M.; Chen, F.; Cook, G.; Nadel, H.; Hoover, K. Effects of temperature on development of Lymantria dispar asiatica and Lymantria dispar japonica (Lepidoptera: Erebidae). Environ. Entomol. 2017, 46, 1012–1023. [Google Scholar] [CrossRef] [PubMed]
  16. Jung, J.K.; Nam, Y.W.; Kim, D.S.; Lee, S.H.; Lim, J.H.; Choi, W.I.; Kim, E.S. Tree-crown defoliation caused by outbreak of forest insect pests in Korea during 2020. Korean J. Appl. Entomol. 2020, 59, 409–410. [Google Scholar] [CrossRef]
  17. Choi, W.I.; Kim, E.S.; Yun, S.J.; Lim, J.H.; Kim, Y.E. Quantification of one-year gypsy moth defoliation extent in Wonju, Korea, Using Landsat Satellite Images. Forests 2021, 12, 545. [Google Scholar] [CrossRef]
  18. GBIF. Available online: https://doi.org/10.15468/dl.3mny9k (accessed on 28 April 2021).
  19. CABI. Available online: https://www.cabi.org/isc/datasheet/31807 (accessed on 3 May 2021).
  20. Boria, R.A.; Olson, L.E.; Goodman, S.M.; Anderson, R.P. Spatial filtering to reduce sampling bias can improve the performance of ecological niche models. Ecol. Model. 2014, 275, 73–77. [Google Scholar] [CrossRef]
  21. Brown, J.L.; Anderson, B. SDMtoolbox: A python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods Ecol. Evol. 2014, 5, 694–700. [Google Scholar] [CrossRef]
  22. Fick, S.E.; Hijmans, R.H. WorldClim 2: New 1km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  23. Schaefer, P.W. Betula platyphylla: The preferred oviposition host of Lymantria dispar japonica in Hokkaido, Japans. Environ. Entomol. 1974, 7, 168–170. [Google Scholar] [CrossRef]
  24. Liebhold, A.M.; Gottschalk, K.W.; Muzika, R.M.; Montgomery, M.E.; Young, R.; O’Day, K.; Kelley, B. Suitability of North American Tree Species to the Gypsy Moth: A Summary of Field and Laboratory Tests; General Technical Report NE-211; USDA Forest Service: Radnor, PA, USA, 1995. [CrossRef]
  25. Pogue, M.; Schaefer, P.W. A Review of Selected Species of Lymantria Hübner (1819) (Lepidoptera: Noctuidae: Lymantriinae) from Subtropical and Temperate Regions of Asia, Including the Descriptions of Three New Species, Some Potentially Invasive to North America; USDA, Forest Health Technology Enterprise Team: Washington, DC, USA, 2007.
  26. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 17 December 2021).
  27. Liebhold, A.; Elkinton, J.; Williams, D.; Muzika, R.M. What causes outbreaks of the gypsy moth in North America? Popul. Ecol. 2000, 42, 257–266. [Google Scholar] [CrossRef]
  28. Keena, M.A. Inheritance and world variation in thermal requirements for egg hatch in Lymantria dispar (Lepidoptera: Erebidae). Environ. Entomol. 2016, 45, 1–10. [Google Scholar] [CrossRef]
  29. Wei, J.; Luo, Y.Q.; Shi, J.; Wang, D.P.; Shen, S.W. Impact of temperature on postdiapause and diapause of the Asian gypsy moth, Lymantria dispar asiatica. J. Insect Sci. 2014, 14, 5. [Google Scholar] [CrossRef] [PubMed]
  30. Hajek, A.E.; Elkinton, J.S.; Witcosky, J.J. Introduction and spread of the fungal pathogen Entomophaga maimaiga (Zygomycetes: Entomophthorales) along the leading edge of gypsy moth (Lepidoptera: Lymantriidae) spread. Environ. Entomol. 1996, 25, 1235–1247. [Google Scholar] [CrossRef]
  31. Alalouni, U.; Schadler, M.; Brandl, R. Natural enemies and environmental factors affecting the population dynamics of the gypsy moth. J. Appl. Entomol. 2013, 137, 721–738. [Google Scholar] [CrossRef]
  32. Berryman, A.A. What causes population cycles of forest Lepidoptera? Trends Ecol. Evol. 1996, 11, 28–32. [Google Scholar] [CrossRef]
  33. Murúa, G.; Molina-Ochoa, J.; Coviella, C. Population dynamics of the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae) and its parasitoids in northwestern Argentina. Fla. Entomol. 2006, 89, 175–182. [Google Scholar] [CrossRef]
  34. Stanley, J.; Chandrasekaran, S.; Preetha, G. Conogethes punctiferalis (Lepidoptera: Pyralidae) its biology and field parasitization. Indian J. Agric. Sci. 2009, 79, 906–909. [Google Scholar]
  35. Mahadi, N.A.; Muhamad, R.; Adam, N.A. Relationship between bagworm Pteroma pendula Joannis (Lepidoptera: Psychidae) populations, parasitoids, and weather parameters in oil palm plantation. J. Agric. Sci. 2012, 4, 13–17. [Google Scholar] [CrossRef]
  36. Paini, D.R.; Mwebaze, P.; Kuhnert, P.M.; Kriticos, D.J. Global establishment threat from a major forest pest via international shipping: Lymantria dispar. Sci. Rep. 2018, 8, 13723. [Google Scholar] [CrossRef]
  37. Kim, S.H.; Kim, D.E.; Lee, H.; Jung, S.; Lee, W.H. Ensemble evaluation of the potential risk areas of yellow-legged hornet distribution. Environ. Monit. Assess. 2021, 193, 601. [Google Scholar] [CrossRef]
Figure 1. Distribution areas for spongy moths with (a) global occurrence coordinates with actual distribution sites (red dots), and (b) field survey areas with high and low spongy moth density applied for climatic analysis.
Figure 1. Distribution areas for spongy moths with (a) global occurrence coordinates with actual distribution sites (red dots), and (b) field survey areas with high and low spongy moth density applied for climatic analysis.
Forests 13 01474 g001
Figure 2. Kernel density distribution of the major climatic factors in global occurrence coordinates obtained for 30 years (1970–2000), and in domestic field surveys in 2020 and 2021. Dark blue areas include 90% of the total distribution with 5% and 95% of the lower and upper threshold, respectively (x-axis has the unit of temperature (°C), and precipitation (mm), while y-axis shows the kernel density).
Figure 2. Kernel density distribution of the major climatic factors in global occurrence coordinates obtained for 30 years (1970–2000), and in domestic field surveys in 2020 and 2021. Dark blue areas include 90% of the total distribution with 5% and 95% of the lower and upper threshold, respectively (x-axis has the unit of temperature (°C), and precipitation (mm), while y-axis shows the kernel density).
Forests 13 01474 g002
Table 1. Spatial analysis of weather data in areas with high and low spongy moth density in 2020.
Table 1. Spatial analysis of weather data in areas with high and low spongy moth density in 2020.
Weather FactorDensityJan.Feb.Mar.Apr.MayJun.Jul.Aug.Sep.Oct.Nov.Dec.Average
Monthly average temperature (°C)High1.1 a2.1 a6.9 a10 a17.1 a22.5 a22.4 a25.3 a19 a11.9 a6.8 a−2 a11.9 a
Low3.9 b4.6 b8.3 a10.9 a17.8 a22.6 a22.7 a27.2 a20.7 a14.8 a9.8 b2 b13.8 a
Monthly maximum temperature (°C)High6.4 a7.7 a13.7 a16.9 a23 a28.1 a26.4 a29.5 a24.1 a18.8 a12.9 a3.9 a17.6 a
Low8.5 a10.2 b14.8 a17.7 a23.4 a28 a26.5 a31.5 b25.7 a21.2 b15.4 b7.7 b19.2 a
Monthly minimum temperature (°C)High−3.3 a−3.1 a0.2 a2.8 a11.6 a17.4 a19 a22.2 a14.9 a6 a1.3 a−7.1 a6.8 a
Low0.2 b−0.3 a2.1 a4.6 a12.9 a18.1 a20.1 a24 a16.8 a9.5 a5 b−2.6 b9.2 a
Precipitation (mm)High59.7 a60.1 a20 a20.3 a93.7 a95.6 a277.9 a516.9 a156.5 a10.5 a36.6 a4.7 a112.7 a
Low79.1 a41.8 b45.6 b53.1 a126.5 a252.3 a463.9 b448.5 a225.5 b20 a32.3 a11.1 b150 a
Data with different letters mean a significant difference in monthly factors between regions based on the Tukey’s test.
Table 2. Temporal weather data analysis in areas with high spongy moth density in 2019, 2020, and 2021.
Table 2. Temporal weather data analysis in areas with high spongy moth density in 2019, 2020, and 2021.
Weather FactorYearJan.Feb.Mar.Apr.MayJun.Jul.Aug.Sep.Oct.Nov.Dec.Average
Monthly average temperature (°C)2019−2.2 a0.6 a5.9 a10.7 a17.6 a20.7 a24.1 a24.9 a20.4 a14 a6.4 a0.6 a12 a
20201.1 b2.1 a6.9 a10 a17.1 a22.5 a22.4 a25.3 a19 a11.9 a6.8 a−2 b11.9 a
2021−3.4 a1.6 a7.6 a12.2 a15.1 a20.7 a25 a23.2 a19.8 a13.2 a6.6 a0.2 a11.8 a
Monthly maximum temperature (°C)20194.3 a6.7 a12.9 a17.7 a25.4 a26.9 a28.8 a30.3 a25.6 a20.1 a13.5 a6.1 a18.2 a
20206.4 b7.7 a b13.7 a16.9 a23 a b28.1 a26.4 a29.5 a24.1 a18.8 a12.9 a3.9 b17.6 a
20212.6 c8.1 b14 a18.5 a20.9 b26 a30.2 a27.7 a24.7 a19.5 a13.2 a5.9 a17.6 a
Monthly minimum temperature (°C)2019−8.3 a−4.8 a−0.5 a3.9 a9.5 a15 a20.2 a20.5 a16 a8.5 a0.6 a−4.1 a6.4 a
2020−3.3 b−3.1 a0.2 a2.8 a11.6 a17.4 a19 a22.2 a14.9 a6 a1.3 a−7.1 a6.8 a
2021−9.2 a−4.4 a1.5 a5.6 a9.5 a16 a20.7 a19.8 a15.5 a8.4 a1.5 a−4.8 a6.7 a
Precipitation (mm)20192.3 a26.8 a34.5 a49.2 a b23.9 a103.2 a189.6 a116.7 a151.4 a81.3 a80.4 a15.2 a72.9 a
202059.7 b60.1 b20 a20.3 a93.7 b95.6 a277.9 a516.9 b156.5 a10.5 b36.6 a4.7 b112.7 a
202112.8 c8.3 c82 b78.6 b156.7 c82.8 a164.8 a210.4 a198 a45.6 a b55.9 a4.4 b91.7 a
Data with different letters mean a significant difference in monthly factors between regions based on the Tukey’s test.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, J.-W.; Jung, J.-M.; Nam, Y.-W.; Byun, H.-M.; Yoon, S.-H.; Jung, S.-H.; Lee, W.-H. Spatiotemporal Statistics for Analyzing Climatic Conditions Influencing Lymantria dispar Outbreaks. Forests 2022, 13, 1474. https://doi.org/10.3390/f13091474

AMA Style

Song J-W, Jung J-M, Nam Y-W, Byun H-M, Yoon S-H, Jung S-H, Lee W-H. Spatiotemporal Statistics for Analyzing Climatic Conditions Influencing Lymantria dispar Outbreaks. Forests. 2022; 13(9):1474. https://doi.org/10.3390/f13091474

Chicago/Turabian Style

Song, Jae-Woo, Jae-Min Jung, Young-Woo Nam, Hye-Min Byun, Sun-Hee Yoon, Sung-Hoon Jung, and Wang-Hee Lee. 2022. "Spatiotemporal Statistics for Analyzing Climatic Conditions Influencing Lymantria dispar Outbreaks" Forests 13, no. 9: 1474. https://doi.org/10.3390/f13091474

APA Style

Song, J. -W., Jung, J. -M., Nam, Y. -W., Byun, H. -M., Yoon, S. -H., Jung, S. -H., & Lee, W. -H. (2022). Spatiotemporal Statistics for Analyzing Climatic Conditions Influencing Lymantria dispar Outbreaks. Forests, 13(9), 1474. https://doi.org/10.3390/f13091474

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop