Gypsy moth invasion in North America: A simulation study of the spatial pattern and the rate of spread

Gypsy moth is regarded as one of the top most harmful invasive species. Its invasion in the northeastern US has led to widespread forest defoliation, wildlife disruption and even a change in biogeochemical conditions over the area of 10⁶ km². Spread of gypsy moth has a few distinct features such as a patchy spatial distribution of the gypsy moth population, which is largely uncorrelated to the environmental heterogeneity, and a high variability (almost over an order of magnitude) in the spread rates. These features are usually explained by human-assisted dispersal, e.g. when masses of gypsy moth eggs are inadvertently transported by cars and vehicles. This theory, however, somewhat disagrees with the existence of the strong Allee effect that tends to wipe out small new colonies. In this paper, we suggest an alternative explanation that the patchy structure can result from the interplay between two natural factors such as wind dispersal and viral infection. In order to check this hypothesis, we describe the gypsy moth spread with a diffusive SI model and study its properties by means of extensive computer simulations. Interestingly, in a certain parameter range our model shows formation of spatial patterns that are qualitatively similar to those observed in the field. To find out the relevant parameter range, we make a careful review of available literature sources. For biologically meaningful parameter values, we then show that the rates of gypsy moth spread predicted by our model are in good agreement with the lower band of the rates observed in nature.