Discrete Eulerian model for population genetics and dynamics under flow

Abstract

Marine species reproduce and compete while being advected by turbulent flows. It is largely unknown, both theoretically and experimentally, how population dynamics and genetics are changed by the presence of fluid flows. Discrete agent-based simulations in continuous space allow for accurate treatment of advection and number fluctuations, but can be computationally expensive for even modest organism densities. In this report, we propose an algorithm to overcome some of these challenges. We first provide a thorough validation of the algorithm in one and two dimensions without flow. Next, we focus on the case of weakly compressible flows in two dimensions. This models organisms such as phytoplankton living at a specific depth in the three-dimensional, incompressible ocean experiencing upwelling and/or downwelling events. We show that organisms born at sources in a two-dimensional time-independent flow experience an increase in fixation probability.

arXiv comment9 pages, 9 figures

References in corpus (10)

• Nonequilibrium Critical Phenomena and Phase Transitions into Absorbing States
• Shift in the velocity of a front due to a cut-off
• Genetic Demixing and Evolutionary Forces in the One-Dimensional Stepping Stone Model
• Integration of Langevin Equations with Multiplicative Noise and Viability of Field Theories for Absorbing Phase Transitions
• Population dynamics at high Reynolds number
• Population genetics in compressible flows
• Growth, competition and cooperation in spatial population genetics
• Fisher equation with turbulence in one dimension
• Population dynamics in compressible flows
• Fixation Probabilities in Weakly Compressible Fluid Flows

Cited by in corpus (2)

• Spatial population genetics with fluid flow
• Oceanic frontal divergence alters phytoplankton competition and distribution