Abstract
Reproduction-Dispersal equations, called reaction-diffusion equations in the physics literature, model the growth and spreading of biological species. Integro-Difference equations were introduced to address the shortcomings of this model, since the dispersal of invasive species is often more widespread than what the classical RD model predicts. In this paper, we extend the RD model, replacing the classical second derivative dispersal term by a fractional derivative of order 1<α ≤ 2. Fractional derivative models are used in physics to model anomalous super-diffusion, where a cloud of particles spreads faster than the classical diffusion model predicts. This paper also establishes a connection between the new RD model and a corresponding ID equation with a heavy tail dispersal kernel. The general theory developed here accommodates a wide variety of infinitely divisible dispersal kernels that adapt to any scale. Each one corresponds to a generalised RD model with a different dispersal operator. The connection established here between RD and ID equations can also be exploited to generate convergent numerical solutions of RD equations along with explicit error bounds.