We study an ecosystem of three keystone species: salmon, bears and vegetation. Bears consume salmon and vegetation for energy and nutrient intake, but the food quality differs significantly due to the nutritional level difference between salmon and vegetation. We propose a stoichiometric predator-prey model that not only tracks the energy flow from one trophic level to another but also nutrient recycling in the system. Analytical results show that bears may coexist with salmon and vegetation at a steady state, but the abundance of salmon may differ under different regimes. Numerical simulations reveal that a smaller vegetation growth rate may drive the vegetation population to extinction, whereas a large vegetation growth rate may drive the salmon population to extinction. Moreover, a large vegetation growth rate may stabilize the system where the bear, salmon and vegetation populations oscillate periodically.
Citation: Christopher Middlebrook, Xiaoying Wang. A mathematical model between keystone species: Bears, salmon and vegetation[J]. Mathematical Biosciences and Engineering, 2023, 20(9): 16628-16647. doi: 10.3934/mbe.2023740
We study an ecosystem of three keystone species: salmon, bears and vegetation. Bears consume salmon and vegetation for energy and nutrient intake, but the food quality differs significantly due to the nutritional level difference between salmon and vegetation. We propose a stoichiometric predator-prey model that not only tracks the energy flow from one trophic level to another but also nutrient recycling in the system. Analytical results show that bears may coexist with salmon and vegetation at a steady state, but the abundance of salmon may differ under different regimes. Numerical simulations reveal that a smaller vegetation growth rate may drive the vegetation population to extinction, whereas a large vegetation growth rate may drive the salmon population to extinction. Moreover, a large vegetation growth rate may stabilize the system where the bear, salmon and vegetation populations oscillate periodically.
[1] | R. T. Paine, A conversation on refining the concept of keystone species, Conserv. Biol., 9 (1995), 962–964. https://doi.org/10.1046/j.1523-1739.1995.09040962.x doi: 10.1046/j.1523-1739.1995.09040962.x |
[2] | J. M. Helfield, R. J. Naiman, Keystone interactions: salmon and bear in riparian forests of Alaska, Ecosystems, 9 (2006), 167–180. https://doi.org/10.1007/s10021-004-0063-5 doi: 10.1007/s10021-004-0063-5 |
[3] | T. Andersen, J. J. Elser, D. O. Hessen, Stoichiometry and population dynamics, Ecol. Lett., 7 (2004), 884–900. https://doi.org/10.1111/j.1461-0248.2004.00646.x doi: 10.1111/j.1461-0248.2004.00646.x |
[4] | J. J. Elser, W. F. Fagan, R. F. Denno, D. R. Dobberfuhl, A. Folarin, A. Huberty, et al., Nutritional constraints in terrestrial and freshwater food webs, Nature, 408 (2000), 578–580. https://doi.org/10.1038/35046058 doi: 10.1038/35046058 |
[5] | E. B. Muller, R. M. Nisbet, S. A. Kooijman, J. J. Elser, E. McCauley, Stoichiometric food quality and herbivore dynamics, Ecol. Lett., 4 (2001), 519–529. https://doi.org/10.1046/j.1461-0248.2001.00240.x doi: 10.1046/j.1461-0248.2001.00240.x |
[6] | I. Loladze, Y. Kuang, Stoichiometry in producer-grazer systems: linking energy flow with element cycling, Bull. Math. Biol., 62 (2000), 1137–1162. https://doi.org/10.1006/bulm.2000.0201 doi: 10.1006/bulm.2000.0201 |
[7] | M. Chen, M. Fan, Y. Kuang, Global dynamics in a stoichiometric food chain model with two limiting nutrients, Math. Biosci., 289 (2017), 9–19. https://doi.org/10.1016/j.mbs.2017.04.004 doi: 10.1016/j.mbs.2017.04.004 |
[8] | X. Li, H. Wang, Y. Kuang, Global analysis of a stoichiometric producer-grazer model with Holling-type functional responses, J. Math. Biol., 63 (2011), 901–932. https://doi.org/10.1007/s00285-010-0392-2 doi: 10.1007/s00285-010-0392-2 |
[9] | I. Loladze, Y. Kuang, J. J. Elser, W. F. Fagan, Competition and stoichiometry: coexistence of two predators on one prey, Theor. Popul Biol., 65 (2004), 1–15. https://doi.org/10.1016/S0040-5809(03)00105-9 doi: 10.1016/S0040-5809(03)00105-9 |
[10] | A. Peace, H. Wang, Y. Kuang, Dynamics of a producer-grazer model incorporating the effects of excess food nutrient content on grazer's growth, Bull. Math. Biol., 76 (2014), 2175–2197. https://doi.org/10.1007/s11538-014-0006-z doi: 10.1007/s11538-014-0006-z |
[11] | H. Stech, B. Peckham, J. Pastor, Enrichment in a general class of stoichiometric producer-consumer population growth models, Theor. Popul. Biol., 81 (2012), 210–222. https://doi.org/10.1016/j.tpb.2012.01.003 doi: 10.1016/j.tpb.2012.01.003 |
[12] | H. Wang, Y. Kuang, I. Loladze, Dynamics of a mechanistically derived stoichiometric producer-grazer model, J. Biol. Dyn., 2 (2008), 286–296. https://doi.org/10.1080/17513750701769881 doi: 10.1080/17513750701769881 |
[13] | T. Xie, X. Yang, X. Li, H. Wang, Complete global and bifurcation analysis of a stoichiometric predator-prey model, J. Dyn. Differ. Equations, 30 (2018), 447–472. https://doi.org/10.1007/s10884-016-9551-5 doi: 10.1007/s10884-016-9551-5 |
[14] | T. Phan, J. J. Elser, Y. Kuang, Rich dynamics of a general producer-grazer interaction model under shared multiple resource limitations, Appl. Sci., 13 (2023), 4150. https://doi.org/10.3390/app13074150 doi: 10.3390/app13074150 |
[15] | C. S. Holling, The functional response of predators to prey density and its role in mimicry and population regulation, The Memoirs of the Entomological Society of Canada, 97 (1965), 5–60. https://doi.org/10.4039/entm9745fv doi: 10.4039/entm9745fv |