skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Unifying deterministic and stochastic ecological dynamics via a landscape-flux approach
The frequency distributions can characterize the population-potential landscape related to the stability of ecological states. We illustrate the practical utility of this approach by analyzing a forest–savanna model. Savanna and forest states coexist under certain conditions, consistent with past theoretical work and empirical observations. However, a grassland state, unseen in the corresponding deterministic model, emerges as an alternative quasi-stable state under fluctuations, providing a theoretical basis for the appearance of widespread grasslands in some empirical analyses. The ecological dynamics are determined by both the population-potential landscape gradient and the steady-state probability flux. The flux quantifies the net input/output to the ecological system and therefore the degree of nonequilibriumness. Landscape and flux together determine the transitions between stable states characterized by dominant paths and switching rates. The intrinsic potential landscape admits a Lyapunov function, which provides a quantitative measure of global stability. We find that the average flux, entropy production rate, and free energy have significant changes near bifurcations under both finite and zero fluctuation. These may provide both dynamical and thermodynamic origins of the bifurcations. We identified the variances in observed frequency time traces, fluctuations, and time irreversibility as kinematic measures for bifurcations. This framework opens the way to characterize ecological systems globally, to uncover how they change among states, and to quantify the emergence of quasi-stable states under stochastic fluctuations.  more » « less
Award ID(s):
1951358 1951394 1951385
PAR ID:
10249282
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Proceedings of the National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
118
Issue:
24
ISSN:
0027-8424
Page Range / eLocation ID:
Article No. e2103779118
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Predator‐prey models, such as the Leslie‐Gower model, are essential for understanding population dynamics and stability within ecosystems. These models help explain the balance between species under natural conditions, but the inclusion of factors like the Allee effect and intraspecific competition adds complexity and realism to these interactions, enhancing our ability to predict system behavior under stress. To detect early indicators of population collapse, this study investigates the intricate dynamics of a modified Leslie‐Gower predator‐prey model with both Allee effect and intraspecific competition. We analyze the existence and stability of equilibria, as well as bifurcation phenomena, including saddle‐node bifurcations of codimension 2, Hopf bifurcations of codimension 2, and Bogdanov‐Takens bifurcations of codimension at least 4. Detailed transitions between bifurcation curves–specifically saddle‐node, Hopf, homoclinic, and limit cycle bifurcations–are also examined. We observe a novel transition phenomenon, where a system jumps from saddle‐node bifurcation to homoclinic and limit cycle bifurcations. This suggests that burst oscillations may serve as an early warning of system collapse rather than simply a tipping point. Our findings indicate that moderate levels of intraspecific competition or Allee effect support coexistence of both populations, while excessive levels may destabilize the entire biological system, leading to collapse. These insights offer valuable implications for ecological management and the early detection of risks in population dynamics. 
    more » « less
  2. In this work, we propose a predator-prey system with a Holling type Ⅱ functional response and study its dynamics when the prey exhibits vigilance behavior to avoid predation and predators exhibit cooperative hunting. We provide conditions for existence and the local and global stability of equilibria. We carry out detailed bifurcation analysis and find the system to experience Hopf, saddle-node, and transcritical bifurcations. Our results show that increased prey vigilance can stabilize the system, but when vigilance levels are too high, it causes a decrease in the population density of prey and leads to extinction. When hunting cooperation is intensive, it can destabilize the system, and can also induce bi-stability phenomenon. Furthermore, it can reduce the population density of both prey and predators and also change the stability of a coexistence state. We provide numerical experiments to validate our theoretical results and discuss ecological implications. 
    more » « less
  3. One of the most intriguing features of biological systems is their ability to regulate the steady-state fluxes of the underlying biochemical reactions; however, the regulatory mechanisms and their physicochemical properties are not fully understood. Fundamentally, flux regulation can be explained with a chemical kinetic formalism describing the transitions between discrete states, with the reaction rates defined by an underlying free energy landscape. Which features of the energy landscape affect the flux distribution? Here we prove that the ratios of the steady-state fluxes of quasi–first-order biochemical processes are invariant to energy perturbations of the discrete states and are only affected by the energy barriers. In other words, the nonequilibrium flux distribution is under kinetic and not thermodynamic control. We illustrate the generality of this result for three biological processes. For the network describing protein folding along competing pathways, the probabilities of proceeding via these pathways are shown to be invariant to the stability of the intermediates or to the presence of additional misfolded states. For the network describing protein synthesis, the error rate and the energy expenditure per peptide bond is proven to be independent of the stability of the intermediate states. For molecular motors such as myosin-V, the ratio of forward to backward steps and the number of adenosine 5′-triphosphate (ATP) molecules hydrolyzed per step is demonstrated to be invariant to energy perturbations of the intermediate states. These findings place important constraints on the ability of mutations and drug perturbations to affect the steady-state flux distribution for a wide class of biological processes. 
    more » « less
  4. Abstract We analyze a spatially extended version of a well-known model of forest-savanna dynamics, which presents as a system of nonlinear partial integro-differential equations, and study necessary conditions for pattern-forming bifurcations. Homogeneous solutions dominate the dynamics of the standard forest-savanna model, regardless of the length scales of the various spatial processes considered. However, several different pattern-forming scenarios are possible upon including spatial resource limitation, such as competition for water, soil nutrients, or herbivory effects. Using numerical simulations and continuation, we study the nature of the resulting patterns as a function of system parameters and length scales, uncovering subcritical pattern-forming bifurcations and observing significant regions of multistability for realistic parameter regimes. Finally, we discuss our results in the context of extant savanna-forest modeling efforts and highlight ongoing challenges in building a unifying mathematical model for savannas across different rainfall levels. 
    more » « less
  5. In this work, we propose and investigate a predator-prey model where the prey population is structured by sex and the predators (unstructured) depredate based on sex-bias. We provide conditions for the existence of equilibrium points and perform local stability analysis on them. We derive global stability conditions for the extinction state. We show the possible occurrence of Hopf and saddle-node bifurcations. Multiple Hopf bifurcations are observed as the sex-biased predation rate is varied. This variation also shows the opposite consequences in the densities of the sex-structured prey. Our results show that sex-biased predation can cause both stabilizing and destabilizing effects for certain parameter choices. It can also cause an imbalanced sex-ratio, which has ecological consequences. Furthermore when intraspecific competition among predators is minimized, it can lead to the extinction of prey. We discuss the ecological implications and application of our results to the biocontrol of invasive species susceptible to sex-biased predation. 
    more » « less