A mathematical model between keystone species: Bears, salmon and vegetation

: 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 di ﬀ ers signiﬁcantly due to the nutritional level di ﬀ erence between salmon and vegetation. We propose a stoichiometric predator-prey model that not only tracks the energy ﬂow 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 di ﬀ er under di ﬀ erent 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.


Introduction
A keystone species is one on which other species in the ecosystem largely depend. So much so that if removed there are drastic effects on the ecosystem [1]. Bears, salmon and vegetation are unique in an ecosystem because all three are keystone species dependent on one another [2]. In this ecosystem, salmon and vegetation are both food resources for bears, but the food quality differs significantly.
Salmon are much more nutrient-rich compared to vegetation, i.e., the nitrogen/phosphorous etc. concentration within salmon's bodies is much higher compared to vegetation. However, salmon must maintain homeostasis of element concentration within the body so that they must excrete excessive nutrients into the environment if the element concentration exceeds a threshold. On the other side, vegetation is less nutrient-rich but can absorb the remaining nutrient in the environment so the nutrient level within vegetation may differ significantly over time at an order of different magnitude [3][4][5].
Bears share common characteristics with salmon where bears are also nutrient-rich but must main-tain a homeostasis of element concentration within their bodies [3][4][5]. When consuming poor-quality food where the element concentration is low, the energy conversion from food consumption to the bear biomass cannot be maximized. In the bear-salmon-vegetation ecosystem, salmon is good-quality food, whereas vegetation is poor-quality food. Therefore, when studying the bear-salmon-vegetation ecosystem, it is important to devise a model that depicts not only the energy transfer but also the nutrient flow between trophic levels. One of the well-known models that study stoichiometric population dynamics is the one proposed by Loladze and Kuang [6]. The model builds on classical predator-prey models of one consumer and one food resource but incorporates a single limiting element phosphorus into modeling to distinguish food quantity and food quality. Following [6], stoichiometric models have been studied extensively, see [7][8][9][10][11][12][13] for example. More recently, Phan, Elser and Kuang extended the previous producer-grazer model framework by including multiple shared limiting elements in the modeling [14]. The formulated model excludes the non-smoothness in the original stoichiometric model but replicates qualitatively similar dynamics in a wide range of parameters [14]. We follow the modeling framework of the aforementioned studies but extend the model to include a specific nutrient-rich food resource, which agrees with the bear-salmon-vegetation ecosystem. In particular, we choose nitrogen as the limiting element of the system because nitrogen is vital for all salmon, bears and vegetation.
We organize the paper in the following. In Section 2, we formulate a stoichiometric model of one consumer and two food resources, where one food resource is of good quality but the other food resource is of poor quality. In Section 3, we analyze the existence and local stability of the steadystate solutions. In Section 4, we conduct numerical simulations that confirm the analytical results. Moreover, numerical simulations also reveal interesting dynamics when certain parameters vary in a range. We end this paper in Section 5 with conclusions and discussions.

Model formulation
In the bear-salmon-vegetation ecosystem, salmon and vegetation are both food resources for bears but bears do not consume certain organs of salmon. Rather, bears leave salmon corpses around vegetation so that the vegetation absorbs any remaining nitrogen that ultimately contributes to the growth of its population. Moreover, bears and salmon are much more nutrient-rich compared to vegetation. However, salmon and bears must maintain homeostasis in nitrogen concentration within their bodies and excrete excessive nitrogen into the environment. Vegetation is less nutrient-rich but can absorb the remaining nitrogen in the environment so the nitrogen concentration in vegetation may vary at an order of different magnitude.
We first formulate a model that includes the energy transfer among trophic levels when nitrogen is abundant among all three species. Denote S (t), V(t), B(t) by the salmon population density (measured in carbon per litre), the vegetation population density (measured in carbon per litre) and the bear population density (measured in carbon per litre) respectively. The model is where α is the growth rate of the salmon population, f (S ) is the consumption rate of salmon by the bear, γ is the intrinsic growth rate of the vegetation, K is the carrying capacity of the environment, g(V) represents the consumption rate of vegetation by the bear, e 1 and e 2 are the energy conversion rate from the consumption of the salmon and the vegetation respectively and λ is the natural death rate of the bear. In (2.1), represents the rate of contribution to the vegetation growth due to the element recycling from the nearby salmon corpses left by bears. Furthermore, the bears produce excrement after digesting the salmon, which provides another source for vegetation to absorb the nitrogen, which contributes positively to vegetation growth at a rate of ψ [2].
The above model (2.1) assumes a linear growth rate for the salmon and a logistic growth for the vegetation. Having spent most of their lives feeding and growing at sea, salmon return to spawn and die. Salmons return to their natal streams carrying marine-derived nutrients in their body tissues. In particular, adult salmon are rich in nitrogen [2]. Salmon do not compete for resources which implies it is safe to assume a linear growth rate of γ. Vegetation, on the other hand, competes for resources limited by solar energy and follows a logistic growth rate. However, the vegetation population is not only limited by light but is also regulated by nitrogen concentration.
Because of the homeostasis of bears and salmon, we can safely assume that the N: C (nitrogen: carbon) ratio within bears and salmon are θ 1 and θ 2 respectively. It follows that the remaining nitrogen in the environment is N − θ 1 B(t) − θ 2 S (t). Moreover, following [6], we assume that the N: C (nitrogen: carbon) ratio within the vegetation never falls below a level N V . Combining the assumptions, we obtain another carrying capacity of the vegetation limited by nitrogen as Together with the carrying capacity limited by solar energy, we obtain an improved carrying capacity for the vegetation population Following [6], we assume that all of the nitrogen within the system gets recycled immediately and there is no pool of free nitrogen in the environment.
It follows that the N: C ratio in vegetation is Vegetation serves as good-quality food for bears if and is poor-quality food otherwise. If vegetation is good-quality food, the biomass conversion from the vegetation to the bear follows the maximum energy intake e 2 . However, if vegetation is poor-quality food, the biomass conversion is reduced by a ratio Hence, biomass conversion efficiency is not a constant but depends on both energetic and nutrient limitations. Note that we can assume that salmon is a good-quality food source for bears since salmon is nutrient-rich. Following [6], we formulate the biomass conversion from vegetation consumption as Taking all the aforementioned evidence into consideration, we arrive at the model below where the parameters are listed in Table 1.
In the following analysis, we assume a linear functional response for both predation of the salmon and the vegetation, i.e., f (S ) = βS and g(V) = δV. We analyze a special case where the element recycling for the vegetation is minimum and can be ignored, i.e., = 0 and ψ = 0. In the following analysis, we restrict the positive invariant set Ω to be where the population densities are bounded by the total element concentration N in the ecosystem.

Steady state solutions
We first analyze steady-state solutions of (2.2) when = 0, ψ = 0. The steady-state solutions are determined by Here, (3.1) implies that S = 0 or B = α/β. If salmon is at extinction, i.e., S = 0, then (3.3) implies that If the bear population is at extinction, i.e., B = 0, then (3.2) simplifies to Therefore, we obtain an extinction equilibrium E 0 (0, 0, 0) and a vegetation-only equilibrium E 1 (0, min{K, N/N V }, 0) that always exist. If Therefore, a boundary equilibrium which is equivalent to where We solve for B and reject the larger root to obtain Therefore, a boundary equilibrium Finally, when (3.5) holds, (3.2) reduces to (3.10) Solving for V and rejecting the trivial solution leads tõ Therefore, a boundary equilibrium Next, we analyze the existence of steady-state solutions when the bear population exists, i.e.,B = α/β. SubstitutingB into (3.3) gives We further substituteS into (3.2) and obtain By solving for V and denoting the positive solution byV, we obtain Therefore, a positive equilibrium Solving for S gives (3.14) We substitute (3.14) into (3.2) and obtain Solving for V and denoting the positive solution byV giveŝ Solving for V and denoting the positive solution by V + leads to Therefore, a positive equilibrium

Stability of the steady-state solutions
Next, we analyze the local stability of the steady-state solutions. Direct calculations lead to the Jacobian matrix Then, the characteristic equation follows as |µI − J| = 0. By substituting each equilibrium into the characteristic equation, we are able to obtain the following theorems that state the stability results. The theorem below shows the stability of the bear-extinction equilibria.
which demonstrates that E 0 is unstable.
Similar calculations give the characteristic roots at E 1 which shows that E 1 is unstable.
The following theorem shows the stability results of the salmon-extinction equilibria E 2 , E 3 , E 4 .
Theorem 3.2. The salmon-extinction equilibrium E 2 is locally asymptotically stable if α/β < B * and is unstable if otherwise. The salmon-extinction equilibrium E 3 is locally asymptotically stable if α/β < B * * < γ/δ and is unstable if otherwise. The salmon-extinction equilibrium E 4 is locally asymptotically stable if α/β <B < γ/δ and is unstable if otherwise.
Finally, direct calculations give the characteristic equation at E 4 We can further simplify J 22 by using (3.10) to obtain J 22 = δB − γ. Moreover, J 33 < 0 is satisfied automatically when E 4 exists. Therefore, E 4 is locally asymptotically stable if α/β <B < γ/δ and is unstable otherwise.
Finally, the following theorem shows the local stabilities of the positive equilibria E 5 , E 6 , E 7 . Proof. Direct calculations give the characteristic equation at E 5 By using (3.12), we can simplify J 33 as J 33 = −e 2 δB < 0. Moreover, by using (3.13), we can simplify J 22 to The above analyses show that E 5 is locally asymptotically stable if J 31 > 0 or equivalently e 1 βθ 1 > e 2 δθ 2 .

