Loading [Contrib]/a11y/accessibility-menu.js
Research article Special Issues

Mathematical analysis of a chemostat system modeling the competition for light and inorganic carbon with internal storage

  • This paper investigates a mathematical model of competition between two species for inorganic carbon and light in a well-mixed water column. The population growth of the species depends on the consumption of two substitutable forms of inorganic carbon, "CO2" (dissolved CO2 and carbonic acid) and "CARB" (bicarbonate and carbonate ions), which are stored internally. Besides, uptake rates also includes self-shading by the phytoplankton population, that is, an increase in population density will reduce light available for photosynthesis, and thereby reducing further carbon assimilation and population growth. We also incorporate the fact that carbon is lost by respiration, and the respiration rate is assumed to be proportional to the size of the transient carbon pool. Then we study the extinction and persistence of a single-species system. Finally, we show that coexistence of the two-species system is possible, depending on parameter values, and both persistence of one population.

    Citation: Fu-Yuan Tsai, Feng-BinWang. Mathematical analysis of a chemostat system modeling the competition for light and inorganic carbon with internal storage[J]. Mathematical Biosciences and Engineering, 2019, 16(1): 205-221. doi: 10.3934/mbe.2019011

    Related Papers:

    [1] Ling Xue, Sitong Chen, Xinmiao Rong . Dynamics of competition model between two plants based on stoichiometry. Mathematical Biosciences and Engineering, 2023, 20(10): 18888-18915. doi: 10.3934/mbe.2023836
    [2] Xiong Li, Hao Wang . A stoichiometrically derived algal growth model and its global analysis. Mathematical Biosciences and Engineering, 2010, 7(4): 825-836. doi: 10.3934/mbe.2010.7.825
    [3] Zhun Han, Hal L. Smith . Bacteriophage-resistant and bacteriophage-sensitive bacteria in a chemostat. Mathematical Biosciences and Engineering, 2012, 9(4): 737-765. doi: 10.3934/mbe.2012.9.737
    [4] Yun Kang, Sourav Kumar Sasmal, Amiya Ranjan Bhowmick, Joydev Chattopadhyay . Dynamics of a predator-prey system with prey subject to Allee effects and disease. Mathematical Biosciences and Engineering, 2014, 11(4): 877-918. doi: 10.3934/mbe.2014.11.877
    [5] Junjing Xiong, Xiong Li, Hao Wang . The survival analysis of a stochastic Lotka-Volterra competition model with a coexistence equilibrium. Mathematical Biosciences and Engineering, 2019, 16(4): 2717-2737. doi: 10.3934/mbe.2019135
    [6] Zhenyao Sun, Da Song, Meng Fan . Dynamics of a stoichiometric phytoplankton-zooplankton model with season-driven light intensity. Mathematical Biosciences and Engineering, 2024, 21(8): 6870-6897. doi: 10.3934/mbe.2024301
    [7] Yuanshi Wang, Hong Wu . Transition of interaction outcomes in a facilitation-competition system of two species. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1463-1475. doi: 10.3934/mbe.2017076
    [8] Nahla Abdellatif, Radhouane Fekih-Salem, Tewfik Sari . Competition for a single resource and coexistence of several species in the chemostat. Mathematical Biosciences and Engineering, 2016, 13(4): 631-652. doi: 10.3934/mbe.2016012
    [9] Yang Kuang, Kaifa Wang . Coexistence and extinction in a data-based ratio-dependent model of an insect community. Mathematical Biosciences and Engineering, 2020, 17(4): 3274-3293. doi: 10.3934/mbe.2020187
    [10] Xinyan Chen, Zhaohui Jiang, Qile Tai, Chunshan Shen, Yuan Rao, Wu Zhang . Construction of a photosynthetic rate prediction model for greenhouse strawberries with distributed regulation of light environment. Mathematical Biosciences and Engineering, 2022, 19(12): 12774-12791. doi: 10.3934/mbe.2022596
  • This paper investigates a mathematical model of competition between two species for inorganic carbon and light in a well-mixed water column. The population growth of the species depends on the consumption of two substitutable forms of inorganic carbon, "CO2" (dissolved CO2 and carbonic acid) and "CARB" (bicarbonate and carbonate ions), which are stored internally. Besides, uptake rates also includes self-shading by the phytoplankton population, that is, an increase in population density will reduce light available for photosynthesis, and thereby reducing further carbon assimilation and population growth. We also incorporate the fact that carbon is lost by respiration, and the respiration rate is assumed to be proportional to the size of the transient carbon pool. Then we study the extinction and persistence of a single-species system. Finally, we show that coexistence of the two-species system is possible, depending on parameter values, and both persistence of one population.


    It has been known that organic nutrients (e.g., nitrogen and phosphorus), light, and inorganic carbon are the important factors that affect the growth of phytoplankton. However, previous competition theory only focused on the interaction between the species and nutrients/light (see, e.g., [7,9]), and neglected the role of inorganic carbon. This is probably due to the complexities including the biochemistry of carbon acquisition by phytoplankton and the geochemistry of inorganic carbon in the ecosystem [10,20]. In the Supplementary Information of [20], the authors proposed a system of ODEs modeling the competition of the species for inorganic carbon and light in a well-mixed water column. Dissolved CO2 and carbonic acid are regarded as one resource (denoted as "CO2"), and bicarbonate and carbonate ions are regarded as another (denoted as "CARB"). The resources "CO2" and "CARB" are stored internally, and they are substitutable in their effects on algal growth [17,20]. On the other hand, uptake rates also includes self-shading by the phytoplankton population, namely, an increase in population density will reduce light available for photosynthesis, and thereby suppressing further carbon assimilation and population growth [20].

    It was known that pH and alkalinity are two main factors in the modeling of inorganic carbon [20]. The consumption terms for "CO2" and "CARB" used in [20] include computations of feedbacks that arise from changes in pH and alkalinity during algal growth. In the recent work [17], the authors ignore these latter feedbacks and assume that the parameters in the system are constants, simplifying the complex processes of "CO2" and "CARB" involved. Incorporating the simplifications used in [17], we modify the model presented in the Supplementary Information of [20] and we shall investigate the following chemostat-type model with internal storage:

    $ \begin{equation}\label{Droop-Co2-light-ori} \begin{cases} \frac{d R}{dt} = (R^{(0)}-R)D-f_{R1}(R, Q_{1})g_1(u_1, u_2)u_{1}-f_{R2}(R, Q_{2})g_2(u_1, u_2)u_{2}\\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ +\gamma_1(Q_1)u_1+\gamma_2(Q_2)u_2-\omega_rR+\omega_sS, \\ \frac{dS}{dt} = (S^{(0)}-S)D-f_{S1}(S, Q_{1})g_1(u_1, u_2)u_{1}-f_{S2}(S, Q_{2})g_2(u_1, u_2)u_{2}+\omega_rR-\omega_sS, \\ \frac{dQ_{1}}{dt} = f_{R1}(R, Q_{1})g_1(u_1, u_2)+f_{S1}(S, Q_{1})g_1(u_1, u_2)-\mu_{1}(Q_{1})Q_{1}-\gamma_1(Q_1), \\ \frac{du_{1}}{dt} = \left[\mu_{1}(Q_{1})-D\right]u_{1}, \\ \frac{dQ_{2}}{dt} = f_{R2}(R, Q_{2})g_2(u_1, u_2)+f_{S2}(S, Q_{2})g_2(u_1, u_2)-\mu_{2}(Q_{2})Q_{2}-\gamma_2(Q_2), \\ \frac{du_{2}}{dt} = \left[\mu_{2}(Q_{2})-D\right]u_{2}, \\ R(0)\ge0, \ S(0)\ge0, \ u_{i}(0)\ge 0, \ Q_{i}(0)\ge Q_{\min, i}, \ i = 1, 2. \end{cases} \end{equation} $ (1.1)

    Here $R(t)$ and $S(t)$ denote the concentrations of "CO2" and "CARB" in the chemostat at time $t$, respectively. $u_{i} (t)$ denotes the concentration of species $i$ at time $t$. $Q_{i}$ represents the amount of cell quota of resources $R$ and $S$ per individual of species $i$ at time $t$. $\mu_{i}(Q_{i})$ is the growth rates of species $i$ as a function of cell quota $Q_{i}$. $f_{Ri}(R, Q_{i})$ ($f_{Si}(S, Q_{i})$) is the per capita uptake rate of species $i$ as a function of resource concentration $R$ ($S$) and cell quota $Q_{i}$. $D$ is the dilution rate of the chemostat. Each nutrient is supplied at the rate $D$, and both input concentrations are $R^{(0)}$ and $S^{(0)}$ respectively. $Q_{\min, i}$ denotes threshold cell quota below which no growth of species $i$ occurs. $\gamma_i(Q_i)$ represents the respiration rate of species $i$ as a function of cell quota $Q_{i}$. $g_i(u_1, u_2)$ stands for the photosynthetic rate of the species $i$ as a function of $u_{1} (t)$ and $u_{2} (t)$. Following the ideas of model simplifications in [17], we also assume that carbonic acid loses a proton to become bicarbonate at the rate $\omega_r$, and the rate of the reverse reaction is denoted by $\omega_s$.

    One type of the photosynthetic rate of the species $i$, $g_i(u_1, u_2)$, takes

    $ \begin{equation}\label{light-g-3} g_i(u_1, u_2) = \frac{1}{z_m}\int_0^{z_m}\frac{m_i\hat{I}(z)}{a_i+\hat{I}(z)}dz, \end{equation} $ (1.2)

    where $m_i$ and $a_i$ are the maximal growth rate and half saturation constant of species $i$, respectively, and

    $ \begin{equation}\label{light-g-4} \hat{I}(z) = I_{\rm in} \exp\left(-k_0 z-k_1 z u _1(t)-k_2 z u _2(t)\right). \end{equation} $ (1.3)

    Here we have assumed that the light intensity at each depth is described by Lambert-Beer law [11,14], which states that the amount of light absorbed is proportional to the light intensity ($I_{\rm in}$), but decreases with the depth in the water column ($z$), the background turbidity of the water itself ($k_0$), the specific light attenuation coefficients of the competing species ($k_i$), and the population densities of the species ($u_i$). Assume that $z_m$ is the total depth of the water column. Substituting (1.3) into (1.2), then the photosynthetic rate $g_i(u_1, u_2)$ becomes

    $ \begin{equation*}\label{light-g-5} g_i(u_1, u_2) = \frac{m_i}{(k_0+k_1u_1+k_2u_2)z_m}\ln\left(\frac{a_i+I_{\rm in}}{a_i+\hat{I}(z_m)}\right) = \frac{m_i}{\ln\left(\frac{I_{\rm in}}{\hat{I}(z_m)}\right)}\ln\left(\frac{a_i+I_{\rm in}}{a_i+\hat{I}(z_m)}\right), \end{equation*} $

    where

    $ \begin{equation*}\label{light-g-6} \hat{I}(z_m) = I_{\rm in} \exp\left(-k_0 z_m-k_1 z_m u _1(t)-k_2 z_m u _2(t)\right). \end{equation*} $

    We also note that the other type of the photosynthetic rate $g_i(u_1, u_2)$ takes the form

    $ \begin{equation}\label{light-g-1} g_i(u_1, u_2) = \frac{m_iI}{a_i + I}, \end{equation} $ (1.4)

    with

    $ \begin{equation*}\label{light-g-2} I(t) = I_{\rm in} \exp\left(-k_0 z_m-k_1 z_m u _1(t)-k_2 z_m u _2(t)\right). \end{equation*} $

    According to [20], we take the growth rate $\mu_{i}(Q_{ i})$ as follows

    $ \begin{equation*}\label{growth} \mu_{i}(Q_{ i}) = \mu_{\max, i}\frac{Q_{ i}-Q_{\min, i}}{Q_{\max, i}-Q_{\min, i}}, \end{equation*} $

    where $\mu_{\max, i}$ is the maximum specific growth rate of species $i$; $Q_{\min, i}$ is the minimum cellular carbon content required for growth of species $i$; $Q_{\max, i}$ is the maximum cellular carbon content of species $i$. From [3,4,5], for $i = 1, 2$, the growth rate $\mu_{i}(Q_{ i})$ can also take the forms :

    $ \begin{eqnarray*}\label{} &&\mu_{i}(Q_{i}) = \mu_{\infty, i}\left(1-\frac{Q_{\min, i}}{Q_{ i}}\right), \text{or}\nonumber \\ \\[-1.5ex] &&\mu_{i}(Q_{ i}) = \mu_{\infty, i} \frac{(Q_{ i}-Q_{\min, i})_+}{A_{ i}+(Q_{ i}-Q_{\min, i})_+}, \nonumber \end{eqnarray*} $

    where $(Q_{ i}-Q_{\min, i})_+$ is the positive part of $(Q_{ i}-Q_{\min, i})$ and $\mu_{\infty, i}$ is the maximal growth rate at infinite quotas (i.e., as $Q_{ i} \rightarrow \infty$) of the species $i$.

    According to [6,16], for $H = R, \ S$ and $i = 1, 2$, the uptake rate $f_{H i}(H, Q_{i})$ takes the form:

    $ \begin{eqnarray*}\label{eqn1.4-1} &&f_{H i}(H, Q_{ i}) = \rho_{H i}(Q_{ i})\frac{H}{K_{H i}+H}. \end{eqnarray*} $

    Here $H = R, \ S$ represents the concentration of the extracellular resource; $\rho_{H i}(Q_{ i})$ represents the maximal uptake rate of the species $i$; $K_{H i}$ is the half-saturation constant, the resource concentration at which uptake rate is half of the maximal rate. The maximal resource uptake rate, $\rho_{H i}(Q_{ i})$, is a decreasing, linear function of quota [8,16], which is defined by

    $ \begin{equation*}\label{uptake-1} \rho_{H i}(Q_{ i}) = \rho_{\max, H i}^{\rm high}-(\rho_{\max, H i}^{\rm high}-\rho_{\max, H i}^{\rm low})\frac{Q_{ i}-Q_{\min, i}}{Q_{\max, i}-Q_{\min, i}}, \end{equation*} $

    for $Q_{\min, i}\leq Q_{ i}\leq Q_{\max, i}$. In other words, the maximal rate of resource uptake, $\rho_{H i}(Q_{ i})$, varies between upper and lower bounds, $\rho_{\max, H i}^{\rm high}$ and $\rho_{\max, H i}^{\rm low}$, respectively, while quota varies between lower and upper bounds, $Q_{\min, i}$ and $Q_{\max, i}$, respectively (see, e.g., [8]). Cunningham and Nisbet [1,2] took $\rho_{H i}(Q_{ i})$ to be a constant. The authors in [20] put $\rho_{\max, H i}^{\rm high}\equiv\rho_{\max, H i}$, $\rho_{\max, H i}^{\rm low} = 0$, then $\rho_{H i}(Q_{ i})$ becomes

    $ \begin{equation*}\label{uptake-2} \rho_{H i}(Q_{ i}) = \rho_{\max, H i}\frac{Q_{\max, i}-Q_{ i}}{Q_{\max, i}-Q_{\min, i}}. \end{equation*} $

    System (1.1) also includes the fact that carbon is lost by respiration. The respiration rate is proportional to the size of the transient carbon pool [20]:

    $ \begin{equation*}\label{respiration} \gamma_{i}(Q_{ i}) = \gamma_{\max, i}\frac{Q_{ i}-Q_{\min, i}}{Q_{\max, i}-Q_{\min, i}}, \end{equation*} $

    where $Q_{\min, i}\leq Q_{ i}\leq Q_{\max, i}$, and $\gamma_{\max, i}$ is the maximum respiration rate of species $i$.

    In this whole paper, we always assume that the photosynthetic rate $g_i(u_1, u_2)$ takes the form in (1.2) or (1.4); the functions $\mu_{ i}(Q_{ i})$, $f_{H i}(H, Q_{ i})$, and $\gamma_{i}(Q_{ i})$ satisfy the following assumptions: ($i = 1, 2$ and $H = R, \ S$)

    $\rm(C1)$ $\mu_{ i}(Q_{i})$ is continuously differentiable for $Q_{ i}\ge Q_{\min, i}$, $ \mu_{ i}(Q_{\min, i}) = 0$, $\mu_{ i}(Q_{ i})\geq0$ and $\mu_{ i}'(Q_{ i})>0$ for $Q_{ i}\ge Q_{\min, i}$.

    $\rm(C2)$ $f_{H i}(H, Q_{ i})$ is continuously differentiable for $H>0$ and $Q_{ i}\ge Q_{\min, i}$, $f_{H i}(0, Q_{i}) = 0$, $f_{H i}(H, Q_{i})\ge 0$, $ \frac{\partial f_{H i}(H, Q_{i})}{\partial H}>0$ and $\frac{\partial f_{H i}(H, Q_{ i})}{\partial Q_{ i}}\leq 0$ for $H>0$ and $Q_{ i}\ge Q_{\min, i}$.

    $\rm(C3)$ $\gamma_{i}(Q_{i})$ is continuously differentiable for $Q_{ i}\ge Q_{\min, i}$, $ \gamma_{i}(Q_{\min, i}) = 0$, $\gamma_{i}(Q_{ i})\geq0$ and $\gamma_{i}'(Q_{ i})>0$ for $Q_{ i}\ge Q_{\min, i}$.

    The rest of the paper is organized as follows. Section 2 is devoted to the study of the single population model. In Section 3, we shall investigate the possibility of coexistence of the two competing species system (1.1). A brief discussion section completes this paper.

    In this section, we first investigate the extinction and persistence of the single population model. Mathematically, it simply means that we remove equations of $Q_2$ and $u_2$ from (1.1). In order to simplify notation, all subscripts are dropped in the remaining equations and the single population model takes the form:

    $ \begin{equation}\label{Droop-Co2-light-single-ori} \begin{cases} \frac{d R}{dt} = (R^{(0)}-R)D-f_{R}(R, Q)G(u)u+\gamma(Q)u-\omega_rR+\omega_sS, \\ \frac{dS}{dt} = (S^{(0)}-S)D-f_{S}(S, Q)G(u)u+\omega_rR-\omega_sS, \\ \frac{dQ}{dt} = f_{R}(R, Q)G(u)+f_{S}(S, Q)G(u)-\mu(Q)Q-\gamma(Q), \\ \frac{du}{dt} = \left[\mu(Q)-D\right]u, \\ R(0)\ge0, \ S(0)\ge0, \ u(0)\ge 0, \ Q(0)\ge Q_{\min}. \end{cases} \end{equation} $ (2.1)

    The specific growth rate $G(u)$ takes the form

    $ \begin{equation*}\label{light-g-single-2} G(u) = \frac{1}{z_m}\int_0^{z_m}\frac{m\tilde{I}(z)}{a+\tilde{I}(z)}dz = \frac{m}{\ln\left(\frac{I_{\rm in}}{\tilde{I}(z_m)}\right)}\ln\left(\frac{a+I_{\rm in}}{a+\tilde{I}(z_m)}\right), \end{equation*} $

    where

    $ \begin{equation*}\label{light-g-single-3} \tilde{I}(z) = I_{\rm in} \exp\left(-k_0 z-k z u (t)\right), \ 0\leq z \leq z_m, \end{equation*} $

    or

    $ \begin{equation*}\label{light-g-single-1} G(u) = \frac{mI}{a + I}, \ \mbox{with}\ I(t) = I_{\rm in} \exp\left(-k_0 z_m-k z_m u (t)\right). \end{equation*} $

    The feasible domain for system (2.1) takes the form

    $ \begin{equation}\label{Droop-Co2-light-single-ori-dom} \mathbf{X} = \{(R, S, Q, u)\in \mathbb{R}_{+}^4:Q \ge Q_{\min}\}. \end{equation} $ (2.2)

    Then it is easy to show that $\mathbf{X}$ is positively invariant for system (2.1). Next, we study the boundedness of solutions of (2.1):

    Lemma 2.1. Every solution $(R(t), S(t), Q(t), u(t))$ of system (2.1) exists globally. Furthermore, solutions of (2.1) are ultimately bounded on $\mathbf{X}$.

    Proof. By the continuation theorem, it suffices to prove that the solution of system (2.1) is bounded on finite time intervals. Let

    $ \begin{equation}\label{Droop-Co2-light-single-ori-glob-ex-1} \Theta(t) = R(t)+S(t)+u(t)Q(t). \end{equation} $ (2.3)

    Then $\Theta(t)$ satisfies

    $ \begin{equation}\label{Droop-Co2-light-single-ori-glob-ex-2} \frac{d \Theta(t)}{dt} = (R^{(0)}+S^{(0)}-\Theta(t))D. \end{equation} $ (2.4)

    Thus, $\Theta(t)$ is bounded on finite time intervals. This fact and the positivity of the solution imply that $R(t)$, $S(t)$, and $u(t)Q(t)$ are bounded on finite time intervals. Since $Q(t)\geq Q_{\min}$, we can also deduce that $u(t)$ is bounded on finite time intervals.

    It remains to show that $Q(t)$ is bounded on finite time intervals. Inspired by the ideas in [15,Proposition 3], we will investigate the dynamics of the variable $V(t) = \frac{1}{2}(Q(t))^2$. By the monotonicity of $f_{N}$ for $N = S, R$, the inequality $x\leq \frac 12(1+x^2)$, and the fact that $Q \ge Q_{\min}$, it is easy to see that $V(t)$ satisfies

    $ \begin{eqnarray}\label{Droop-Co2-light-single-ori-V-inequality} \frac{d V(t)}{dt}& = &\frac{d Q(t)}{dt}Q\nonumber\\ &\leq& G(0)[f_{R}(R(t), Q_{\min})+f_{S}(S(t), Q_{\min})]Q\nonumber\\ &\leq& \frac 12G(0)[f_{R}(R(t), Q_{\min})+f_{S}(S(t), Q_{\min})][1+Q^2] \\ & = & G(0)[f_{R}(R(t), Q_{\min})+f_{S}(S(t), Q_{\min})][\frac 12+V].\nonumber \end{eqnarray} $ (2.5)

    Since $R(t)$ and $S(t)$ in (2.5) are bounded on finite time intervals, we see that $V(t)$ is bounded on finite time intervals, and hence, so is $Q(t)$. Thus, every solution of system (2.1) exists globally.

    From (2.4), it is easy to see that

    $ \begin{equation}\label{Droop-Co2-light-single-ori-glob-ex-3} \lim\limits_{t\rightarrow \infty}\Theta(t) = R^{(0)}+S^{(0)}, \end{equation} $ (2.6)

    and hence, $\Theta(t)$ is ultimately bounded. This together with the positivity of the solution, we deduce that $R(t)$, $S(t)$, and $u(t)Q(t)$ are ultimately bounded. Using the fact $Q(t)\geq Q_{\min}$, we see that $u(t)$ is ultimately bounded. Finally, we show that $Q(t)$ is ultimately bounded. From (2.3), (2.6), and the third equation of (2.1), it follows that there exist $\tau_0>0$ and $\eta_0>0$ such that

    $ \frac{d Q}{dt}\leq G(0)[f_{R}(R^{(0)}+S^{(0)}, Q)+f_{S}(R^{(0)}+S^{(0)}, Q)]+\eta_0-\mu(Q)Q-\gamma(Q), \ t\geq \tau_0. $

    Then

    $ \lim\limits_{t\rightarrow \infty}Q(t)\leq Q^{\eta_0}, $

    where $Q^{\eta_0}$ is the positive root uniquely determined by

    $ G(0)[f_{R}(R^{(0)}+S^{(0)}, Q)+f_{S}(R^{(0)}+S^{(0)}, Q)]+\eta_0-\mu(Q)Q-\gamma(Q) = 0. $

    Thus, solutions of (2.1) are ultimately bounded on $\mathbf{X}$ since $\eta_0>0$ is independent of initial values.

    In order to find the species-free equilibrium of system (2.1), which corresponds to the absence of species, we put $u = 0$ in (2.1). Then we consider the following system

    $ \begin{equation}\label{Co2-only} \begin{cases} \frac{d R}{dt} = (R^{(0)}-R)D-\omega_rR+\omega_sS, \\ \frac{dS}{dt} = (S^{(0)}-S)D+\omega_rR-\omega_sS, \\ R(0)\ge 0, \ S(0)\ge 0. \end{cases} \end{equation} $ (2.7)

    It is easy to see that (2.7) is a cooperative/monotone system (see, e.g., [18]), and

    $ (R^{*}, S^{*}): = (\frac{DR^{(0)}+\omega_sR^{(0)}+\omega_sS^{(0)}}{D+\omega_r+\omega_s}, \frac{DS^{(0)}+\omega_rS^{(0)}+\omega_rR^{(0)}}{D+\omega_r+\omega_s}) $

    is the unique equilibrium for (2.7). For a monotone dynamical system, the unique steady state is globally asymptotically stable if and only if every forward orbit has compact closure (see [13,Theorem D]). By Lemma 2.1 and the above discussions, we have the following results:

    Lemma 2.2. The unique equilibrium $(R^{*}, S^{*})$ is globally asymptotically stable for (2.7) in $\mathbb{R}^2$.

    From Lemma 2.2, the species-free equilibrium of system (2.1), which we label $E_0$, is given by

    $ E_0 = (R, S, Q, u) = (R^{*}, S^{*}, Q^{*}, 0), $

    where $Q^{*}$ satisfies

    $ \begin{equation}\label{single-Qstar} f_{R}(R^{*}, Q^{*})G(0)+f_{S}(S^{*}, Q^{*})G(0)-\mu(Q^{*})Q^{*}-\gamma(Q^{*}) = 0 \end{equation} $ (2.8)

    The local stability of $E_0$ is determined by the Jacobian matrix of (2.1) at $E_0$, denoted by

    $ \begin{eqnarray*} J_0 = \left(\begin{array}{cccc} -D-\omega_r & \omega_s & 0 &-f_R(R^{*}, Q^{*})G(0)+\gamma(Q^{*}) \\ \omega_r& -D-\omega_s& 0 &-f_S(S^{*}, Q^{*})G(0) \\ \frac{\partial f_R(R^{*}, Q^{*})}{\partial R}G(0) & \frac{\partial f_S(S^{*}, Q^{*})}{\partial S}G(0) & j_{33} &[f_R(R^{*}, Q^{*})+f_S(S^{*}, Q^{*})]G'(0) \\ 0& 0 &0&\mu(Q^{*})-D \end{array}\right), \end{eqnarray*} $

    where

    $ \begin{equation*} j_{33} = [\frac{\partial f_R(R^{*}, Q^{*})}{\partial Q}+\frac{\partial f_S(S^{*}, Q^{*})}{\partial Q}]G(0)-[\mu(Q^{*})+\mu'(Q^{*})Q^{*}]-\gamma'(Q^{*}) \lt 0. \end{equation*} $

    It is easy to see that the eigenvalues of $J_0$ are $j_{33}$, $\mu(Q^{*})-D$, together with the eigenvalues of the following matrix

    $ \begin{eqnarray*} \tilde{J}_0 = \left(\begin{array}{cc} -D-\omega_r & \omega_s \\ \omega_r& -D-\omega_s \end{array}\right). \end{eqnarray*} $

    Since $j_{33}<0$ and the eigenvalues of $\tilde{J}_0$ are both negative, we see that the sign of $\mu(Q^{*})-D$ determines the stability of $E_0$. That is, $E_0$ is locally asymptotically stable if $\mu(Q^{*})-D<0$, and unstable if $\mu(Q^{*})-D>0$. Thus, we have proved the following result concerning with the local stability of $E_0$:

    Lemma 2.3. $E_0$ is locally asymptotically stable if $\mu(Q^{*})-D<0$, and unstable if $\mu(Q^{*})-D>0$.

    This subsection is devoted to the investigations of persistence of system (2.1).

    Theorem 2.4. Assume that $(R(t), S(t), Q(t), u(t))$ is the unique nonnegative solution of system (2.1), for all $t\in[0, \infty)$, with the initial value $(R^0, S^0, Q^0, u^0)\in \mathbf{X}$. If $\mu(Q^{*})-D>0$, then system (2.1) is uniformly persistent in the sense that there exists a $\xi>0$ such that

    $ \liminf\limits_{t\rightarrow\infty}u(t)\geq \xi, \ provided~ that\ u^0\neq 0. $

    Further, system (2.1) admits at least one positive equilibrium $(\hat{R}, \hat{S}, \hat{Q}, \hat{u})$.

    Proof. Recall that $\mathbf{X}$ is defined in (2.2). Let

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-dom-interior} \mathbf{X}_0 = \{(R, S, Q, u)\in \mathbf{X}:u \gt 0\}, \end{equation*} $

    and

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-dom-bd} \partial \mathbf{X}_0: = \mathbf{X}\backslash\mathbf{X}_0: = \{(R, S, Q, u)\in \mathbf{X}: u = 0 \}. \end{equation*} $

    It is easy to see that both $\mathbf{X}$ and $\mathbf{X}_0$ are positively invariant for system (2.1), and $\partial \mathbf{X}_0$ is relatively closed in $\mathbf{X}$. Furthermore, system (2.1) is point dissipative (see Lemma 2.1). Let $\Phi(t):\mathbf{X}\rightarrow \mathbf{X}$ be the solution maps associated with system (2.1). Set $\tilde{M}_\partial: = \{P\in \partial\mathbf{X}_0: \Phi(t)P\in \partial\mathbf{X}_0, \ \forall \ t\geq0\}$ and $\tilde{\omega}(P)$ be the omega limit set of the orbit $\tilde{O}^{+}(P): = \{\Phi(t)P:t\geq 0\}$. We show the following claim.

    Claim: $\tilde{\omega}(P) = \{E_0\}, \ \forall \ P\in \tilde{M}_{\partial}$.

    Since $P\in \tilde{M}_{\partial}$, we have $\Phi(t)P\in \partial\mathbf{X}_0, \ \forall \ t\geq 0$. Thus, $u(t, P) = 0, \ \forall \ t\geq 0.$ In view of the first two equations of (2.1), it follows that $(R(t, P), S(t, P))$ satisfies (2.7), $\forall \ t\geq 0$. Then Lemma 2.2 implies that

    $ \lim\limits_{t\rightarrow \infty}(R(t, P), S(t, P)) = (R^{*}, S^{*}). $

    Then, the equation for $Q(t)$ in (2.1) is asymptotic to

    $ \frac{dQ}{dt} = f_{R}(R, Q)G(0)+f_{S}(S, Q)G(0)-\mu(Q)Q-\gamma(Q). $

    From the theory for asymptotically autonomous semiflows (see, e.g., [19,Corollary 4.3]), it follows that $\lim_{t\rightarrow \infty}Q(t, P) = Q^{*}$, where $Q^{*}$ is given in (2.8). Hence, the claim is proved.

    Let $\eta_1: = \frac{1}{2}\left(\mu(Q^{*})-D\right)>0$. Then it follows from the continuity of $\mu(Q)$ that there exists $\sigma_1>0$ such that

    $ \begin{equation}\label{Droop-Co2-light-single-ori-weak-repeller-1} \mu(Q) \gt \mu(Q^{*})-\eta_1, \ \forall \ |Q-Q^{*}| \lt \sigma_1. \end{equation} $ (2.9)

    Claim: $E_0$ is a uniform weak repeller for $\mathbf{X}_0$ in the sense that

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-weak-repeller-2} \limsup\limits_{t\rightarrow\infty}|\Phi_t(P)-E_0|\geq \sigma_1, \ \forall\ P\in \mathbf{X}_0. \end{equation*} $

    Suppose not. Then there exists a $P\in \mathbf{X}_0$ such that

    $ \limsup\limits_{t\rightarrow\infty}|\Phi_t(P)-E_0| \lt \sigma_1. $

    Thus, there exists a $\tau_1>0$ such that

    $ |Q(t, P)-Q^{*}| \lt \sigma_1, \ \forall \ t\geq \tau_1. $

    This and (2.9) imply that

    $ \mu(Q(t, P))-D \gt \mu(Q^{*})-D-\eta_1 = \eta_1, \ \forall \ t\geq \tau_1. $

    From this inequality and the fourth equation of (2.1), we have

    $ \begin{equation*}\label{} \frac{du(t, P)}{dt} \gt \eta_1u(t, P), \ \forall \ t\geq \tau_1, \end{equation*} $

    which shows that $\lim_{t\rightarrow\infty}u(t, P) = \infty$, a contradiction.

    Therefore, $E_0$ is isolated in $\mathbf{X}$ and $\tilde{W}^s(E_0)\cap \mathbf{X}_{0} = \emptyset$, where $\tilde{W}^{s}(E_0)$ is the stable set of $E_0$ (see [21]). Since $\Phi_t:\mathbf{X} \rightarrow \mathbf{X}$ is point dissipative and compact, we conclude from [21,Theorem 1.1.3] that there exists a global attractor $\mathcal{A}$ for $\Phi_t$ in $\mathbf{X}$. By [21,Theorem 1.3.1] on strong repellers, $\Phi_t:\mathbf{X} \rightarrow \mathbf{X}$ is uniformly persistent with respect to $(\mathbf{X}_0, \partial\mathbf{X}_0)$. It follows from [21,Theorem 1.3.6] that there exists a global attractor $\mathcal{A}_0$ for $\Phi_t$ in $\mathbf{X}_0$ and $\Phi_t$ admits at least one fixed point

    $ (\hat{R}, \hat{S}, \hat{Q}, \hat{u})\in \mathbf{X}_0. $

    Then $\hat{Q}\geq Q_{\min}>0$, $\hat{u}>0$, and $(\hat{R}, \hat{S})$ satisfies

    $ \begin{equation}\label{Droop-Co2-light-single-ori-weak-repeller-3} \begin{cases} (R^{(0)}-\hat{R})D-f_{R}(\hat{R}, \hat{Q})g(\hat{u})\hat{u}+\gamma(\hat{Q})\hat{u}-\omega_r\hat{R}+\omega_s\hat{S} = 0, \\ (S^{(0)}-\hat{S})D-f_{S}(\hat{S}, \hat{Q})g(\hat{u})\hat{u}+\omega_r\hat{R}-\omega_s\hat{S} = 0. \end{cases} \end{equation} $ (2.10)

    In view of (2.10), we deduce that $\hat{R}>0, \ \hat{S}>0$. Thus, $(\hat{R}, \hat{S}, \hat{Q}, \hat{u})$ is a positive steady-state solution for (2.1). We complete the proof.

    In this subsection, we neglect the effect of respiration and investigate the extinction of system (2.1). Putting the respiration rate to be zero,

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extin-ass} \gamma(Q)\equiv0, \ \forall\ Q\geq Q_{\min}. \end{equation} $ (2.11)

    Then we have the following result:

    Theorem 2.5. Suppose $(R(t), S(t), Q(t), u(t))$ is the unique nonnegative solution of system (2.1), for all $t\in[0, \infty)$, with the initial value $(R^0, S^0, Q^0, u^0)\in \mathbf{X}$. Assume that (2.11) holds. If $\mu(Q^{*})-D<0$, then system (2.1) is washed out, that is,

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-extinction-eq} \lim\limits_{t\rightarrow\infty}(R(t), S(t), Q(t), u(t)) = (R^{*}, S^{*}, Q^{*}, 0). \end{equation*} $

    Proof. For $\eta>0$, we assume that $Q^{\eta}$ is the unique root of

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extinction-Q} G(0)[f_{R}(R^{*}, Q)+f_{S}(S^{*}, Q)]+\eta-\mu(Q)Q = 0. \end{equation} $ (2.12)

    Recall that $Q^{*}$ is the unique root of (2.8). Since

    $ \lim\limits_{\eta\rightarrow 0}[\mu(Q^{\eta})-D] = \mu(Q^{*})-D\ \mbox{and}\ \frac{-1}{3}[\mu(Q^{*})-D] \gt 0, $

    we may find an $\eta_2>0$ such that

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extinction-QQ0} \mu(Q^{\eta_2})-D \lt [\mu(Q^{*})-D]+\frac{-1}{3}[\mu(Q^{*})-D] = \frac{2}{3}[\mu(Q^{*})-D]. \end{equation} $ (2.13)

    On the other hand, by the continuity, we may find a $\sigma_2>0$ such that

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extinction-f-ineq} \begin{cases} f_{R}(R^{*}+\sigma_2, Q) \lt f_{R}(R^{*}, Q)+\frac{\eta_2}{2G(0)}, \\ f_{S}(S^{*}+\sigma_2, Q) \lt f_{S}(S^{*}, Q)+\frac{\eta_2}{2G(0)}, \\ \mu(Q^{\eta_2}+\sigma_2) \lt \mu(Q^{\eta_2})+\frac{-1}{3}[\mu(Q^{*})-D] \end{cases} \end{equation} $ (2.14)

    In view of the assumption (2.11) and the first two equations of system (2.1), it follows that

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-extinction-RS} \begin{cases} \frac{d R}{dt}\leq (R^{(0)}-R)D-\omega_rR+\omega_sS, \\ \frac{dS}{dt}\leq (S^{(0)}-S)D+\omega_rR-\omega_sS. \end{cases} \end{equation*} $

    By the comparison arguments and Lemma 2.2, we have

    $ \lim\limits_{t\rightarrow \infty}(R(t), S(t))\leq (R^{*}, S^{*}). $

    Then there exists a $\tau_2>0$ such that

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-extinction-RS-long} R(t)\leq R^{*}+\sigma_2, \ S(t)\leq S^{*}+\sigma_2, \ \forall \ t\geq \tau_2. \end{equation*} $

    Then it follows from the third equation of (2.1) that

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extinction-Q-long} \frac{dQ}{dt}\leq G(0)[f_{R}(R^{*}+\sigma_2, Q)+f_{S}(S^{*}+\sigma_2, Q)]-\mu(Q)Q, \ \forall \ t\geq \tau_2. \end{equation} $ (2.15)

    In view of the first two inequalities of (2.14) and (2.15), we have

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-extinction-Q-long-1} \frac{dQ}{dt}\leq G(0)[f_{R}(R^{*}, Q)+f_{S}(S^{*}, Q)]+\eta_2-\mu(Q)Q, \ \forall \ t\geq \tau_2. \end{equation*} $

    Using the comparison arguments, we have

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extinction-Q-long-2} \lim\limits_{t\rightarrow \infty}Q(t)\leq Q^{\eta_2}, \end{equation} $ (2.16)

    where $Q^{\eta_2}$ is the unique root of (2.12) with $\eta = \eta_2$. Thus, we may find a $\tau_3>0$ such that

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-extinction-Q-long-3} Q(t)\leq Q^{\eta_2}+\sigma_2, \ \forall \ t\geq \tau_3, \end{equation*} $

    and hence,

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extinction-Q-long-4} \mu(Q(t))\leq \mu(Q^{\eta_2}+\sigma_2), \ \forall \ t\geq \tau_3. \end{equation} $ (2.17)

    In view of the third inequality of (2.14) and (2.17), we have

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extinction-Q-long-5} \mu(Q(t))\leq \mu(Q^{\eta_2})+\frac{-1}{3}[\mu(Q^{*})-D], \ \forall \ t\geq \tau_3. \end{equation} $ (2.18)

    By (2.13), (2.18) together with the fourth equation of (2.1), it follows that

    $ \begin{equation}\label{Droop-Co2-light-single-ori-extinction-u-long} \frac{du}{dt} = \left[\mu(Q(t))-D\right]u\leq \frac{1}{3}[\mu(Q^{*})-D] u, \ \forall \ t\geq \tau_3. \end{equation} $ (2.19)

    Since $\mu(Q^{*})-D<0$ and (2.19), we have

    $ \begin{equation*}\label{Droop-Co2-light-single-ori-extinction-u-long-1} \lim\limits_{t\rightarrow \infty}u(t) = 0. \end{equation*} $

    Then $(R(t), S(t))$ is asymptotic to (2.7). By Lemma 2.2 and the theory for asymptotically autonomous semiflows (see, e.g., [19,Corollary 4.3]), it follows that

    $ \lim\limits_{t\rightarrow \infty}(R(t), S(t)) = (R^{*}, S^{*}). $

    Similarly, $Q(t)$ is asymptotic to

    $ \frac{dQ}{dt} = f_{R}(R^{*}, Q)G(0)+f_{S}(S^{*}, Q)G(0)-\mu(Q)Q, $

    and $\lim_{t\rightarrow \infty}Q(t) = Q^{*}$, where $Q^{*}$ is given in (2.8) with $\gamma(Q)\equiv 0$. We complete the proof.

    In this section, we shall concentrate on the study of coexistence of system (1.1). The subsequent discussions will reveal that two semi-trivial steady-state solutions of system (1.1), corresponds to the presence of one of the species and the absence of the other species, are not necessarily unique. This makes our analysis more difficult. Fortunately, we can adopt the ideas developed in [12,Section 4] to overcome this difficulty.

    The trivial steady-state solution of (1.1), labeled $\mathcal{E}_0$, corresponds to the absence of both species. It is given by

    $ \mathcal{E}_0 = (R, S, Q_1, u_1, Q_2, u_2) = (R^{*}, S^{*}, Q_1^{*}, 0, Q_2^{*}, 0), $

    where $(R^{*}, S^{*})$ is given in Lemma 2.2, and $Q_i^{*}$ satisfies

    $ \begin{equation*}\label{Qstar-i} f_{Ri}(R^{*}, Q_i^{*})G_i(0)+f_{Si}(S^{*}, Q_i^{*})G_i(0)-\mu_i(Q_i^{*})Q_i^{*}-\gamma_i(Q_i^{*}) = 0, \ i = 1, 2, \end{equation*} $

    where $G_1(u_1) = g_1(u_1, 0)$ and $G_2(u_2) = g_2(0, u_2)$.

    In order to determine the semi-trivial steady-state solutions of (1.1), we need the following single population system associated with the growth of species $i$:

    $ \begin{equation}\label{Droop-Co2-light-ori-species-i} \begin{cases} \frac{d R}{dt} = (R^{(0)}-R)D-f_{Ri}(R, Q_{i})G_i(u_i)u_{i}+\gamma_i(Q_i)u_i-\omega_rR+\omega_sS, \\ \frac{dS}{dt} = (S^{(0)}-S)D-f_{Si}(S, Q_{i})G_i(u_i)u_{i}+\omega_rR-\omega_sS, \\ \frac{dQ_{i}}{dt} = f_{Ri}(R, Q_{i})G_i(u_i)+f_{Si}(S, Q_{i})G_i(u_i)-\mu_{i}(Q_{i})Q_{i}-\gamma_i(Q_i), \\ \frac{du_{i}}{dt} = \left[\mu_{i}(Q_{i})-D\right]u_{i}, \\ R(0)\ge0, \ S(0)\ge0, \ u_{i}(0)\ge 0, \ Q_{i}(0)\ge Q_{\min, i}, \ i = 1, 2. \end{cases} \end{equation} $ (3.1)

    Next, we shall summarize the result of system (3.1). By Theorem 2.4, for $i = 1, 2$, system (3.1) admits at least one positive equilibrium and we may assume that $A_i^{0}\subset$$\rm{Int}$$\mathbb{R}_{+}^4$ is a global attractor of the semiflows generated by system (3.1), under the condition that $\mu_i(Q_i^{*})-D>0$. One of the semi-trivial steady-state solutions of (1.1), labeled $\mathcal{E}_1$, corresponds to the presence of species 1 and the absence of species 2. It is given by

    $ \mathcal{E}_1 = (R, S, Q_1, u_1, Q_2, u_2) = (\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1, \hat{Q}_2, 0), $

    where $(\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1)\in A_1^{0}$ is a positive equilibrium of system (3.1) with $i = 1$, which is not necessarily unique. Here, $\hat{Q}_2 = \hat{Q}_2(\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1)$ is the root of

    $ \begin{equation}\label{Droop-Co2-light-ori-species-E1-Q} f_{R2}(\hat{R}_1, Q_{2})g_2(\hat{u}_1, 0)+f_{S2}(\hat{S}_1, Q_{2})g_2(\hat{u}_1, 0)-\mu_{2}(Q_{2})Q_{2}-\gamma_2(Q_2) = 0. \end{equation} $ (3.2)

    Inspired by the arguments in [12,Section 4], we assume that

    $ \begin{equation}\label{Droop-Co2-light-ori-species-E1-Q-min} \hat{Q}_2^{\min} = \inf\{\hat{Q}_2(\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1):(\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1)\in A_1^{0} \}. \end{equation} $ (3.3)

    The other semi-trivial steady-state solution of (1.1), labeled $\mathcal{E}_2$, corresponds to the presence of species 2 and the absence of species 1. It is given by

    $ \mathcal{E}_2 = (R, S, Q_1, u_1, Q_2, u_2) = (\check{R}_2, \check{S}_2, \check{Q}_1, 0, \check{Q}_2, \check{u}_2), $

    where $(\check{R}_2, \check{S}_2, \check{Q}_2, \check{u}_2)\in A_2^{0}$ is a positive equilibrium of system (3.1) with $i = 2$, which is not necessarily unique. Here, $\check{Q}_1 = \check{Q}_1(\check{R}_2, \check{S}_2, \check{Q}_2, \check{u}_2)$ is the root of

    $ \begin{equation*}\label{Droop-Co2-light-ori-species-E2-Q} f_{R1}(\check{R}_2, Q_{1})g_1(0, \check{u}_2)+f_{S1}(\check{S}_2, Q_{1})g_1(0, \check{u}_2)-\mu_{1}(Q_{1})Q_{1}-\gamma_1(Q_1) = 0. \end{equation*} $

    Similarly, we assume

    $ \begin{equation}\label{Droop-Co2-light-ori-species-E2-Q-min} \check{Q}_1^{\min} = \inf\{\check{Q}_1(\check{R}_2, \check{S}_2, \check{Q}_2, \check{u}_2):(\check{R}_2, \check{S}_2, \check{Q}_2, \check{u}_2)\in A_2^{0} \}. \end{equation} $ (3.4)

    The feasible domain for system (1.1) takes the form

    $ \begin{equation*}\label{Droop-Co2-light-ori-dom} \mathbf{Y} = \{(R, S, Q_1, u_1, Q_2, u_2)\in \mathbb{R}_{+}^6:Q \ge Q_{\min, i}, \ i = 1, 2\}. \end{equation*} $

    Then it is easy to show that $\mathbf{Y}$ is positively invariant for system (1.1). By the similar arguments in Lemma 2.1, we can show the following result:

    Lemma 3.1. Every solution $(R(t), S(t), Q_1(t), u_1(t), Q_2(t), u_2(t))$ of system (1.1) exists globally. Furthermore, solutions of (1.1) are ultimately bounded on $\mathbf{Y}$.

    Assume that $\Psi(t):\mathbf{Y}\rightarrow\mathbf{Y}$ is the semiflow associated with system (1.1). Let

    $ \begin{equation*}\label{Droop-Co2-light-ori-dom-interior} \mathbf{Y}_0 = \{(R, S, Q_1, u_1, Q_2, u_2)\in \mathbf{Y}:u_1 \gt 0\ \mbox{and}\ u_2 \gt 0\}, \end{equation*} $

    and

    $ \begin{equation*}\label{Droop-Co2-light-ori-dom-bd} \partial \mathbf{Y}_0: = \mathbf{Y}\backslash\mathbf{Y}_0: = \{(R, S, Q_1, u_1, Q_2, u_2)\in \mathbf{Y}: u_1 = 0\ \mbox{or}\ u_2 = 0\}. \end{equation*} $

    Following the ideas in [12,Section 4], we assume that $\mathcal{M}_0 = \{\mathcal{E}_0\}$,

    $ \begin{equation*}\label{} \mathcal{M}_1 = \{(\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1, \hat{Q}_2, 0)\in \mathbf{Y}: (\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1)\in A_1^{0}\ \mbox{and}\ \hat{Q}_2\ \ \mbox{is defined by} \ (3.2)\}, \end{equation*} $

    and

    $ \begin{equation*}\label{} \mathcal{M}_2 = \{(\check{R}_2, \check{S}_2, \check{Q}_1, 0, \check{Q}_2, \check{u}_2)\in \mathbf{Y}: (\check{R}_2, \check{S}_2, \check{Q}_2, \check{u}_2)\in A_2^{0}\ \mbox{and}\ \check{Q}_1\ \ \mbox{is defined by} \ (3.4)\}. \end{equation*} $

    One can easily to use "the method of proof by contradiction" to deduce the following result:

    Lemma 3.2. Let $\mu_i(Q_i^{*})-D>0$, for some $i\in \{1, 2\}$. Then $\mathcal{M}_0$ is a uniform weak repeller in the sense that there exists a $\delta_0>0$ such that

    $ \begin{equation*}\label{Droop-Co2-light-ori-M0-repeller-eq} \limsup\limits_{t\rightarrow\infty}|\Psi(t) \mathbf{v}^0-\mathcal{M}_0|\geq \delta_0, \ for~all\ \mathbf{v}^0\in \mathbf{Y}_0. \end{equation*} $

    Next, we shall use the strategy in [12,Lemma 4.2] to show the following result:

    Lemma 3.3. Let $\mu_i(Q_i^{*})-D>0$, for each $i\in \{1, 2\}$. If $\mu_2(\hat{Q}_2^{\min})-D>0$, then $\mathcal{M}_1$ is a uniform weak repeller in the sense that there exists a $\delta_1>0$ such that

    $ \begin{equation}\label{Droop-Co2-light-ori-M1-repeller-eq} \limsup\limits_{t\rightarrow\infty}\text{dist}(\Psi(t) \mathbf{v}^0, \mathcal{M}_1)\geq \delta_1, \ for ~all\ \mathbf{v}^0\in \mathbf{Y}_0. \end{equation} $ (3.5)

    Proof. Let

    $ \begin{equation*}\label{} B_1 = \{\hat{Q}_2 = \hat{Q}_2(\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1): (\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1)\in A_1^{0}\ \mbox{and}\ \hat{Q}_2\ \ \mbox{is defined by} \ (3.2)\}. \end{equation*} $

    Then

    $ \hat{Q}_2^{\min} = \inf\{\hat{Q}_2(\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1):\hat{Q}_2(\hat{R}_1, \hat{S}_1, \hat{Q}_1, \hat{u}_1)\in B_1\}. $

    Setting

    $ \epsilon_1 = \frac{1}{2}[\mu_2(\hat{Q}_2^{\min})-D] \gt 0. $

    Define $\mathcal{G}:B_1\rightarrow \mathbb{R}$ by

    $ \mathcal{G}(\phi) = \mu_2(\phi), \ \phi\in B_1. $

    We may find a $\delta_1>0$ such that

    $ \text{dist}(\mathcal{G}(\phi), \mathcal{G}(B_1)) \lt \epsilon_1, $

    whenever $\phi\in \mathbb{R}$ with $\text{dist}(\phi, B_1)<\delta_1$. Since $B_1$ is compact, it follows that for any $\phi\in \mathbb{R}$ with $\text{dist}(\phi, B_1)<\delta_1$, there exists $\phi_*\in B_1$ with $\phi_*$ depending on $\phi$ such that

    $ |\mathcal{G}(\phi)-\mathcal{G}(\phi_*)| = \text{dist}(\mathcal{G}(\phi), \mathcal{G}(B_1)) \lt \epsilon_1. $

    Thus, we have

    $ \begin{equation}\label{Droop-Co2-light-ori-M1-repeller-eq1} |\mu_2(\phi)-\mu_2(\phi_*)| = |\mathcal{G}(\phi)-\mathcal{G}(\phi_*)| \lt \epsilon_1, \end{equation} $

    whenever $\phi\in \mathbb{R}$ with $\text{dist}(\phi, B_1)<\delta_1$.

    Suppose that (3.5) is not true. Then there exists $\mathbf{v}^0\in\mathbf{Y}_0$ such that

    $ \begin{equation*}\label{} \limsup\limits_{t\rightarrow\infty}\text{dist}(\Psi(t) \mathbf{v}^0, \mathcal{M}_1) \lt \delta_1. \end{equation*} $

    This implies that

    $ \begin{equation}\label{Droop-Co2-light-ori-M1-repeller-eq2} \limsup\limits_{t\rightarrow\infty}\text{dist}(Q_2(t), B_1) \lt \delta_1\ \ \mbox{and}\ \ \limsup\limits_{t\rightarrow\infty}|u_2(t)| \lt \delta_1. \end{equation} $ (3.7)

    From the first inequality of (3.7), we can choose $t_1>0$ such that

    $ \begin{equation*}\label{Droop-Co2-light-ori-M1-repeller-eq3} \text{dist}(Q_2(t), B_1) \lt \delta_1, \ \forall \ t\geq t_1. \end{equation*} $

    By (3.6), it follows that there exists $\phi_*^t\in B_1$ such that

    $ \begin{equation*}\label{Droop-Co2-light-ori-M1-repeller-eq4} |\mu_2(Q_2(t))-\mu_2(\phi_*^t)| \lt \epsilon_1, \ \forall \ t\geq t_1, \end{equation*} $

    which implies that

    $ \begin{equation*}\label{Droop-Co2-light-ori-M1-repeller-eq5} \mu_2(Q_2(t))-D \gt \mu_2(\phi_*^t)-D-\epsilon_1\geq \mu_2(\hat{Q}_2^{\min})-D-\epsilon_1 = \epsilon_1, \ \forall \ t\geq t_1. \end{equation*} $

    From the sixth equation of (1.1), we have

    $ \begin{equation*}\label{Droop-Co2-light-ori-M1-repeller-eq6} \frac{du_{2}(t)}{dt} = \left[\mu_{2}(Q_{2}(t))-D\right]u_{2}(t) \gt \epsilon_1u_{2}(t), \ \forall \ t\geq t_1. \end{equation*} $

    We deduce that $\lim_{t\rightarrow \infty}u_{2}(t) = \infty$ since $\epsilon_1>0$ and $u_{2}(t_1)>0$. This contradicts the second inequality of (3.7) and we complete the proof.

    By the same arguments in Lemma 3.3, the following result holds:

    Lemma 3.4. Let $\mu_i(Q_i^{*})-D>0$, for each $i\in \{1, 2\}$. If $\mu_1(\check{Q}_1^{\min})-D>0$, then $\mathcal{M}_2$ is a uniform weak repeller in the sense that there exists a $\delta_2>0$ such that

    $ \begin{equation*}\label{Droop-Co2-light-ori-M2-repeller-eq} \limsup\limits_{t\rightarrow\infty}\text{dist}(\Psi(t) \mathbf{v}^0, \mathcal{M}_2)\geq \delta_2, \ for~ all\ \mathbf{v}^0\in \mathbf{Y}_0. \end{equation*} $

    Now we are in a position to prove the main result of this paper.

    Theorem 3.5. Assume that $(R(t), S(t), Q_1(t), u_1(t), Q_2(t), u_2(t))$ is the unique solution of (1.1) with the initial value $(R(0), S(0), Q_1(0), u_1(0), Q_2(0), u_2(0))\in \mathbf{Y}$. Let $\mu_i(Q_i^{*})-D>0, \ \forall \ i = 1, 2$, $\mu_2(\hat{Q}_2^{\min})-D>0$, and $\mu_1(\check{Q}_1^{\min})-D>0$. Then system (1.1) is uniformly persistent with respect to $(\mathbf{Y}_0, \partial \mathbf{Y}_0)$ in the sense that there is a positive constant $\zeta>0$ such that if $u_1(0)\neq 0$ and $u_2(0)\neq 0$, we have

    $ \begin{equation*}\label{roop-Co2-light-ori-co-eq} \liminf\limits_{t\rightarrow \infty}u_i(t)\geq \zeta, \ i = 1, 2. \end{equation*} $

    Furthermore, system (1.1) admits at least one (componentwise) positive equilibrium.

    Proof. Recall that $\Psi(t):\mathbf{Y}\rightarrow\mathbf{Y}$ is the semiflow associated with system (1.1). It is easy to see that $\Psi(t)\mathbf{Y}_0\subset\mathbf{Y}_0$. Since solutions of the system (1.1) are ultimately bounded (see Lemma 3.1), it follows that $\Psi(t)$ is point dissipative and compact, and hence, $\Psi(t)$ admits a global attractor (see, e.g., [21,Theorem 1.1.3]). Let

    $ M_{\partial}: = \{\mathbf{v}^0 \in \partial \mathbf{Y}_{0}:\Psi(t)\mathbf{v}^0\in \partial \mathbf{Y}_{0}, \forall \ t\geq 0\}, $

    and $\omega(\mathbf{v}^0)$ be the omega limit set of the orbit $O^{+}(\mathbf{v}^0): = \{\Psi(t)\mathbf{v}^0:t\geq 0\}$.

    Claim: $\bigcup_{\mathbf{v}^0\in M_{\partial}}\omega(\mathbf{v}^0)\subset \mathcal{M}_0\cup \mathcal{M}_1 \cup\mathcal{M}_2$.

    For any given $\mathbf{v}^0: = (R^0, S^0, Q_1^0, u_1^0, Q_2^0, u_2^0)\in M_{\partial}$, we have $\mathbf{v}^0\in \partial \mathbf{Y}_{0}$ and $\Psi(t)\mathbf{v}^0\in \partial \mathbf{Y}_{0}, \ \forall \ t\geq 0$. We discuss the following three subcases:

    (ⅰ) If $u_1^0 = 0, \ u_2^0 = 0$, then we have $u_1(t, \mathbf{v}^0) = 0$ and $u_2(t, \mathbf{v}^0) = 0$, $\forall\ t\geq0$. Thus, it is easy to see that $\lim_{t\rightarrow \infty}\Psi(t)\mathbf{v}^0 = \mathcal{E}_0$.

    (ⅱ) If $u_1^0 \neq 0, \ u_2^0 = 0$, then we have $u_1(t, \mathbf{v}^0)> 0$ and $u_2(t, \mathbf{v}^0) = 0$, $\forall\ t\geq0$. Then $(R(t, \mathbf{v}^0), S(t, \mathbf{v}^0), Q_1(t, \mathbf{v}^0), u_1(t, \mathbf{v}^0))$ satisfies system (3.1) with $i = 1$. Since $\mu_1(Q_1^{*})-D>0$, it follows from Theorem 2.4 that

    $ (R(t, \mathbf{v}^0), S(t, \mathbf{v}^0), Q_1(t, \mathbf{v}^0), u_1(t, \mathbf{v}^0)) $

    will eventually enter the global attractor $A_1^{0}\subset$$\rm{Int}$$\mathbb{R}_{+}^4$, and hence, $\Psi(t)\mathbf{v}^0$ will eventually enter $\mathcal{M}_1$.

    (ⅲ) If $u_1^0 = 0, \ u_2^0\neq 0$, then we have $u_1(t, \mathbf{v}^0) = 0$ and $u_2(t, \mathbf{v}^0)> 0$, $\forall\ t\geq0$. Then $(R(t, \mathbf{v}^0), S(t, \mathbf{v}^0), Q_2(t, \mathbf{v}^0), u_2(t, \mathbf{v}^0))$ satisfies system (3.1) with $i = 2$. Since $\mu_2(Q_2^{*})-D>0$, it follows from Theorem 2.4 that

    $ (R(t, \mathbf{v}^0), S(t, \mathbf{v}^0), Q_2(t, \mathbf{v}^0), u_2(t, \mathbf{v}^0)) $

    will eventually enter the global attractor $A_2^{0}\subset$$\rm{Int}$$\mathbb{R}_{+}^4$, and hence, $\Psi(t)\mathbf{v}^0$ will eventually enter $\mathcal{M}_2$.

    The proof of the claim is complete.

    By Lemma 3.2, Lemma 3.3 and Lemma 3.4, it follows that for $i = 0, 1, 2$, $\mathcal{M}_i$ is a uniform weak repeller for $\mathbf{Y}_{0}$ in the sense that there exists $\delta_i>0$ such that

    $ \begin{equation*}\label{} \limsup\limits_{t\rightarrow\infty}\text{dist}(\Psi(t) \mathbf{v}^0, \mathcal{M}_i)\geq \delta_i, \ \mbox{for all}\ \mathbf{v}^0\in \mathbf{Y}_0. \end{equation*} $

    Note that $\mathcal{M}_0$, $\mathcal{M}_1$, and $\mathcal{M}_2$ are pairwise disjoint, compact and isolated invariant sets for $\Psi(t)$ in $\partial\mathbf{Y}_0$. Further, each $\mathcal{M}_i$ is isolated in $\mathbf{Y}$ and $\mathcal{W}^s(\mathcal{M}_i)\cap \mathbf{Y}_{0} = \emptyset$, where $\mathcal{W}^{s}(\mathcal{M}_i)$ is the stable set of $\mathcal{M}_i$ (see [21]). It is easy to see that no subsets of $\mathcal{M}_0$, $\mathcal{M}_1$, and $\mathcal{M}_2$ forms a cycle in $\partial\mathbf{Y}_0$. By [21,Theorem 1.3.1] on strong repellers, $\Psi(t):\mathbf{Y}\rightarrow\mathbf{Y}$ is uniformly persistent with respect to $(\mathbf{Y}_0, \partial\mathbf{Y}_0)$. It follows from [21,Theorem 1.3.6] that there exists a global attractor $\hat{\mathcal{A}}_0$ for $\Psi(t)$ in $\mathbf{Y}_0$ and $\Psi(t)$ has at least one fixed point

    $ (\tilde{R}, \tilde{S}, \tilde{Q}_1, \tilde{u}_1, \tilde{Q}_2, \tilde{u}_2)\in \mathbf{Y}_0. $

    Thus, $(\tilde{R}, \tilde{S}, \tilde{Q}_1, \tilde{u}_1, \tilde{Q}_2, \tilde{u}_2)$ is a positive steady-state solution for system (1.1). This completes the proof.

    In this paper, we study the chemostat-type system (1.1) modeling the interactions of two species competing for "CO2" (dissolved CO2 and carbonic acid), "CARB" (bicarbonate and carbonate ions), and light in a spatially homogeneous water column. Our mathematical model presented in this paper is inspired by the recent works [17,20]. In fact, system (1.1) is a modified version of the model in the Supplementary Information of [20], where the specific growth rate of the competing species $i$ depends on their stored cellular carbon content (quota) $Q_i$. The dynamics of $Q_i$ can be affected by the uptake rates of inorganic carbon, "CO2" ($R(t)$) and "CARB" ($S(t)$), photosynthetic activity ($g_i(u_1, u_2)$), and respiration ($\gamma_i(Q_i)$). The resources "CO2" and "CARB" are substitutable in their effects on algal growth, which also involves a very complex processes. In order to make our system (1.1) analytically tractable, we have adopted the ideas in [17] to assume that carbonic acid loses a proton to become bicarbonate at the rate $\omega_r$, and the rate of the reverse reaction is denoted by $\omega_s$, which simplifies the complex processes of "CO2" and "CARB" involved.

    Solutions of both the two-species system (1.1) and its single-species sub-system (2.1) follow mass conservation laws, and are eventually bounded (see Lemma 2.1 and Lemma 3.1). Persistence of a single species depends on the sign of $\mu(Q^{*})-D$ (see Theorem 2.4), where $Q^{*}$ is given in (2.8). Biologically, $Q^*$ represents the quota that a species can obtain when the inorganic carbon concentration is at its long-term upper bound $(R^{*}, S^{*})$, which is the unique equilibrium for system (2.7) governing the available inorganic carbon in a species-free habitat. Then Theorem 2.4 states that the species can persist if the quota $Q^{*}$ exceeds the quota $\hat{Q}$ required for growth to balance losses (i.e., $\mu(\hat{Q}) = D$). Thus, the persistence criterion (i.e., $\mu(Q^{*})-D>0$) summarizes the characteristics of carbon uptake, the growth rate, light availability, and the respiration rate. If the quota $Q^{*}$ is less than the quota $\hat{Q}$, then we can show that the species population is washed out of the habitat (see Theorem 2.5), where we have ignored the effect of respiration (see the specific assumption (2.11)), due to a technical reason.

    In Theorem 2.4, we only show that the single-species model (2.1) admits at least one positive equilibrium if the species can persist by using the theory of uniform persistence. The uniqueness and global stability of positive equilibrium for (2.1) are still open if no extra assumptions are imposed. Thus, two semi-trivial steady-state solutions of the two-species system (1.1), corresponds to the presence of one of the species and the absence of the other species, are not necessarily unique. This makes the investigation of coexistence for the two-species system (1.1) more difficult. Inspired by [12,Section 4], we first define two suitable parameters, $\hat{Q}_2^{\min}$ and $\check{Q}_1^{\min}$ (see 3.3 and 3.4), then we are able to show that the compact attractor $\mathcal{M}_1$ (on the boundary $u_2 = 0$), and the compact attractor $\mathcal{M}_2$ (on the boundary $u_1 = 0$) are uniform weak repellers for two-species system (1.1) (see Lemma 3.3 and Lemma 3.4) under appropriate conditions depending on $\hat{Q}_2^{\min}$ and $\check{Q}_1^{\min}$, respectively. Finally, we are able to show that system (1.1) is uniformly persistent, and it admits at least one coexistence (componentwise positive) steady-state solution (see Theorem 3.5) when the trivial steady-state solution, the compact set $\mathcal{M}_1$, and the compact set $\mathcal{M}_2$ are all invasible. From biological viewpoints, invasibility will depend on whether the missing competitor obtains sufficiently large quotas ($\hat{Q}_2^{\min}$ or $\check{Q}_1^{\min}$) to permit a growth rate that exceeds the loss to dilution ($D$). Then robust coexistence occurs when there is mutual invasibility of both $\mathcal{M}_1$ and $\mathcal{M}_2$.

    Research of F.-B. Wang is supported in part by Ministry of Science and Technology, Taiwan; and National Center for Theoretical Sciences (NCTS), National Taiwan University; and Chang Gung Memorial Hospital (CRRPD3H0011, BMRPD18 and NMRPD5F0543). We would like to express our thanks to Professor Sze-Bi Hsu for the helpful discussions in this paper.

    All authors declare no conflicts of interest in this paper.



    [1] A. Cunningham and R. M. Nisbet, Time lag and co-operativity in the transient growth dynamics of microalgae, J. Theoret. Biol., 84 (1980), 189–203.
    [2] A. Cunningham and R. M. Nisbet, Transient and Oscillation in Continuous Culture, in Mathematics in Microbiology, M. J. Bazin, ed., Academic ress, New York, 1983.
    [3] M. Droop, Vitamin B12 and marine ecology. IV. The kinetics of uptake, growth and inhibition in Monochrysis Lutheri, J. Mar. Biol. Assoc. UK, 48 (1968), 689–733.
    [4] M. Droop, Some thoughts on nutrient limitation in algae, J. Phycol., 9 (1973), 264–272.
    [5] M. Droop, The nutrient status of algal cells in continuous culture, J. Mar. Biol. Assoc. UK, 54 (1974), 825–855.
    [6] J. P. Grover, Constant- and variable-yield models of population growth: Responses to environmental variability and implications for competition, J. Theoret. Biol., 158 (1992), 409–428.
    [7] J. P. Grover, Resource Competition, Chapman and Hall, London, 1997.
    [8] J. P. Grover, Resource storage and competition with spatial and temporal variation in resource availability, The American Naturalist, 178 (2011), E 124–E 148.
    [9] S.-B Hsu and C. J. Lin, Dynamics of two phytoplankton Species Competing for light and nutrient with internal storage, Discrete Cont. Dyn. S, 7 (2014), 1259–1285.
    [10] S. B. Hsu, K. Y. Lam and F. B. Wang, Single species growth consuming inorganic carbon with internal storage in a poorly mixed habitat, J. Math. Biol., 75 (2017), 1775–1825.
    [11] J. Huisman, P. v. Oostveen and F. J.Weissing, Species dynamics in phytoplankton blooms: incomplete mixing and competition for light, The American Naturalist, 154 (1999), 46–67.
    [12] S. B. Hsu, F. B.Wang, and X. Q. Zhao, A reaction-diffusion model of harmful algae and zooplankton in an ecosystem, J. Math. Anal. Appl., 451 (2017), 659–677.
    [13] J. Jiang, On the global stability of cooperative systems, Bull London Math. Soc., 26 (1994), 455– 458.
    [14] J. T. O. Kirk, Light and photosynthesis in aquatic ecosystems, 2nd edition, Cambridge University Press, Cambridge, 1994.
    [15] P. D. Leenheer, S. A. Levin, E. D. Sontag and C. A. Klausmeier, Global stability in a chemostat with multiple nutrients, J. Math. Biol., 52 (2006), 419–438.
    [16] F. M. M. Morel, Kinetics of nutrient uptake and growth in phytoplankton, J. Phycol., 23 (1987), 137–150.
    [17] H. Nie, S. B. Hsu and J. P. Grover, Algal Competition in a water column with excessive dioxide in the atmosphere, J. Math. Biol., 72 (2016), 1845–1892.
    [18] H. L. Smith, Monotone Dynamical Systems:An Introduction to the Theory of Competitive and Cooperative Systems, Math. Surveys Monogr 41, American Mathematical Society Providence, RI, 1995.
    [19] H. R. Thieme, Convergence results and a Poincare-Bendixson trichotomy for asymptotically autonomous differential equations, J. Math. Biol., 30 (1992), 755–763.
    [20] D. B. V. deWaal, J. M. H. Verspagen, J. F. Finke, V. Vournazou, A. K. Immers,W. E. A. Kardinaal, L. Tonk, S. Becker, E. V. Donk, P. M. Visser and J. Huisman, Reversal in competitive dominance of a toxic versus non-toxic cyanobacterium in response to rising CO2, ISME J., 5 (2011), 1438–1450.
    [21] X. Q. Zhao, Dynamical Systems in Population Biology, Springer, New York, 2003.
  • Reader Comments
  • © 2019 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4458) PDF downloads(644) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog