•Modeling and analysis of predator–prey system with disease in prey.•Prey individuals close to herd’s boundary are subject to hunting by predators.•Total collapse caused by heteroclinic connection ...between two prey-only saddle equilibria.•Torus bifurcation distorted by saddle limit cycle.•Prey can exist for larger mortality rates of predators by weakening of prey by infection.
In this paper we analyse a predator–prey model where the prey population shows group defense and the prey individuals are affected by a transmissible disease. The resulting model is of the Rosenzweig–MacArthur predator–prey type with an SI (susceptible-infected) disease in the prey. Modeling prey group defense leads to a square root dependence in the Holling type II functional for the predator–prey interaction term. The system dynamics is investigated using simulations, classical existence and asymptotic stability analysis and numerical bifurcation analysis. A number of bifurcations, such as transcritical and Hopf bifurcations which occur commonly in predator–prey systems will be found. Because of the square root interaction term there is non-uniqueness of the solution and a singularity where the prey population goes extinct in a finite time. This results in a collapse initiated by extinction of the healthy or susceptible prey and thereafter the other population(s). When also a positive attractor exists this leads to bistability similar to what is found in predator–prey models with a strong Allee effect. For the two-dimensional disease-free (i.e. the purely demographic) system the region in the parameter space where bistability occurs is marked by a global bifurcation. At this bifurcation a heteroclinic connection exists between saddle prey-only equilibrium points where a stable limit cycle together with its basin of attraction, are destructed. In a companion paper (Gimmelli et al., 2015) the same model was formulated and analysed in which the disease was not in the prey but in the predator. There we also observed this phenomenon. Here we extend its analysis using a phase portrait analysis. For the three-dimensional ecoepidemic predator–prey system where the prey is affected by the disease, also tangent bifurcations including a cusp bifurcation and a torus bifurcation of limit cycles occur. This leads to new complex dynamics. Continuation by varying one parameter of the emerging quasi-periodic dynamics from a torus bifurcation can lead to its destruction by a collision with a saddle-cycle. Under other conditions the quasi-periodic dynamics changes gradually in a trajectory that lands on a boundary point where the prey go extinct in finite time after which a total collapse of the three-dimensional system occurs.
The majority of studies on environmental change focus on the response of single species and neglect fundamental biotic interactions, such as mutualism, competition, predation, and parasitism, which ...complicate patterns of species persistence. Under global warming, disruption of community interactions can arise when species differ in their sensitivity to rising temperature, leading to mismatched phenologies and/or dispersal patterns. To study species persistence under global climate change, it is critical to consider the ecology and evolution of multispecies interactions; however, the sheer number of potential interactions makes a full study of all interactions unfeasible. One mechanistic approach to solving the problem of complicated community context to global change is to (i) define strategy groups of species based on life-history traits, trophic position, or location in the ecosystem, (ii) identify species involved in key interactions within these groups, and (iii) determine from the interactions of these key species which traits to study in order to understand the response to global warming. We review the importance of multispecies interactions looking at two trait categories: thermal sensitivity of metabolic rate and associated life-history traits and dispersal traits of species. A survey of published literature shows pronounced and consistent differences among trophic groups in thermal sensitivity of life-history traits and in dispersal distances. Our approach increases the feasibility of unraveling such a large and diverse set of community interactions, with the ultimate goal of improving our understanding of community responses to global warming.
We introduce a compartmental host-vector model for dengue with two viral strains, temporary cross-immunity for the hosts, and possible secondary infections. We study the conditions on existence of ...endemic equilibria where one strain displaces the other or the two virus strains co-exist. Since the host and vector epidemiology follow different time scales, the model is described as a slow-fast system. We use the geometric singular perturbation technique to reduce the model dimension. We compare the behaviour of the full model with that of the model with a quasi-steady approximation for the vector dynamics. We also perform numerical bifurcation analysis with parameter values from the literature and compare the bifurcation structure to that of previous two-strain host-only models.
Group defense in prey and hunting cooperation in predators are two important ecological phenomena and can occur concurrently. In this article, we consider cooperative hunting in generalist predators ...and group defense in prey under a mathematical framework to comprehend the enormous diversity the model could capture. To do so, we consider a modified Holling-Tanner model where we implement Holling type IV functional response to characterize grazing pattern of predators where prey species exhibit group defense. Additionally, we allow a modification in the attack rate of predators to quantify the hunting cooperation among them. The model admits three boundary equilibria and up to three coexistence equilibrium points. The geometry of the nontrivial prey and predator nullclines and thus the number of coexistence equilibria primarily depends on a specific threshold of the availability of alternative food for predators. We use linear stability analysis to determine the types of hyperbolic equilibrium points and characterize the non-hyperbolic equilibrium points through normal form and center manifold theory. Change in the model parameters leading to the occurrences of a series of local bifurcations from non-hyperbolic equilibrium points, namely, transcritical, saddle-node, Hopf, cusp and Bogdanov-Takens bifurcation; there are also occurrences of global bifurcations such as homoclinic bifurcation and saddle-node bifurcation of limit cycles. We observe two interesting closed ‘bubble’ form induced by global bifurcations due to change in the strength of hunting cooperation and the availability of alternative food for predators. A three dimensional bifurcation diagram, concerning the original system parameters, captures how the alternation in model formulation induces gradual changes in the bifurcation scenarios. Our model highlights the stabilizing effects of group or gregarious behaviour in both prey and predator, hence supporting the predator-herbivore regulation hypothesis. Additionally, our model highlights the occurrence of “saltatory equilibria" in ecological systems and capture the dynamics observed for lion-herbivore interactions.
We study a predator–prey model with different characteristic time scales for the prey and predator populations, assuming that the predator dynamics is much slower than the prey one. Geometrical ...Singular Perturbation theory provides the mathematical framework for analyzing the dynamical properties of the model. This model exhibits a Hopf bifurcation and we prove that when this bifurcation occurs, a canard phenomenon arises. We provide an analytic expression to get an approximation of the bifurcation parameter value for which a maximal canard solution occurs. The model is the well-known Rosenzweig–MacArthur predator–prey differential system. An invariant manifold with a stable and an unstable branches occurs and a geometrical approach is used to explicitly determine a solution at the intersection of these branches. The method used to perform this analysis is based on Blow-up techniques. The analysis of the vector field on the blown-up object at an equilibrium point where a Hopf bifurcation occurs with zero perturbation parameter representing the time scales ratio, allows to prove the result. Numerical simulations illustrate the result and allow to see the canard explosion phenomenon.
The motivation for the research reported in this paper comes from modeling the spread of vector-borne virus diseases. To study the role of the host versus vector dynamics and their interaction we use ...the susceptible-infected-removed (SIR) host model and the susceptible-infected (SI) vector model. When the vector dynamical processes occur at a faster scale than those in the host-epidemics dynamics, we can use a time-scale argument to reduce the dimension of the model. This is often implemented as a quasi steady-state assumption (qssa) where the slow varying variable is set at equilibrium and an ode equation is replaced by an algebraic equation. Singular perturbation theory will appear to be a useful tool to perform this derivation. An asymptotic expansion in the small parameter that represents the ratio of the two time scales for the dynamics of the host and vector is obtained using an invariant manifold equation. In the case of a susceptible-infected-susceptible (SIS) host model this algebraic equation is a hyperbolic relationship modeling a saturated incidence rate. This is similar to the Holling type II functional response (Ecology) and the Michaelis-Menten kinetics (Biochemistry). We calculate the value for the force of infection leading to an endemic situation by performing a bifurcation analysis. The effect of seasonality is studied where the force of infection changes sinusoidally to model the annual fluctuations of the vector population. The resulting non-autonomous system is studied in the same way as the autonomous system using bifurcation analysis.
•Ecoepidemic model for the herd behavior of prey attacked by specialist predators.•Dynamics includes diseased prey and feeding saturation of predators.•Several types of bifurcations are discovered.
...In this paper we consider a model for the herd behavior of prey, that are subject to attacks by specialist predators. The latter are affected by a transmissible disease. With respect to other recently introduced models of the same nature, we focus here our attention to the possible feeding satiation phenomenon. The system dynamics is thoroughly investigated, to show the occurrence of several types of bifurcations. In addition to the transcritical and Hopf bifurcation that occur commonly in predator–prey system also a zero-Hopf and a global bifurcation occur. The Hopf and the global bifurcation occur only in the disease-free (so purely demographic) system. The latter is a heteroclinic connection for the between saddle equilibrium points where a stable limit cycle is disrupted and where the system disease-free collapses while in a parameter space region the endemic system exists stably.
In many countries in Asia and South-America dengue fever (DF) and dengue hemorrhagic fever (DHF) has become a substantial public health concern leading to serious social-economic costs. Mathematical ...models describing the transmission of dengue viruses have focussed on the so-called antibody-dependent enhancement (ADE) effect and temporary cross-immunity trying to explain the irregular behavior of dengue epidemics by analyzing available data. However, no systematic investigation of the possible dynamical structures has been performed so far. Our study focuses on a seasonally forced (non-autonomous) model with temporary cross-immunity and possible secondary infection, motivated by dengue fever epidemiology. The notion of at least two different strains is needed in a minimalistic model to describe differences between primary infections, often asymptomatic, and secondary infection, associated with the severe form of the disease. We extend the previously studied non-seasonal (autonomous) model by adding seasonal forcing, mimicking the vectorial dynamics, and a low import of infected individuals, which is realistic in the dynamics of dengue fever epidemics. A comparative study between three different scenarios (non-seasonal, low seasonal and high seasonal with a low import of infected individuals) is performed. The extended models show complex dynamics and qualitatively a good agreement between empirical DHF monitoring data and the obtained model simulation. We discuss the role of seasonal forcing and the import of infected individuals in such systems, the biological relevance and its implications for the analysis of the available dengue data. At the moment only such minimalistic models have a chance to be qualitatively understood well and eventually tested against existing data. The simplicity of the model (low number of parameters and state variables) offer a promising perspective on parameter values inference from the DHF case notifications.
► We extend the previously studied non-seasonal model by adding seasonal forcing and import of infected. ► A comparative study between three different scenarios (non-seasonal, seasonal and seasonal with import) is performed. ► The complete analysis of the extended models shows complex dynamics. ► The simplicity of the model offer a promising perspective on parameter values inference from empirical data.
Mathematical models have a long history in epidemiological research, and as the COVID-19 pandemic progressed, research on mathematical modeling became imperative and very influential to understand ...the epidemiological dynamics of disease spreading.
Mathematical models describing dengue fever epidemiological dynamics are found back from 1970. Dengue fever is a viral mosquito-borne infection caused by four antigenically related but distinct serotypes (DENV-1 to DENV-4). With 2.5 billion people at risk of acquiring the infection, it is a major international public health concern. Although most of the cases are asymptomatic or mild, the disease immunological response is complex, with severe disease linked to the antibody-dependent enhancement (ADE) - a disease augmentation phenomenon where pre-existing antibodies to previous dengue infection do not neutralize but rather enhance the new infection. Here, we present a 10-year systematic review on mathematical models for dengue fever epidemiology. Specifically, we review multi-strain frameworks describing host-to-host and vector-host transmission models and within-host models describing viral replication and the respective immune response.
Following a detailed literature search in standard scientific databases, different mathematical models in terms of their scope, analytical approach and structural form, including model validation and parameter estimation using empirical data, are described and analyzed.
Aiming to identify a consensus on infectious diseases modeling aspects that can contribute to public health authorities for disease control, we revise the current understanding of epidemiological and immunological factors influencing the transmission dynamics of dengue. This review provide insights on general features to be considered to model aspects of real-world public health problems, such as the current epidemiological scenario we are living in.
•10-year systematic review on mathematical models for dengue fever epidemiology.•Detailed review of multi-strain frameworks describing host-to-host dengue transmission models.•Detailed review of vector-host dengue transmission.•Detailed review of within-host dengue models.•Comprehensive description of dengue and COVID-19 co-infection modeling.
•Four predator–prey models in chemostat are built from Dynamic Energy Budget theory.•Structural sensitivity to functional response formulation is investigated.•The less detailed metabolic model ...(Monod) leads to higher structural sensitivity.•Structural sensitivity is lower in mass-balanced models including maintenance.
Many current issues in ecology require predictions made by mathematical models, which are built on somewhat arbitrary choices. Their consequences are quantified by sensitivity analysis to quantify how changes in model parameters propagate into an uncertainty in model predictions. An extension called structural sensitivity analysis deals with changes in the mathematical description of complex processes like predation. Such processes are described at the population scale by a specific mathematical function taken among similar ones, a choice that can strongly drive model predictions. However, it has only been studied in simple theoretical models. Here, we ask whether structural sensitivity is a problem of oversimplified models. We found in predator–prey models describing chemostat experiments that these models are less structurally sensitive to the choice of a specific functional response if they include mass balance resource dynamics and individual maintenance. Neglecting these processes in an ecological model (for instance by using the well-known logistic growth equation) is not only an inappropriate description of the ecological system, but also a source of more uncertain predictions.