Patterns for Competing Populations with Species Specific Nonlocal Coupling

A. Bayliss¹ ; V. A. Volpert¹

¹ Department of Engineering Sciences and Applied Mathematics, Northwestern University, Evanston, IL 60208-3125, USA

Mathematical modelling of natural phenomena, Tome 10 (2015) no. 6, pp. 30-47

Résumé

We study a reaction-diffusion population model for two competing species, u and v, which also includes nonlocal integral terms that represent competition for resources. The integral terms have the form of a convolution of a given kernel function and the solution. They manifest the fact that consumption of resources by the species at a spatial location x depends not just on the populations at point x but rather on weighted averages of the populations in an interval about x. The kernel functions that we employ are characterized by two parameters, δ, which gives a spatial scale of the nonlocality, with δ = 0 corresponding to the local case, and α, a parameter associated with the extent of the asymmetry of the kernel, where α = 0 corresponds to a symmetric kernel, i.e., to a kernel that is an even function of its argument. We consider the parameter regime where for δ = 0 the system admits a stable coexistence equilibrium, which is destabilized for sufficiently large δ. Deeply in the instability region, the ensuing patterns involve arrays of islands, regions of nonzero population, separated by deadzones where the populations are essentially extinct. These structures are stationary if α = 0 and propagate if α ≠ 0. Unlike previous work, we study the patterns when the kernel parameters are different for the two species, focusing first on different δ with α = 0 and then on different α with the same δ. In the first case we numerically find stationary patterns consisting of (i) plateau regions for the more local species (v), where the plateau value is essentially the v −component of the unstable equilibrium point, (ii) oscillatory patterns, where oscillations appear along the plateau and along the u −island and (iii) split patterns, where instead of an island, there is an island cluster of two or more neighbors separated by deadzones or near deadzones. In the second case we consider weakly coupled systems and show that there are two kinds of patterns, (i) bound waves, when each species propagates with the same speed and (ii) unbound patterns where each species propagates with its own speed, interacting with the other species when they intersect, i.e., when the two species occupy the same or closely spaced regions. We show that whether the patterns are bound or unbound depends on the stability of the decoupled system, i.e., the system when the coupling parameters are zero. When one of the species for the decoupled system is stable, then the ensuing patterns for weak coupling is bound, while if both species are unstable for the uncoupled systems, the resulting patterns are unbound for weak coupling.

Détail

DOI : 10.1051/mmnp/201510604

Bibliographie

[1] N. F. Britton J. Theoret. Biol. 1989 57 66

[2] N. F. Britton SIAM. J. Appl. Math. 1990 1663 1688

[3] M. A. Fuentes, M. N. Kuperman, V. M. Kenkre Phys. Rev. Letters 2003 158104 1

[4] M. A. Fuentes, M. N. Kuperman, V. M. Kenkre J. Phys. Chem. B 2004 10505 10508

[5] S. Genieys, N. Bessonov, V. Volpert Mathematical and Computer Modelling 2009 2109 2115

[6] S. Genieys, V. Volpert, P. Auger Math. Model. Nat. Phenom. 2006 63 80

[7] S. Genieys, V. Volpert, P. Auger C. R. Biologies 2006 876 879

[8] S. S. Gourley, M. A. J. Chaplain, F. A. Davidson Dynamical Systems 2001 173 192

[9] S. M. Merchant, W. Nagata Theoretical Population Biology 2011 289 297

[10] B. Perthame, S. Genieys Math. Model. Nat. Phenom. 2007 135 151

[11] B. L. Segal, V. A. Volpert, A. Bayliss Physica D 2013 12 22

[12] M. C. Tanzy, V. A. Volpert, A. Bayliss, M. E. Nehrkorn Mathematical Biosciences 2013 14 26

[13] M. C. Tanzy, V. A. Volpert, A. Bayliss, M. E. Nehrkorn Mathematical Biosciences 2015 70 82