IMSP schemes for spatially explicit models of cyclic populations and metapopulation dynamics

We examine spatially explicit models described by reaction–diffusion partial differential equations for the study of predator–prey population dynamics. The numerical methods we propose are based on the coupling of a finite difference/element spatial discretization and a suitable partitioned Runge–Kutta scheme for the approximation in time. The RK scheme here implemented uses an implicit scheme for the stiff diffusive term and a partitioned RK symplectic scheme for the reaction term (IMSP schemes). We revisit some results provided in the literature for the classical Lotka–Volterra system and the Rosenzweig–MacArthur model. We then extend the approach to metapopulation dynamics in order to numerically investigate the effect of migration through a corridor connecting two habitat patches. Moreover, we analyze the synchronization properties of subpopulation dynamics, when the migration occurs through corridors of variable sizes.