Macroscopic limit from a structured population model to the Kirkpatrick-Barton model
Gaël Raoul
TL;DR
The work analyzes a population density $n(t,x,y)$ structured by space $x$ and trait $y$, governed by a diffusion in $x$ and a reproducing, trait-diffusing operator in $y$ (the SIM). It first establishes global existence and uniqueness of SIM solutions with uniform tail and regularity controls, then leverages a Tanaka-type Wasserstein contraction to show that for large $\gamma$ the trait distribution concentrates near a Gaussian with mean $Z(t,x)$ and variance $A$, i.e. $n\approx N\,\Gamma_A(\cdot-Z)$. The macroscopic fields $(N,Z)$ satisfy a closed Kirkpatrick-Barton system (KBM) up to explicit corrections $\varphi_{N,\gamma}$ and $\varphi_{Z,\gamma}$, with quantitative bounds $W_2(\tilde n, \Gamma_A(\cdot-Z)) \le C/\gamma^{\theta}$ after an initial layer, and convergence of $(N,Z)$ to the KBM solution as $\gamma\to\infty$. This provides a rigorous link between a mesoscopic, sexually reproducing structured population model and the classical ecological KBM, with implications for understanding climate-driven trait dynamics and range shifts.
Abstract
We consider an ecology model in which the population is structured by a spatial variable and a phenotypic trait. The model combines a parabolic operator on the spatial variable with a kinetic operator on the trait variable. We prove the existence of solutions to that model, and show that these solutions are unique. The kinetic operator present in the model, that represents the effect of sexual reproductions, satisfies a Tanaka-type inequality: it implies a contraction of the Wasserstein distance in the space of phenotypic traits. We combine this contraction argument with parabolic estimates controlling the spatial regularity of solutions to prove the convergence of the population size and the mean phenotypic trait to solutions of the Kirkpatrick-Barton model, which is a well-established model in evolutionary ecology. Specifically, at high reproductive rates, we provide explicit convergence estimates for the moments of solutions of the kinetic model.