By the Routh-Hurwitz criterion, E 7 is locally asymptotically stable if
Here it is obvious that −(J 23 J 32 + J 13 J 31 ) > 0. Next, by substituting (3.14) and (3.17) into J 22 , we can simplify Here J 31 J 22 − J 21 J 32 < 0 can be simplified to 0 < −γθ 2 N V e 2 δV 2 by substituting (3.14) and (3.17), which leads to the contradiction. Therefore, E 7 is always unstable when exists. Thus, completes the proof.

Simulation of (2.2) with a linear predation
Now we explore the dynamics of (2.2) numerically. For the numerical simulation, we will be assuming that the salmon-bear-vegetation ecological system lies within the riparian forests of Alaska. More specifically, we will be looking at a 100km 2 region of Lynx Creek, a tributary of the Wood River Lakes system in the Bristol Bay region of southwestern Alaska, USA (59 • 29 N, 158 • 55 W) [2]. This is a well-documented region where we will be able to get accurate values for our parameters. Based on data from [2], we explore the following scenarios by using biologically realistic parameters. Figure 1 shows that the salmon-extinction boundary equilibrium is locally asymptotically stable. In particular, 1(a) demonstrates that E 2 is locally asymptotically stable, whereas both E 3 , E 4 do not exist under the parameter set. Figure 1(b) indicates that E 3 is locally asymptotically stable, whereas both E 2 , E 4 do not exist under the parameter set. Moreover, Figure 1(c) shows that E 4 is locally asymptotically stable, whereas both E 2 , E 3 do not exist under the parameter set. Via our extensive numerical experiments, we find that E 2 cannot coexist with either/both E 3 , E 4 . Figure 2 demonstrates that all salmon, vegetation and bear populations coexist together. Figure 2(a) shows that the positive equilibrium E 5 is locally asymptotically stable, whereas E 6 does not exist under the parameter set. Figure 2(b) shows that the positive equilibrium E 6 is locally asymptotically stable, whereas E 5 does not exist under the parameter set.  Next, we explore how the vegetation growth rate impacts the long-term dynamics of (2.2). We generate Figure 3 by using the same set of parameter as Figure 2(b) but varying γ. Figure 3(a) demonstrates that a small vegetation growth rate drives the vegetation population to extinction, whereas salmon coexist with bears, but the populations oscillate periodically. On the other hand, the salmon population goes to extinction, but the vegetation coexists with bears at a steady state when γ is large. The results are reasonable biologically because a small vegetation growth rate cannot sustain the persistence of the vegetation population. However, a large vegetation growth rate facilitates the nitrogen recycling of the vegetation, which supports the bear population via consumption but drives the salmon population to extinction due to the scarcity of nitrogen.
Finally, Figure 4(a) demonstrates another scenario where the salmon, bear and vegetation coexist but the populations oscillate periodically. However, Figure 4(b) shows that a larger vegetation growth rate may stabilize the oscillating populations to a steady state.

Simulation of (2.2) with a logistic growth of the salmon population and Holling type II predation
Next, we explore the extended model of (2.2) where the salmon population follows a logistic growth We consider that the predation of the salmon and the vegetation now follow the Holling type II functional response [15] where h 1 , h 2 represent the half-saturating constants respectively. Figure 5(a) demonstrates the bifurcation diagram of the salmon population with respect to the salmon carrying capacity K 1 . Moreover, Figure 5(b) demonstrates the bifurcation diagram of the vegetation population with respect to K 1 . The bifurcation diagram of the bear population with respect to K 1 is similar to Figure 5(a) and is thus omitted. Figure 5 shows that the salmon, vegetation and bear populations coexist at a steady state if the salmon carrying capacity is relatively small. However, when K 1 increases and passes 4.1, the vegetation population goes to extinction, whereas the salmon population coexists with the bear population at a steady state. Moreover, when K 1 further increases and passes 5.7, the vegetation population remains at extinction, but the salmon and bear populations oscillate periodically. The results are biologically reasonable because if the salmon carrying capacity is relatively small, the salmon, vegetation and bear compete for the nutrients but coexist due to the scarcity of the salmon population. However, if the salmon carrying capacity is at the intermediate range, the competition for nutrients drives the vegetation population to extinction because the salmon and bear are more nutrient-rich. Finally, a large salmon carrying capacity destabilizes the coexistence steady-state of the salmon and bear due to the enrichment of salmon. Next, to compare the dynamics of the salmon-vegetation-bear model with the grazer-consumer model in [6], we run the simulation by using the same set of parameters as in [6], except the salmon carrying capacity K 1 and the half-saturating constant of the salmon consumption h 1 . Figure 6 shows the bifurcation diagram of the salmon population/vegetation population/bear population with respect to the vegetation carrying capacity K respectively. If K is relatively small, Figure 6 shows that the salmon, bear and vegetation populations coexist at a steady state. However, when K increases and passes 0.24, the salmon population goes to extinction, whereas the vegetation and bear populations coexist at a steady state. If K further increases and passes 0.55, the vegetation and bear coexist but the populations oscillate periodically. If K further increases and passes 0.97, the vegetation and bear populations coexist but return to the steady state status. Finally, if K becomes relatively large and passes 1.5, the salmon population is no longer extinct but coexists with the vegetation and bear populations at a steady state.
Overall, when K is within the intermediate range, i.e., the salmon population is at extinction, Figures 6(a), 6(b) and 6(c) demonstrate similar dynamics with the producer-grazer model in [6]. The vegetation-bear coexistent steady state loses stability due to energy enrichment. However, a larger car-rying capacity limited by the light energy drives the oscillating vegetation and bear populations to a steady state because of the competition for nutrients. Different from the results in [6], the bear population does not go to extinction but coexists with either/both salmon or/and vegetation populations when K is relatively small or large. The results suggest that a low or high carrying capacity limited by light energy facilitates the persistence of all salmon, vegetation and bear populations. Moreover, an intermediate carrying capacity limited by light energy drives the salmon population to extinction.

Conclusions and discussion
In this paper, we study an ecosystem of bears, salmon and vegetation, where bears consume salmon and vegetation for survival. Because salmon return to their natal streams to spawn and carry marinederived nutrients, it is important to characterize the nutritional level difference between salmon, bears, and vegetation. In general, salmon and bears are more nutrient-rich compared to vegetation but must maintain homeostasis of the nutritional level within their bodies. On the other side, vegetation is less nutrient-rich but can recycle the remaining nutrient in the environment and the nutritional level within the body may differ significantly.
We propose a stoichiometric predator-prey model that tracks both the energy flow and the nutrient recycling from one trophic level to another. Analytical results show that boundary equilibria E 0 , E 1 where bears are extinct exist but are always unstable. Moreover, boundary equilibria E 2 , E 3 , E 4 where salmon are extinct but bears persist may exist and remain locally asymptotically stable if certain conditions are satisfied. Positive equilibria E 5 , E 6 , E 7 where salmon, bears and vegetation coexist may exist if certain conditions are satisfied. Analyses show that E 5 and E 6 may remain locally asymptotically stable under certain conditions, but E 7 is always unstable.
Numerical simulations demonstrate that a small vegetation growth rate may drive the vegetation population to extinction where the salmon population and the bear population coexist in the periodic setting. Moreover, a large vegetation growth rate may drive the salmon population to extinction but the vegetation coexists with the bears at a steady state. Alternatively, the salmon, bears and vegetation populations may coexist periodically. In this scenario, a large vegetation growth rate may stabilize the system and drive the salmon, bear and vegetation to coexist at a steady state.
In this paper, for the analytical analysis, we adopt the linear functional response for the predation of the salmon and the predation of the vegetation in (2.2) for simplicity. However, a linear functional response has its limitation and is more suitable for an ecosystem with a sparse population density. A Holling type II functional response has a saturating effect when population density becomes large and therefore may be suitable for a broader regime [15]. Figures 5 and 6 in the simulation also confirm that rich dynamics may occur if the Holling type II functional response is adopted. Because of the nonlinearity of the Holling type II functional response, we expect that the analytical analysis becomes more challenging but on the other hand may deepen our understanding of the ecosystem of keystone species, which leaves as future works.

Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.