Stationary moments, diffusion limits, and extinction times for logistic growth with random catastrophes

Abstract

A central problem in population ecology is understanding the consequences of stochastic fluctuations. Analytically tractable models with Gaussian driving noise have led to important, general insights, but they fail to capture rare, catastrophic events, which are increasingly observed at scales ranging from global fisheries to intestinal microbiota. Due to mathematical challenges, growth processes with random catastrophes are less well characterized and it remains unclear how their consequences differ from those of Gaussian processes. In the face of a changing climate and predicted increases in ecological catastrophes, as well as increased interest in harnessing microbes for therapeutics, these processes have never been more relevant. To better understand them, I revisit here a differential equation model of logistic growth coupled to density-independent catastrophes that arrive as a Poisson process, and derive new analytic results that reveal its statistical structure. First, I derive exact expressions for the model's stationary moments, revealing a single effective catastrophe parameter that largely controls low order statistics. Then, I use weak convergence theorems to construct its Gaussian analog in a limit of frequent, small catastrophes, keeping the stationary population mean constant for normalization. Numerically computing statistics along this limit shows how they transform as the dynamics shifts from catastrophes to diffusions, enabling quantitative comparisons. For example, the mean time to extinction increases monotonically by orders of magnitude, demonstrating significantly higher extinction risk under catastrophes than under diffusions. Together, these results provide insight into a wide range of stochastic dynamical systems important for ecology and conservation.

Higher-Order Interactions Dampen Pairwise Competition in the Zebrafish Gut Microbiome

“…Because of strong aggregation observed in nearly all bacterial species, most individual bacteria residing in the bulk of clusters will not directly interact with other species, leading to interactions that are sublinear in population size, suggesting a logarithmic or α < 1 power law functional form. Furthermore, stochastic dynamics can be mapped onto robust average properties for populations ( 39 , 61 ). Therefore, it is reasonable to make use of simple models not as rigorous descriptions of the system but as approximations whose parameters characterize effective behaviors.…”

Experimental bacterial adaptation to the zebrafish gut reveals a primary role for immigration

“…We note, however, that additional processes are present and could be elaborated in future models if warranted by relevant data. Bacterial growth, death, and expulsion from the gut [32] are inherently stochastic, and future stochastic models could quantitatively link these processes to measurable population statistics, either through brute force computation or potentially by making use of recently developed analytic tools [63]. Additional processes neglected in this work include the transition time between lag phase and exponential growth upon entering the gut environment (itself stochastic) and interindividual variation in growth rates.…”

Cost‐efficient monitoring of continuous‐time stochastic processes based on discrete observations

“…We use the simplest population dynamics model: where X t is the algae biomass per unit area of riverbed and r > 0 (1/day) is the population growth rate. This is the simplest logistic growth model subject to sudden population decrease 75 …”