How adaptation to food resources and death rates shape oscillatory dynamics in a microbial population

Microbes constantly interact with their environment by depleting and transforming food sources. Theoretical studies have mainly focused on Lotka-Volterra models, which do not account for food source dynamics. In contrast, consumer-resource models, which consider food source dynamics, are less explored. In particular, it is still unclear what physical mechanisms control oscillatory dynamics at a single population level, a phenomenon which can only be captured by a consumer-resource model. Here, we present a minimalistic consumer-resource model of a single microbial population with growth and death dynamics, consuming a continuously replenishing substrate. Our model reveals that decaying oscillations can occur around steady state if and only if the timescale of microbial adaptation to food supply changes exceeds the death timescale. This interplay of timescales allows us to rationalize the emergence of oscillatory dynamics when adding various biophysical ingredients to the model. We find that microbial necromass recycling or complementary use of multiple food sources reduces the parameter range for oscillations and increases the decay rate of oscillations. Requiring multiple simultaneous food sources has the opposite effect. Essentially, facilitating growth reduces the likelihood of oscillations around a fixed point. We further demonstrate that such damped oscillatory behavior is correlated with persistent oscillatory behavior in a noisy environment. We hope our work will motivate further investigations of consumer-resource models to improve descriptions of environments where food source distributions vary in space and time.

💡 Research Summary

The paper presents a minimalist consumer‑resource (CR) model that captures the dynamics of a single microbial population feeding on a continuously supplied abiotic substrate. Traditional Lotka‑Volterra (LV) approaches treat all species as biotic and ignore resource dynamics, which limits their ability to describe phenomena such as resilience, niche expansion, or nutrient cycling. By contrast, the authors formulate the system with two state variables: microbial density X(t) and substrate concentration S(t). Growth follows a rate r µ(S) where µ(S) is a monotonic increasing function (e.g., Monod), while death occurs at a constant rate λ. Substrate is supplied at rate n K_S and is consumed proportionally to microbial growth via a conversion factor γ.

Through nondimensionalization (x = X K_S/γ, s = S/K_S, τ = λt) the equations reduce to a system governed solely by the dimensionless reproduction‑to‑death ratio r/λ and the dimensionless substrate inflow n/λ. The unique non‑trivial steady state is given by (x*, s*) = (n/λ, µ⁻¹(λ/r)). Linear stability analysis yields a Jacobian whose eigenvalues are λ̃± = −β ± √(β²−4β) ⁄ 2, where the key parameter β = (r/λ) µ′(s*) x* encapsulates the ratio of two intrinsic timescales: the adaptation time to changes in resource availability (τ_food) and the recovery time set by death or predation (τ_death). When τ_food ≤ τ_death (β ≤ 1) the system quickly tracks resource fluctuations and converges monotonically. When τ_food ≥ τ_death (β ≥ 1) the two processes become out‑of‑phase, producing damped oscillations around the steady state. A more precise classification emerges: β ≤ 2 yields under‑damped oscillations with observable cycles, 2 < β ≤ 4 gives a “damped” regime where oscillations are heavily suppressed, and β > 4 leads to an overdamped, monotonic relaxation.

The authors then explore how various biologically realistic extensions modify β and thus the propensity for oscillations:

1. Necromass recycling – dead cells are decomposed back into substrate, adding an extra source term to the substrate equation. This increases µ′(s*) and effectively shortens τ_food, lowering β, shrinking the oscillatory region, and accelerating decay of any transients.

2. Multiple resources – two scenarios are considered. When resources are complementary (either can support growth), the effective resource supply is higher, again reducing τ_food and β, suppressing oscillations. Conversely, when resources are simultaneously required (both must be present), the system’s response to fluctuations slows, τ_food grows, β rises, and the oscillatory window expands.

3. Substrate removal/dilution – introducing a chemostat‑like outflow for substrate shortens the resource residence time, decreasing τ_food and β, thereby damping oscillations.

Stochastic simulations demonstrate that in the deterministic regime where β ≤ 4, adding noisy inflow can sustain persistent oscillations, akin to limit‑cycle behavior observed experimentally in continuous‑culture bioreactors and natural microbial ecosystems (soil, marine). This links the deterministic timescale ratio to observable variability under realistic environmental fluctuations.

Overall, the study provides a clear mechanistic insight: the emergence and strength of damped oscillations in a single‑species consumer‑resource system are governed by a single dimensionless parameter β that captures the competition between resource‑adaptation speed and mortality‑driven recovery. By adjusting measurable quantities—growth rate r, death rate λ, and substrate inflow n—researchers can predict whether a microbial community will quickly settle to equilibrium or exhibit prolonged transient oscillations. The framework also clarifies how biophysical processes such as necromass recycling, multi‑resource dependence, and external dilution modulate these timescales, offering a quantitative guide for both experimental design and ecological modeling of microbial dynamics in spatially and temporally heterogeneous environments.