W. P. Petersen,S. Callegari,N. Tkachenko,J. D. Weissmann,Ch. P. E. Zollikofer
Abstract:We describe a general Godunov type splitting for numerical simulations of the Fisher/Kolmogorov-Petrovski-Piskunov growth and diffusion equation in two spatial dimensions. In particular, the method is appropriate for modeling population growth and dispersal on a terrestrial map. The procedure is semi-implicit, hence quite stable, and approximately second order accurate, excluding boundary condition complications. It also has low memory requirements and shows good performance. We illustrate an application of this solver: global human dispersal in the late Pleistocene, modeled via growth and diffusion over geographical maps of paleovegetation and paleoclimate.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a stable and efficient finite - difference solver for simulating the population growth and diffusion phenomena described by the Fisher - Kolmogorov - Petrovski - Piskunov (Fisher/KPP) equation on a two - dimensional geographical map. Specifically, this method is suitable for modeling the growth and diffusion of population on land maps, especially the global human diffusion in the late Pleistocene.
### Main Problems and Solutions
1. **Requirement for Numerical Simulation**:
- The paper points out that for population dynamics modeling on large - scale spatio - temporal scales (such as the diffusion of modern humans out of Africa or the diffusion and extinction of Neanderthals), reliable numerical methods are required to explain the local and global patterns of genetic, phenotypic, and cultural variation.
2. **Model Selection**:
- The Fisher/KPP equation is a classical reaction - diffusion equation, initially used to describe the spread of favorable gene mutations in an existing population in a one - dimensional case. Kolmogorov et al. extended it to the two - dimensional case and it has been widely used in the research of population diffusion and growth.
3. **Improvement of Numerical Method**:
- The author proposes a semi - implicit finite - difference scheme based on Godunov splitting, which has the following characteristics:
- High stability: The method is semi - implicit, so it has good stability.
- Relatively high precision: It is approximately \(O(\Delta x^2)+O(\Delta t^2)\), excluding the influence of boundary condition complexity.
- Low memory requirement: It is suitable for large - grid computing and has good performance.
4. **Application Example**:
- The author shows an application example of this solver: by performing growth and diffusion simulations on maps of ancient vegetation and ancient climate, the global human diffusion process in the late Pleistocene is reproduced.
### Mathematical Model
The Fisher/KPP equation can be expressed as:
\[
\frac{\partial p}{\partial t}=\lambda\left(1 - \frac{p}{K}\right)p + c\nabla^{2}p
\]
where:
- \(p\) is the population density,
- \(\lambda\) is the growth rate,
- \(K\) is the environmental carrying capacity,
- \(c\) is the diffusion coefficient.
Under certain simplified conditions, the equation can be further simplified to:
\[
\frac{\partial u}{\partial t}=(1 - u)u+\frac{1}{2}\nabla^{2}u
\]
### Conclusion
The paper successfully develops a stable and efficient finite - difference solver, which is suitable for simulating the Fisher/KPP equation on a two - dimensional domain of any shape (such as a geographical map). This method can not only handle the spatially - dependent environmental carrying capacity \(K(x)\), but also be extended to handle \(K(x, t)\) which is dependent on both time and space. Through practical applications, the effectiveness and rationality of this method in simulating the global diffusion process of humans out of Africa are verified.