Coexistence and Extinction in Flow-Kick Systems: An invasion growth rate approach

TL;DR

This paper develops an invasion growth rate theory for flow-kick systems, providing criteria for species persistence or extinction in models with continuous and discrete dynamics, with applications to microbial, spatial, and Lotka-Volterra systems.

Contribution

It introduces a novel invasion growth rate framework for flow-kick models with state-dependent timing, linking Lyapunov exponents to species persistence and extinction.

Findings

Invasion growth rates characterize species persistence in flow-kick systems.

Permanence is determined by the signs of invasion growth rates and the structure of invasion graphs.

The theory is demonstrated through applications to microbial, spatial, and Lotka-Volterra models.

Abstract

Populations experience a complex interplay of continuous and discrete processes: continuous growth and interactions are punctuated by discrete reproduction events, dispersal, and external disturbances. These dynamics can be modeled by impulsive or flow-kick systems, where continuous flows alternate with instantaneous discrete changes. To study species persistence in these systems, an invasion growth rate theory is developed for flow-kick models with state-dependent timing of kicks and auxiliary variables. The invasion growth rates are Lyapunov exponents characterizing the average per-capita growth of species when rare. Two theorems are proven to characterize permanence i.e. the extinction set is a repellor. The first theorem uses Morse decompositions of the extinction set and requires that there exists a species with a positive invasion growth rate for every invariant measure supported on a component of the Morse decomposition. The second theorem uses invasion growth rates to define invasion graphs whose vertices correspond to communities and directed edges to potential invasions. Provided the invasion graph is acyclic, permanence is fully characterized by the signs of the invasion growth rates. Invasion growth rates are also used to identify the existence of extinction-bound trajectories and attractors that lie on the extinction set. The results are illustrated with three applications: (i) a microbial serial transfer model, (ii) a spatially structured consumer-resource model, and (iii) an empirically parameterized Lotka-Volterra model. Mathematical challenges and promising biological applications are discussed.