The Allee Effect in Compressible Flows

Technical Summary: The Allee Effect in Compressible Flows

Problem Statement
Microbial populations in marine environments often inhabit thin near-surface layers where they are advected by turbulent flows. While these flows are fundamentally three-dimensional and incompressible, the constraint of organisms to a two-dimensional interface renders the effective flow field compressible. Previous studies have established that for populations with simple logistic growth, this effective compressibility reduces the carrying capacity. However, many marine microbes exhibit an Allee effect, where per-capita growth rates increase at low densities due to cooperative mechanisms (e.g., quorum sensing). This paper investigates how the interplay between compressible turbulent transport and cooperative growth (specifically strong Allee effects) influences population survival, stability, and extinction thresholds.

Methodology
The authors model the population density $c(x,t)$ using a two-dimensional convection–diffusion equation with a local Allee growth term:
$$\partial_t c + \nabla \cdot (uc) = D\nabla^2 c + \mu c(1-c)(c-a)$$
where $u$ is a compressible flow field, $D$ is diffusivity, $\mu$ is the growth rate, and $a$ represents the Allee strength.

To generate the flow field, the authors solve the three-dimensional incompressible Navier–Stokes equations with random large-scale forcing. They extract a single mid-plane layer to construct a projected 2D flow field with a measured effective compressibility $\kappa \approx 0.17$. The system is analyzed across varying Damköhler numbers ($Da = \tau_{eddy}/\tau_{rep}$), which compares the eddy turnover time to the organism's reproduction time.

The study employs:

1. Direct Numerical Simulations: Tracking concentration fields and long-time averaged carrying capacities $\langle Z \rangle$ across the $(a, Da)$ parameter space.
2. Perturbative Analysis: Expanding the concentration field for large $Da$ (fast reproduction) to derive corrections to the carrying capacity.
3. Moment Analysis: For small $Da$ (slow reproduction), the population density is approximated as a passive tracer density $P$. The authors derive analytical expressions for the stationary carrying capacity by balancing growth and decay terms using the first three moments of the probability density function ($\langle P \rangle, \langle P^2 \rangle, \langle P^3 \rangle$).
4. Compressibility Variation: A continuous family of flow fields is constructed by mixing incompressible and compressible components to isolate the effect of the compressibility parameter $\kappa$.

Key Results

• Large Damköhler Limit ($Da \gg 1$): When reproduction is fast relative to turbulent mixing, the population generally tracks the stable fixed point $c=1$. However, compressibility introduces a perturbative reduction in carrying capacity proportional to $\tau_{rep}^2$ and the Allee strength $a$. Crucially, for $a > 1/2$, the empty state ($c=0$) becomes energetically favorable. In this regime, compressibility can nucleate depleted domains that expand as "pushed waves," leading to total population extinction even if the initial condition is uniform and high.

• Small Damköhler Limit ($Da \ll 1$): When reproduction is slow, the population density becomes slaved to the clustering statistics of passive tracers in the compressible flow. The authors derive an analytical condition for the stationary carrying capacity:
  $$\langle Z \rangle_{\pm} = \frac{\langle P \rangle \langle P^2 \rangle}{2\langle P^3 \rangle} \left[ (1+a) \pm \sqrt{(1+a)^2 - 4a\gamma} \right]$$
  where $\gamma = \langle P \rangle \langle P^3 \rangle / \langle P^2 \rangle^2 \geq 1$.
  This analysis reveals a critical Allee strength $a_{crit}^-$ below which the population survives and above which it goes extinct. Unlike the well-mixed case where $a_{crit}=1$, compressibility ($\gamma > 1$) drastically lowers this threshold. For the simulated flow ($\gamma \approx 186$), the critical Allee strength drops to $a_{crit}^- \approx 0.1$.

• Extinction Transition: The transition from survival to extinction is characterized as a first-order phase transition. There is a discontinuous jump in the carrying capacity at the critical Allee strength. The extinction boundary is determined by the statistics of the compressible flow (specifically the multifractal clustering of passive particles) rather than just the intrinsic growth parameters.

• Role of Compressibility: Increasing the effective compressibility $\kappa$ systematically suppresses survival. Even modest increases in $\kappa$ significantly increase the higher moments of the tracer distribution, thereby lowering the critical Allee strength required for extinction.

Significance and Claims
The paper demonstrates that compressible turbulent transport qualitatively alters the stability landscape of populations exhibiting cooperative growth. While compressibility merely reduces the carrying capacity for logistic growth, it can become fatal for species with strong Allee effects by inducing extinction at Allee strengths that would otherwise be sustainable in well-mixed or incompressible environments.

The authors claim that for such populations, turbulent transport cannot be viewed merely as a renormalization of mixing or diffusivity. Instead, the compressibility of the flow sets an upper limit on the Allee strength a species can tolerate. This suggests that the survival of marine microbes relying on quorum sensing or cooperative metabolism may be critically constrained by the compressible nature of their transport in thin surface layers. The work provides explicit analytical criteria for extinction thresholds based on flow statistics, linking the physics of passive scalar clustering directly to biological population dynamics.