Optimal harvesting of prey–predator system with interval biological parameters: A bioeconomic model

The paper presents the study of one prey one predator harvesting model with imprecise biological parameters. Due to the lack of precise numerical information of the biological parameters such as prey population growth rate, predator population decay rate and predation coefficients, we consider the model with imprecise data as form of an interval in nature. Many authors have studied prey–predator harvesting model in different form, here we consider a simple prey–predator model under impreciseness and introduce parametric functional form of an interval and then study the model. We identify the equilibrium points of the model and discuss their stabilities. The existence of bionomic equilibrium of the model is discussed. We study the optimal harvest policy and obtain the solution in the interior equilibrium using Pontryagin’s maximum principle. Numerical examples are presented to support the proposed model.     

Dynamics analysis of a predator–prey system with harvesting prey and disease in prey species

In this paper, a predator–prey system with harvesting prey and disease in prey species is given. In the absence of time delay, the existence and stability of all equilibria are investigated. In the presence of time delay, some sufficient conditions of the local stability of the positive equilibrium and the existence of Hopf bifurcation are obtained by analysing the corresponding characteristic equation, and the properties of Hopf bifurcation are given by using the normal form theory and centre manifold theorem. Furthermore, an optimal harvesting policy is investigated by applying the Pontryagin's Maximum Principle. Numerical simulations are performed to support our analytic results.     

Bioeconomic harvesting of a prey–predator fishery

This paper deals with the problem of non-selective harvesting of a prey–predator system by using a reasonable catch-rate function instead of usual catch-per-unit-efforthypothesis. Here both the prey and the predator species obey the law of logistic growth. We have taken the predator functional response to prey density in such a form that each predator's functional response to the prey density approaches a constant as the prey population increases. Boundedness of the exploited system is examined. The existence of its steady states and their stability (local and global) are studied using Eigenvalue analysis. The existence of bionomic equilibria has been illustrated using a numerical example. The problem of determining the optimal harvesting policy is then solved by using Pontryagin's maximum principle.     

Mathematical analysis of a delayed stage-structured predator–prey model with impulsive diffusion between two predators territories

In most models of population dynamics, diffusion between two patches is assumed to be either continuous or discrete, but in reality, many species diffuse only during a single period, and diffusion often occurs in regular pulses. Further, in forest habitats, the highest-level predator species are restricted to a specific territory, but prey can impulsively move between territories. Therefore, in this paper, we consider a delayed stage-structured predator–prey model with impulsively diffusive prey between two patches; in the model, patches represent the territories of two different predator populations. Here, we analytically obtain the global attractivity condition of predator-extinction periodic solutions for the system by using the concepts of Hui and Chen (2005); a numerical simulation is also included to illustrate this result. Further, we establish permanence conditions for the coexistence of the species using the theory of impulsive delayed differential equations. Finally, we explore the possibilities of the permanence of the system by using the growth rates of immature predators and the impulse period as critical parameters, and we also obtain the parameters’ threshold limits using numerical experimentation.     

A nonlocal kinetic model for predator–prey interactions

We extend the aggregation model from Fetecau (2011) by adding a field of vision to individuals and by including a second species. The two species, assumed to have a predator–prey relationship, have dynamics governed by nonlocal kinetic equations that include advection and turning. The latter is the main mechanism for aggregation and orientation, which results from interactions among individuals of the same species as well as predator–prey relationships. We illustrate numerically a diverse set of predator–prey behaviors that can be captured by this model. We show that a prey’s escape outcome depends on the social interactions between its group members, the prey’s field of vision and the sophistication of the predator’s hunting strategies.     

Time delay can enhance spatio-temporal chaos in a prey–predator model

In this paper we explore how the time delay induced Hopf-bifurcation interacts with Turing instability to determine the resulting spatial patterns. For this study, we consider a delayed prey–predator model with Holling type-II functional response and intra-specific competition among the predators. Analytical criteria for the delay induced Hopf-bifurcation and for the delayed spatio-temporal model are provided with numerical example to validate the analytical results. Exhaustive numerical simulation reveals the appearance of three types of stationary patterns, cold spot, labyrinthine, mixture of stripe-spot and two non-stationary patterns, quasi-periodic and spatio-temporal chaotic patterns. The qualitative features of the patterns for the non-delayed and the delayed spatio-temporal model are the same but their occurrence is solely controlled by the temporal parameters, rate of diffusivity and magnitude of the time delay. It is evident that the magnitude of time delay parameter beyond the Hopf-bifurcation threshold mostly produces spatio-temporal chaotic patterns.     

Turing patterns in a predator–prey model with seasonality

Many ecological systems show striking non-homogeneous population distributions. Diffusion-driven instabilities are commonly studied as mechanisms of pattern formation in many fields of biology but only rarely in ecology, in part because some of the conditions seem quite restrictive for ecological systems. Seasonal variation is ubiquitous in temperate ecosystems, yet its effect on pattern formation has not yet been explored. We formulate and analyze an impulsive reaction–diffusion system for a resource and its consumer in a two-season environment. While the resource grows throughout the ‘summer’ season, the consumer reproduces only once per year. We derive conditions for diffusion-driven instability in the system, and we show that pattern formation is possible with a Beddington–DeAngelis functional response. More importantly, we find that a low overwinter survival probability for the resource enhances the propensity for pattern formation: diffusion-driven instability occurs even when the diffusion rates of prey and predator are comparable (although not when they are equal).

     

Global analysis of a predator–prey model with variable predator search rate

We consider a modified Holling-type II predator–prey model, based on the premise that the search rate of predators is dependent on the prey density, rather than constant. A complete analysis of the global behavior of the model is presented, and shows that the model exhibits a dichotomy similar to the classical Holling-type II model: either the coexistence steady state is globally stable; or it is unstable, and then a unique, globally stable limit cycle exists. We discuss the similarities, but also important differences between our model and the Holling-type II model. The main differences are that: 1. The paradox of enrichment which always occurs in the Holling-type II model, does not always occur here, and 2. Even when the paradox of enrichment occurs, predators can adapt by lowering their search rate, and effectively stabilize the system.

     

Structural sensitivity and resilience in a predator–prey model with density-dependent mortality

Numerous formulations with the same mathematical properties can be relevant to model a biological process. Different formulations can predict different model dynamics like equilibrium vs. oscillations even if they are quantitatively close (structural sensitivity). The question we address in this paper is: does the choice of a formulation affect predictions on the number of stable states? We focus on a predator–prey model with predator competition that exhibits multiple stable states. A bifurcation analysis is realized with respect to prey carrying capacity and species body mass ratio within range of values found in food web models. Bifurcation diagrams built for two type-II functional responses are different in two ways. First, the kind of stable state (equilibrium vs. oscillations) is different for 26.0–49.4% of the parameter values, depending on the parameter space investigated. Using generalized modelling, we highlight the role of functional response slope in this difference. Secondly, the number of stable states is higher with Ivlev's functional response for 0.1–14.3% of the parameter values. These two changes interact to create different model predictions if a parameter value or a state variable is altered. In these two examples of disturbance, Holling's disc equation predicts a higher system resilience. Indeed, Ivlev's functional response predicts that disturbance may trap the system into an alternative stable state that can be escaped from only by a larger alteration (hysteresis phenomena). Two questions arise from this work: (i) how much complex ecological models can be affected by this sensitivity to model formulation? and (ii) how to deal with these uncertainties in model predictions?     

Permanence of a predator–prey discrete system with Holling-IV functional response and distributed delays

A predator–prey discrete-time model with Holling-IV functional response and distributed delays is investigated in this paper. By using the comparison theorem of the difference equation and some analysis technique, some sufficient conditions are obtained for the permanence of the discrete predator–prey system. Two examples are given to illustrate the feasibility of the obtained result.     

Rich dynamics in a predator–prey model with both noise and periodic force

A spatial version of the predator–prey model with Holling III functional response, which includes some important factors such as external periodic forces, noise, and diffusion processes is investigated. For the model only with diffusion, it exhibits spiral waves in the two-dimensional space. However, combined with noise, it has the feature of chaotic patterns. Moreover, the oscillations become more obvious when the noise intensity is increased. Furthermore, the spatially extended system with external periodic forces and noise exhibits a resonant pattern and frequency-locking phenomena. These results may help us to understand the effects arising from the undeniable susceptibility to random fluctuations in the real ecosystems.     

Combined harvesting of a stage structured prey–predator model incorporating cannibalism in competitive environment

In this paper, we propose a prey–predator system with stage structure for predator. The proposed system incorporates cannibalism for predator populations in a competitive environment. The combined fishing effort is considered as control used to harvest the populations. The steady states of the system are determined and the dynamical behavior of the system is discussed. Local stability of the system is analyzed and sufficient conditions are derived for the global stability of the system at the positive equilibrium point. The existence of the Hopf bifurcation phenomenon is examined at the positive equilibrium point of the proposed system. We consider harvesting effort as a control parameter and subsequently, characterize the optimal control parameter in order to formulate the optimal control problem under the dynamic framework towards optimal utilization of the resource. Moreover, the optimal system is solved numerically to investigate the sustainability of the ecosystem using an iterative method with a Runge–Kutta fourth-order scheme. Simulation results show that the optimal control scheme can achieve sustainable ecosystem. Results are analyzed with the help of graphical illustrations.     

The complex dynamics of a diffusive prey–predator model with an Allee effect in prey

This paper investigates complex dynamics of a predator–prey interaction model that incorporates: (a) an Allee effect in prey; (b) the Michaelis–Menten type functional response between prey and predator; and (c) diffusion in both prey and predator. We provide rigorous mathematical results of the proposed model including: (1) the stability of non-negative constant steady states; (2) sufficient conditions that lead to Hopf/Turing bifurcations; (3) a prior estimates of positive steady states; (4) the non-existence and existence of non-constant positive steady states when the model is under zero-flux boundary condition. We also perform completed analysis of the corresponding ODE model to obtain a better understanding on effects of diffusion on the stability. Our analytical results show that the small values of the ratio of the prey's diffusion rate to the predator's diffusion rate are more likely to destabilize the system, thus generate Hopf-bifurcation and Turing instability that can lead to different spatial patterns. Through numerical simulations, we observe that our model, with or without Allee effect, can exhibit extremely rich pattern formations that include but not limit to strips, spotted patterns, symmetric patterns. In addition, the strength of Allee effects also plays an important role in generating distinct spatial patterns.     

Parametric analysis of the ratio-dependent predator–prey model

We present a complete parametric analysis of stability properties and dynamic regimes of an ODE model in which the functional response is a function of the ratio of prey and predator abundances. We show the existence of eight qualitatively different types of system behaviors realized for various parameter values. In particular, there exist areas of coexistence (which may be steady or oscillating), areas in which both populations become extinct, and areas of "conditional coexistence" depending on the initial values. One of the main mathematical features of ratio-dependent models, distinguishing this class from other predator-prey models, is that the Origin is a complicated equilibrium point, whose characteristics crucially determine the main properties of the model. This is the first demonstration of this phenomenon in an ecological model. The model is investigated with methods of the qualitative theory of ODEs and the theory of bifurcations. The biological relevance of the mathematical results is discussed both regarding conservation issues (for which coexistence is desired) and biological control (for which extinction is desired).     

Community-driven dispersal in an individual-based predator–prey model

We present a spatial, individual-based predator–prey model in which dispersal is dependent on the local community. We determine species suitability to the biotic conditions of their local environment through a time and space varying fitness measure. Dispersal of individuals to nearby communities occurs whenever their fitness falls below a predefined tolerance threshold. The spatiotemporal dynamics of the model is described in terms of this threshold. We compare this dynamics with the one obtained through density-independent dispersal and find marked differences. In the community-driven scenario, the spatial correlations in the population density do not vary in a linear fashion as we increase the tolerance threshold. Instead we find the system to cross different dynamical regimes as the threshold is raised. Spatial patterns evolve from disordered, to scale-free complex patterns, to finally becoming well-organized domains. This model therefore predicts that natural populations, the dispersal strategies of which are likely to be influenced by their local environment, might be subject to complex spatiotemporal dynamics.     

The Beverton–Holt model with periodic and conditional harvesting

In this theoretical study, we investigate the effect of different harvesting strategies on the discrete Beverton–Holt model in a deterministic environment. In particular, we make a comparison between the constant, periodic and conditional harvesting strategies. We find that for large initial populations, constant harvest is more beneficial to both the population and the maximum sustainable yield. However, periodic harvest has a short-term advantage when the initial population is low, and conditional harvest has the advantage of lowering the risk of depletion or extinction. Also, we investigate the periodic character under each strategy and show that periodic harvesting drives population cycles to be multiples (period-wise) of the harvesting period.     

The role of noise in a predator–prey model with Allee effect

The existence and implications of alternative stable states in ecological systems have been investigated extensively within deterministic models. However, it is known that natural systems are undeniably subject to random fluctuations, arising from either environmental variability or internal effects. Thus, in this paper, we study the role of noise on the pattern formation of a spatial predator–prey model with Allee effect. The obtained results show that the spatially extended system exhibits rich dynamic behavior. More specifically, the stationary pattern can be induced to be a stable target wave when the noise intensity is small. As the noise intensity is increased, patchy invasion emerges. These results indicate that the dynamic behavior of predator–prey models may be partly due to stochastic factors instead of deterministic factors, which may also help us to understand the effects arising from the undeniable susceptibility to random fluctuations of real ecosystems.     

Global behaviour of a predator–prey like model with piecewise constant arguments

The present study deals with the analysis of a predator–prey like model consisting of system of differential equations with piecewise constant arguments. A solution of the system with piecewise constant arguments leads to a system of difference equations which is examined to study boundedness, local and global asymptotic behaviour of the positive solutions. Using Schur–Cohn criterion and a Lyapunov function, we derive sufficient conditions under which the positive equilibrium point is local and global asymptotically stable. Moreover, we show numerically that periodic solutions arise as a consequence of Neimark-Sacker bifurcation of a limit cycle.     

Multiple periodic solutions of a delayed predator–prey model with non-monotonic functional response and stage structure

The paper studies a periodic and delayed predator–prey system with non-monotonic functional responses and stage structure. In the system, both the predator and prey are divided into immature individuals and mature individuals by two fixed ages. It is assumed that the immature predators cannot attack preys, and the case of the mature predators attacking the immature preys is also ignored. Based on Mawhin's coincidence degree, sufficient conditions are obtained for the existence of two positive periodic solutions of the system. An example is presented to illustrate the feasibility of the main results.