Modeling Study of Factors Determining Efficacy of Biological Control of Adventive Weeds

We model the spatiotemporal dynamics of a community consisting of competing weed and cultivated plant species and a population of specialized phytophagous insects used as the weed biocontrol agent. The model is formulated as a PDE system of taxis–diffusion–reaction type and computer-implemented for one-dimensional and two-dimensional cases of spatial habitat for the Neumann zero-flux boundary condition. In order to discretize the original continuous system, we applied the method of lines. The obtained system of ODEs is integrated using the Runge–Kutta method with a variable time step and control of the integration accuracy. The numerical simulations provide insights into the mechanism of formation of solitary population waves (SPWs) of the phytophage, revealing the factors that determine the efficacy of combined application of the phytophagous insect (classical biological method) and cultivated plant (phytocenotic method) to suppress weed foci. In particular, the presented results illustrate the stabilizing action of cultivated plants, which fix the SPW effect by occupying the free area behind the wave front so that the weed remains suppressed in the absence of a phytophage.