Fisher-KPP Wave Simulator: Population Invasion Fronts

simulator intermediate ~10 min
Loading simulation...
c* = 2.0 km/yr — minimum Fisher-KPP wave speed

With D=1 km²/yr and r=1/yr, the minimum traveling wave speed is c*=2 km/yr. The population front advances steadily, converting empty habitat into fully occupied territory at carrying capacity.

Formula

∂u/∂t = D ∂²u/∂x² + ru(1 − u/K)
c* = 2√(Dr) (minimum wave speed)
u(x,t) → (1 + Ce^(−λ(x−ct)))⁻¹ (traveling wave profile)

Waves of Advance

In 1937, R.A. Fisher asked a fundamental question: how fast does a beneficial gene spread through a spatially distributed population? His answer — the Fisher equation combining diffusion with logistic growth — predicted traveling wavefronts advancing at a speed determined by just two parameters: the diffusion coefficient and the growth rate. The same year, Kolmogorov, Petrovsky, and Piskunov independently proved rigorous mathematical properties of this equation, establishing it as a cornerstone of mathematical biology.

The Minimum Speed Selection

The Fisher-KPP equation admits traveling wave solutions for any speed c ≥ c* = 2√(Dr), but nature consistently selects the minimum speed. This remarkable result arises because the dynamics at the leading edge of the front — where population density is low and growth is approximately exponential — determine how fast the wave advances. The minimum speed represents the balance between diffusive spreading and local growth at the frontier.

Invasion Biology

The Fisher wave provides a quantitative framework for understanding biological invasions. The spread of muskrats across Europe after their introduction in 1905, the advance of cane toads across Australia, and the northward migration of tree species after the last ice age all show approximately constant-speed advancing fronts. By measuring D and r independently, ecologists can predict invasion speeds before they occur — critical for managing invasive species.

Beyond Biology

The Fisher-KPP equation appears far beyond ecology. In epidemiology, plague wavefronts during the Black Death advanced across Europe at speeds consistent with Fisher wave predictions. In oncology, glioblastoma tumor margins advance as Fisher waves through brain tissue. In population genetics, 'gene surfing' at expansion fronts can fix normally rare alleles. The equation even appears in combustion theory and chemical kinetics, wherever diffusion competes with autocatalytic growth.

FAQ

What is the Fisher-KPP equation?

The Fisher-KPP equation (∂u/∂t = D∇²u + ru(1−u)) combines diffusion with logistic growth. Independently proposed by R.A. Fisher (1937) for gene frequency spread and by Kolmogorov, Petrovsky, and Piskunov (1937), it predicts traveling wavefronts that advance at a minimum speed c* = 2√(Dr).

What determines the wave speed?

The minimum (and selected) wave speed is c* = 2√(Dr), depending only on the diffusion coefficient D and the linear growth rate r. Remarkably, the carrying capacity K does not affect the speed — only how diffusion and reproduction compete at the leading edge of the front matters.

Where is the Fisher equation applied?

Applications include biological invasions (spread of muskrats across Europe, cane toads in Australia), epidemic wavefronts (Black Death spreading across medieval Europe at ~2-5 km/day), tumor growth, gene surfing in expanding populations, and the spread of agriculture during the Neolithic revolution.

Why is the minimum speed selected?

For compactly supported initial conditions (population localized in space), the wavefront always approaches the minimum speed c* = 2√(Dr). This is because the leading edge of the population front decays exponentially, and steeper initial profiles select faster waves, but compact initial data select the minimum.

Sources

Embed

<iframe src="https://homo-deus.com/lab/mathematical-biology/fisher-wave/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub