# American Institute of Mathematical Sciences

January  2021, 26(1): 107-120. doi: 10.3934/dcdsb.2020264

## Dynamics of a discrete-time stoichiometric optimal foraging model

 1 School of Science, Dalian Maritime University, 1 Linghai Road, Dalian, Liaoning, 116026, China 2 Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, Alberta T6G 2G1, Canada

* Corresponding Author (hao8@ualberta.ca)

Received  April 2020 Revised  June 2020 Published  August 2020

Fund Project: Partially supported by NSFC-11801052, NSFLP-2019-ZD-1056, NSERC RGPIN-2020-03911 and NSERC RGPAS-2020-00090

In this paper, we discretize and analyze a stoichiometric optimal foraging model where the grazer's feeding effort depends on the producer's nutrient quality. We systematically make comparisons of the dynamical behaviors between the discrete-time model and the continuous-time model to study the robustness of model predictions to time discretization. When the maximum growth rate of producer is low, both model types admit similar dynamics including bistability and deterministic extinction of the grazer caused by low nutrient quality of the producer. Especially, the grazer is benefited from optimal foraging similarly in both discrete-time and continuous-time models. When the maximum growth rate of producer is high, dynamics of the discrete-time model are more complex including chaos. A phenomenal observation is that under extremely high light intensities, the grazer in the continuous-time model tends to perish due to poor food quality, however, the grazer in the discrete-time model persists in regular or irregular oscillatory ways. This significant difference indicates the necessity of studying discrete-time models which naturally include species' generations and are thus more popular in theoretical biology. Finally, we discuss how the shape of the quality-based feeding function regulates the beneficial or restraint effect of optimal foraging on the grazer population.

Citation: Ming Chen, Hao Wang. Dynamics of a discrete-time stoichiometric optimal foraging model. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 107-120. doi: 10.3934/dcdsb.2020264
##### References:

show all references

##### References:
Attractor of discrete-time optimal foraging model (4) in phase plane for different light intensities in two cases. Panels of ($a_i$) describe the case when the producer's growth rate is low $b = 1.2$, $i = 1,2,3$. Panels of ($b_i$) describe the case when the producer's growth rate is high $b = 3$, $i = 1,2,3$. The red dashed curves are defined by $F(x,y) = 0$, which denote the producer nullclines. The blue dotted curves are defined by $G(x,y) = 0$, which denote the grazer nullclines. Solid bullets denote stable equilibria while circles represent unstable equilibria
Bifurcation diagram of the population densities with respect to $K$ (light intensity) for the discrete-time model (4) ($a_i$), i = 1, 2, and the continuous-time model (1) ($b_i$), i = 1, 2. Shaded regions with $+$ represent the parameter regions of the optimal foraging behaviors benefiting the grazers. All parameters are provided in Table 1 with $b = 1.2$
Solution curves for system (4) and (1). ($a_i$) and ($b_i$) denote the dynamics of (4) and (1) with increasing $K$, respectively. Producer and grazer's densities ($\mathrm{mg\; C/L}$) are plotted by dashed and solid lines, respectively. All parameters are provided in Table 1 with $b = 1.2$
Solution curves for system (4) and (1). ($a_i$) and ($b_i$) denote the dynamics of (4) and (1) with increasing $K$, respectively. Producer and grazer densities ($\mathrm{mg\; C/L}$) are described by dashed and solid lines, respectively. All parameters are provided in Table 1 with $b = 3$
The bifurcation curves with respect to $K$ for the discrete-time model ($a_i$), i = 1, 2, and continuous-time model ($b_i$), i = 1, 2. All parameters are provided in Table 1 with $b = 3$
Spectrum of the maximum Lyapunov exponent (MLE) with respect to $K$ for the discrete-time model. All parameters are provided in Table 1 with $b = 3$
A two-parameter bifurcation diagram for varying light level $K$ and varying maximal growth rate of producer $b$ for the discrete-time model (a) and continuous-time model (b). All other parameter values are listed in Table 1 and the initial point is $x(0) = 0.2\; \mathrm{mg C L}^{-1}$ and $y(0) = 0.2\; \mathrm{mg C L}^{-1}$. Discrete-time model (4) exhibits periodic oscillations in blue region and chaotic behaviors in red region. Outside these regions, model (4) has stable equilibria
Bifurcation diagram of the grazer densities with respect to $K$ (light intensity) for the discrete-time model (4) ($a_i$), i = 1, 2, and the continuous-time model (1) ($b_i$), i = 1, 2. Specially, ($a_1$) and ($b_1$) denote the case with the low growth rate of producer ($b = 1.2$); ($a_2$) and ($b_2$) denote the case with the high growth rate of producer ($b = 3$). Light ($+$) and dark ($-$) shaded regions represent the parameter regions of the optimal foraging behaviors benefiting and restraining the grazers, respectively. All parameters are provided in Table 1 except the parameter $a_1 = 3.5$
Parameters of model (4) with default values and units
 Par. Description Value Unit $P_T$ Total phosphorus $0.02$ $\mathrm{mg P L}^{-1}$ $K$ Producer carrying capacity determined by light $0-3.5$ $\mathrm{mg C L}^{-1}$ $b$ Maximal growth rate of the producer $1.2$ or $3$ $\mathrm{day}^{-1}$ $\delta$ Grazer loss rate $0.12$ $\mathrm{day}^{-1}$ $\theta$ Grazer constant $\mathrm{P:C}$ 0.03 $\mathrm{mgP/mgC}$ $q$ Producer minimal $\mathrm{P:C}$ $0.0038$ $\mathrm{mgP/mgC}$ $e$ Maximal production efficiency in carbon terms for grazer $0.8$ $\alpha$ Phosphorus half saturation constant of the producer $0.008$ $\mathrm{mg C L}^{-1}$ $\mu$ Water cleared/mg C invested to generate filtering energy $700$ $\mathrm{L/mg C}$ $\tau$ Handling time (-inverse of max feeding rate) $1.23$ $\mathrm{day}$ $\xi(Q)$ Feeding cost, function for optimal foraging model $a_0=0.01$, $a_1=5.17$ $\xi(Q)=\min\{a_0,a_1Q^2+a_2Q+a_3\}$ $a_2=-0.31$, $a_3=0.007$
 Par. Description Value Unit $P_T$ Total phosphorus $0.02$ $\mathrm{mg P L}^{-1}$ $K$ Producer carrying capacity determined by light $0-3.5$ $\mathrm{mg C L}^{-1}$ $b$ Maximal growth rate of the producer $1.2$ or $3$ $\mathrm{day}^{-1}$ $\delta$ Grazer loss rate $0.12$ $\mathrm{day}^{-1}$ $\theta$ Grazer constant $\mathrm{P:C}$ 0.03 $\mathrm{mgP/mgC}$ $q$ Producer minimal $\mathrm{P:C}$ $0.0038$ $\mathrm{mgP/mgC}$ $e$ Maximal production efficiency in carbon terms for grazer $0.8$ $\alpha$ Phosphorus half saturation constant of the producer $0.008$ $\mathrm{mg C L}^{-1}$ $\mu$ Water cleared/mg C invested to generate filtering energy $700$ $\mathrm{L/mg C}$ $\tau$ Handling time (-inverse of max feeding rate) $1.23$ $\mathrm{day}$ $\xi(Q)$ Feeding cost, function for optimal foraging model $a_0=0.01$, $a_1=5.17$ $\xi(Q)=\min\{a_0,a_1Q^2+a_2Q+a_3\}$ $a_2=-0.31$, $a_3=0.007$
 [1] Veena Goswami, Gopinath Panda. Optimal customer behavior in observable and unobservable discrete-time queues. Journal of Industrial & Management Optimization, 2021, 17 (1) : 299-316. doi: 10.3934/jimo.2019112 [2] Cuicui Li, Lin Zhou, Zhidong Teng, Buyu Wen. The threshold dynamics of a discrete-time echinococcosis transmission model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020339 [3] Angelica Pachon, Federico Polito, Costantino Ricciuti. On discrete-time semi-Markov processes. Discrete & Continuous Dynamical Systems - B, 2021, 26 (3) : 1499-1529. doi: 10.3934/dcdsb.2020170 [4] Peter Giesl, Zachary Langhorne, Carlos Argáez, Sigurdur Hafstein. Computing complete Lyapunov functions for discrete-time dynamical systems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 299-336. doi: 10.3934/dcdsb.2020331 [5] Mahir Demir, Suzanne Lenhart. A spatial food chain model for the Black Sea Anchovy, and its optimal fishery. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 155-171. doi: 10.3934/dcdsb.2020373 [6] Tien-Yu Lin, Bhaba R. Sarker, Chien-Jui Lin. An optimal setup cost reduction and lot size for economic production quantity model with imperfect quality and quantity discounts. Journal of Industrial & Management Optimization, 2021, 17 (1) : 467-484. doi: 10.3934/jimo.2020043 [7] Haixiang Yao, Ping Chen, Miao Zhang, Xun Li. Dynamic discrete-time portfolio selection for defined contribution pension funds with inflation risk. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020166 [8] Guangjun Shen, Xueying Wu, Xiuwei Yin. Stabilization of stochastic differential equations driven by G-Lévy process with discrete-time feedback control. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 755-774. doi: 10.3934/dcdsb.2020133 [9] Lars Grüne, Matthias A. Müller, Christopher M. Kellett, Steven R. Weller. Strict dissipativity for discrete time discounted optimal control problems. Mathematical Control & Related Fields, 2020  doi: 10.3934/mcrf.2020046 [10] Xu Zhang, Chuang Zheng, Enrique Zuazua. Time discrete wave equations: Boundary observability and control. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 571-604. doi: 10.3934/dcds.2009.23.571 [11] Reza Lotfi, Zahra Yadegari, Seyed Hossein Hosseini, Amir Hossein Khameneh, Erfan Babaee Tirkolaee, Gerhard-Wilhelm Weber. A robust time-cost-quality-energy-environment trade-off with resource-constrained in project management: A case study for a bridge construction project. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020158 [12] Chao Xing, Jiaojiao Pan, Hong Luo. Stability and dynamic transition of a toxin-producing phytoplankton-zooplankton model with additional food. Communications on Pure & Applied Analysis, 2021, 20 (1) : 427-448. doi: 10.3934/cpaa.2020275 [13] Claudio Arancibia-Ibarra, José Flores, Michael Bode, Graeme Pettet, Peter van Heijster. A modified May–Holling–Tanner predator-prey model with multiple Allee effects on the prey and an alternative food source for the predator. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 943-962. doi: 10.3934/dcdsb.2020148 [14] Simone Göttlich, Elisa Iacomini, Thomas Jung. Properties of the LWR model with time delay. Networks & Heterogeneous Media, 2020  doi: 10.3934/nhm.2020032 [15] Stefan Siegmund, Petr Stehlík. Time scale-induced asynchronous discrete dynamical systems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1011-1029. doi: 10.3934/dcdsb.2020151 [16] Ting Liu, Guo-Bao Zhang. Global stability of traveling waves for a spatially discrete diffusion system with time delay. Electronic Research Archive, , () : -. doi: 10.3934/era.2021003 [17] Theresa Lange, Wilhelm Stannat. Mean field limit of ensemble square root filters - discrete and continuous time. Foundations of Data Science, 2021  doi: 10.3934/fods.2021003 [18] Guihong Fan, Gail S. K. Wolkowicz. Chaotic dynamics in a simple predator-prey model with discrete delay. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 191-216. doi: 10.3934/dcdsb.2020263 [19] Ran Zhang, Shengqiang Liu. On the asymptotic behaviour of traveling wave solution for a discrete diffusive epidemic model. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1197-1204. doi: 10.3934/dcdsb.2020159 [20] Matúš Tibenský, Angela Handlovičová. Convergence analysis of the discrete duality finite volume scheme for the regularised Heston model. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 1181-1195. doi: 10.3934/dcdss.2020226

2019 Impact Factor: 1.27

## Tools

Article outline

Figures and Tables