Spatial Predator-Prey Metapopulation

Patch dynamics on a lattice — local Lotka-Volterra with dispersal
∂P/∂t = rP − αPQ + D∇²P  |  ∂Q/∂t = βPQ − dQ + D∇²Q
0.80
0.50
0.30
0.50
0.05
Prey
Predator
Empty
Total Prey
Total Predators
Occupied patches