We investigate a novel model of coupled stochastic differential equations modeling the interaction of mussel and algae in a random environment, in which combined effect of white noises and telegraph noises formulated under regime switching are incorporated. We derive sufficient condition of extinction for mussel species. Then with the help of stochastic Lyapunov functions, a well-grounded understanding of the existence of ergodic stationary distribution is obtained. Meticulous numerical examples are also employed to visualize our theoretical results in detail. Our analytical results indicate that dynamic behaviors of the stochastic mussel-algae model are intimately associated with two kinds of random perturbations.
Citation: Yan Xie, Zhijun Liu, Ke Qi, Dongchen Shangguan, Qinglong Wang. A stochastic mussel-algae model under regime switching[J]. Mathematical Biosciences and Engineering, 2022, 19(5): 4794-4811. doi: 10.3934/mbe.2022224
Related Papers:
[1]
Dengxia Zhou, Meng Liu, Ke Qi, Zhijun Liu .
Long-time behaviors of two stochastic mussel-algae models. Mathematical Biosciences and Engineering, 2021, 18(6): 8392-8414.
doi: 10.3934/mbe.2021416
[2]
Hui Jiang, Ling Chen, Fengying Wei, Quanxin Zhu .
Survival analysis and probability density function of switching heroin model. Mathematical Biosciences and Engineering, 2023, 20(7): 13222-13249.
doi: 10.3934/mbe.2023590
[3]
Wenjie Yang, Qianqian Zheng, Jianwei Shen, Linan Guan .
Bifurcation and pattern dynamics in the nutrient-plankton network. Mathematical Biosciences and Engineering, 2023, 20(12): 21337-21358.
doi: 10.3934/mbe.2023944
[4]
Lin Li, Wencai Zhao .
Deterministic and stochastic dynamics of a modified Leslie-Gower prey-predator system with simplified Holling-type Ⅳ scheme. Mathematical Biosciences and Engineering, 2021, 18(3): 2813-2831.
doi: 10.3934/mbe.2021143
[5]
Xuehui Ji, Sanling Yuan, Tonghua Zhang, Huaiping Zhu .
Stochastic modeling of algal bloom dynamics with delayed nutrient recycling. Mathematical Biosciences and Engineering, 2019, 16(1): 1-24.
doi: 10.3934/mbe.2019001
[6]
Xueqing He, Ming Liu, Xiaofeng Xu .
Analysis of stochastic disease including predator-prey model with fear factor and Lévy jump. Mathematical Biosciences and Engineering, 2023, 20(2): 1750-1773.
doi: 10.3934/mbe.2023080
[7]
H. J. Alsakaji, F. A. Rihan, K. Udhayakumar, F. El Ktaibi .
Stochastic tumor-immune interaction model with external treatments and time delays: An optimal control problem. Mathematical Biosciences and Engineering, 2023, 20(11): 19270-19299.
doi: 10.3934/mbe.2023852
[8]
Yuanpei Xia, Weisong Zhou, Zhichun Yang .
Global analysis and optimal harvesting for a hybrid stochastic phytoplankton-zooplankton-fish model with distributed delays. Mathematical Biosciences and Engineering, 2020, 17(5): 6149-6180.
doi: 10.3934/mbe.2020326
[9]
Daniele Cappelletti, Badal Joshi .
Transition graph decomposition for complex balanced reaction networks with non-mass-action kinetics. Mathematical Biosciences and Engineering, 2022, 19(8): 7649-7668.
doi: 10.3934/mbe.2022359
[10]
Shengyu Huang, Hengguo Yu, Chuanjun Dai, Zengling Ma, Qi Wang, Min Zhao .
Dynamic analysis of a modified algae and fish model with aggregation and Allee effect. Mathematical Biosciences and Engineering, 2022, 19(4): 3673-3700.
doi: 10.3934/mbe.2022169
Abstract
We investigate a novel model of coupled stochastic differential equations modeling the interaction of mussel and algae in a random environment, in which combined effect of white noises and telegraph noises formulated under regime switching are incorporated. We derive sufficient condition of extinction for mussel species. Then with the help of stochastic Lyapunov functions, a well-grounded understanding of the existence of ergodic stationary distribution is obtained. Meticulous numerical examples are also employed to visualize our theoretical results in detail. Our analytical results indicate that dynamic behaviors of the stochastic mussel-algae model are intimately associated with two kinds of random perturbations.
1.
Introduction
The invasion of alien species has caused significant economic and ecological problems in a diverse array of ecosystems. One of the most instructive examples was the zebra mussel, a highly invasive bivalve in the Great Lakes in 1986, which has high fecundity and strong ability, not only causing ecological modifications and economic damages, but also bringing great inconvenience to human life [1,2], such as blocking pipes, polluting water sources and crowding out native species. In recent years, governments in many countries have invested lots of capital in mussel control and damage repair [3], and an increasing number of scholars have begun to study this hot issue. Many significant existing literatures indicate that the growth and survival of mussel are mainly dependent on algae, which means that the food supply of algae will limit the intake of mussel (see [4,5,6,7]). To date, quite a few ecologists and applied mathematical researchers have paid close attention to the development of mussel with various types of mathematical mussel-algae models, and numerous results have been presented about mussel-algae models on account of their extensive applications [8,9,10,11,12]. For example, Turing pattern [8], Turing-Hopf bifurcation [9,10], wavetrain solution [11] and so on.
Reviewing many of the above investigations, we can see that deterministic model can simplify complex system, which is conducive to theoretical research and analysis. However, when the external disturbance is large, the population number may be significantly reduced, the law of large numbers is not on trial, and the deterministic model is not objective. If the noises are ignored, the model will be inaccurate. So in order to describe the species in the real world more objectively, the environmental fluctuations should be taken into account[13,14,15,16]. It has been found that environmental stress, salinity and so on have a remarkable effect on the microbiota associated with zebra mussel, therefore influence the survival analysis of mussel [17]. Usually, one of the most common and indispensable environmental fluctuations in an ecosystem, white noises, which can be used to describe these continuous and small fluctuations, have recently become quite attractive (see [18,19,20,21,22]).
However, there are several moderate fluctuations which can not be characterised by white noises [23], let alone some sudden environmental abrupt fluctuations. Take temperature as an example, organisms have a certain average range of optimal growth temperature, and the change of temperature during individuals of the same species or different physiological periods of the same species, may have an effect on population separation and behaviour differences [24]. With the exception of birds and mammals, which are thermostats, all other animals are poikilotherms, namely, their body temperature always varies in response to the changes in external temperature [25]. When the effect of temperature change exceeds the threshold of environment or ecosystem, results can range from mild individual death to species turnover or ecosystem change [25]. There is one vivid example in Finland during the 1950s, the pikeperch population crumbled just because of a trifling cooling together with the grown fishing mortality [26], manifesting that even a slight change of mean temperature, like 0.5−1.0∘C, may appreciably affect populations living nearby [27]. Although the temperature in annual specific season is relatively constant in general, the stochastic effects on temperature still may lead to tremendous variations in larval fish survival [28], and hence none of the white noises or the period describing form may be applicable for this fluctuation. Coincidentally, it has been observed that temporal patterns in the life cycle of the zebra mussel vary on increasing temperature in experiment [29]. Thus, when the system jumps back and forth caused by moderate fluctuations or sudden environmental distributions, a further step to investigate telephone noises, another common noises type, can be used to elucidate a switching between different (two or more) environmental regimes, is noteworthy [30]. Because this switching is memoryless and the waiting time of two different switches obeys an exponential distribution [31], using state switching of a continuous-time Markov chain to drive the changes in the main parameters of system can satisfy our needs. Numerous literatures focus on these issues, for example [32,33,34] and the reference therein.
The necessary environmental fluctuations exist everywhere in various forms. In recent years, as two most common and representative noises, white noises and telephone noises are widely used in stochastic population systems, which has been studied with the efforts of many researchers, see [35,36,37,38,39]. Especially, a stochastic phytoplankton allelopathy model formulated by Zhao et al. [35], which derives the stationary distribution and its statistical characteristics (i.e., the mean and variance). Considering a Holling-II stochastic predator-prey model, Jiang et al. [36] established a set of sufficient criteria for the strong persistence in the mean, extinction as well as stationary distribution. Very recently, Cai et al. [37] dealt with the long-time properties of a stochastic delay foraging arena predator-prey model, including stochastic ultimate boundedness, pathwise estimation and extinction.
As is well known, there are few papers related to stochastic mussel models with white noises and telephone noises. There is a doubt whether or not we can derive that the extinction and the existence of unique stationary distribution of such regime-switching differential system. Fortunately, the answer is affirmative and the related detailed proofs will be given later.
Focusing on two essentials, the extinction of mussel species and ergodic stationary distribution, the paper is bent on exploring the long-time properties of model (2.4). The present paper has the following structure: in Section 2 we introduce a stochastic mussel-algae model simultaneously affected by white noises and telephone noises. Section 3 ensures the existence and positivity of global solution. Some dynamic behaviors about extinction are carried out in Section 4 and ergodic stationary distribution are proved in Section 5, respectively. Several numerical examples are set to demonstrate our mathematical results in Section 6. Section 7 ends with some conclusions and future directions. A necessary theory listed in Appendix. Last but not least, several literatures inspired the present work are listed.
2.
Model and preliminaries
For the final export of the model we will discuss, let us first introduce the following nondiffusive mussel-algae model which corresponds to a nondiffusive version proposed by Koppel et al.[12].
With the idea of the influence of intraspecific competition, the above model (2.1) is recently extended by Zhou et al. [40] to the following deterministic version
where B1(t) and B2(t) are the standard independent Brownian motions, σ21 and σ22 represent their intensities. Thanks to the analysis of the dynamical behaviors of model (2.3), not only the survival of mussel as well as ESD (ergodic stationary distribution) are obtained, but an interesting result is discovered, namely the lower amplitude noises could cause mussel species outbreak, while the higher amplitude noises could make mussel species become extinct.
As an extension of model (2.3), motivated by the above studies we further introduce the telephone noises into model (2.3), and suppose that Markov chain r(t) take values in state space S={1,2,...,N} controlling the switching between the environmental regimes [41]. Then the resulting stochastic mussel-algae model simultaneously affected by white noises and telephone noises can be expressed as
r(0)∈S, X(0)>0, Y(0)>0. Here, suppose that Markov chain r(t) is independent of the Brownian Motions B1(t) and B2(t). The coefficients f(k) are all nonnegative constants for any k∈S, where f=g, b, c, β, l, σ1, σ2, Yup.
For the convenience of subsequent sections, we list the following preliminaries. Define probability space by (Ω,F,{Ft}t≥0,P) with a filtration {Ft}t≥0 satisfying the usual conditions. Assign Rn+={z∈Rn|zi>0,1≤i≤n}. Denote right-continuous Markov chain r(t)(t≥0) on the complete probability space whose domain is a finite state space S. Let ˆh=mink∈S{h(k)} and ˇh=maxk∈S{h(k)}, where h=(h(1),...,h(N)). The generator ˜Γ=(μij)N×N is given, for △>0,μij≥0, by
where μij implies the transition rate from state i to state j, and keeps nonnegative if i≠j satisfying ∑Ni=1μij=0. When r(t) is irreducible, then it admits a unique stationary distribution π={π1,π2,...,πN} such that π˜Γ=0 satisfying ∑Nk=1πk=1(πk>0) for any k∈S.
By the help of r(⋅), define the diffusion process (q(t),r(t)) as below:
where B(⋅) signifies Brownian motion. Assign ς(⋅,⋅):Rn×S→Rn,G(⋅,⋅):Rn×S→Rn×n such that G(q,k)GT(q,k)=(dij(q,k)). For each k∈S, let W(q,k) be twice continuously differentiable with respect to q, and define the operator L as:
(˜X(0),Y(0))=(lnX(0),Y(0)). It's clear that the coefficients of model (3.1) obey locally Lipschitz continuous conditions. Thus, a unique solution (˜X(t),Y(t)) on [0,τe) to model (3.1) which belongs to τe≤+∞ exists. In other words, the presence of a unique local positive solution of model (2.4) is guaranteed on [0,τe), denoted by (X(t),Y(t))=(e˜X(t),Y(t)). The only thing we need to do is to verify τe=+∞a.s. Denote n0>0 be enough large such that 1/n0≤X(0),Y(0)≤n0. For each n≥n0, assign
Set τ∞=limn→+∞τn, then τ∞≤τe. Testify τ∞=+∞ can finish this proof. If this assertion is false, then assume that τe<∞. There are two constants T>0 and ϵ∈(0,1) such that P{τ∞≤T}>ϵ. Consequently there exists an integer n1≥n0 satisfying
P{τn≤T}≥ϵ.
(3.2)
Define
U(X,Y,k)=(X−1−lnX)+ap(Y−1−lnY),
from which we conclude, together with Itô's formula, that
where λ1>0 is sufficiently large. According to Eqs (5.2)–(5.3), we know that
lim infℓ1→+∞,(X,Y,k)∈(R2+∖Uℓ1)×SW(X,Y,k)=+∞,
where Uℓ1=(1ℓ1,ℓ1)×(1ℓ1,ℓ1), ℓ1>1 is a sufficiently large number. Continuous function W(X,Y,k) guarantees that it has a minimum point (X0,Y0,k) in R2+×S, we further define the following C2-function
Summarizing the above discussions, the conditions (i)-(iii) of Lemma A.1 have been checked, which means that model (2.4) has a ESD. The proof is completed.
6.
Numerical simulations
Here, it's time to verify our analytical results by some numerical simulations. Give the initial values of model (2.4) by (X(0),Y(0))=(0.3,0.2). Suppose that r(t) ranges in state space S={1,2} with the generator
Γ=(−664−4),
and the corresponding stationary distribution π=(π1,π2)=(0.4,0.6).
Example 6.1 We fix a=0.2, p=0.1. When r(t)=1, we choose a set of parameters with c(1)=0.02, g(1)=0.01, b(1)=0.1, Yup(1)=1, l(1)=0.4, β(1)=150, σ1(1)=0.3, σ2(1)=0.3. When r(t)=2, we select the other group parameters with c(2)=0.3, g(2)=0.02, b(2)=0.2, Yup(2)=1.5, l(2)=0.6, β(2)=200, σ1(2)=0.4, σ2(2)=0.4. By calculating, we get ∑Nk=1πk(aˇbˆll(k)Yup(k)−12σ21(k))=−0.005<0, which together with Theorem 4.1 shows that mussel X(t) of model (2.4) tends to be extinct, depicted in Figure 1.
Figure 1.
a. The sample paths of the Markov chain r(t). b. The mussel population density X(t).
Example 6.2 Fix a=0.7, p=0.1. When r(t)=1, select a set of parameters with c(1)=0.03, g(1)=0.02, b(1)=0.1, Yup(1)=1.2, l(1)=0.6, β(1)=150, σ1(1)=0.05, σ2(1)=0.05. And choose the other set of parameters with c(2)=0.04, g(2)=0.01, b(2)=0.2, Yup(2)=2.2, l(2)=0.7, β(2)=200, σ1(2)=0.1, σ2(2)=0.1 when r(t)=2. We then compute R0≈6.18>1 and ˇσ22−ˆl=−0.59<0. As a result, we know from Theorem 5.1 that model (2.4) has a unique ESD. Simulations shown in Figure 2 can confirm this clearly.
Figure 2.
a. The sample paths of the Markov chain r(t), the mussel population density X(t) and algae population density Y(t), respectively. b. The probability density function of r(t), X(t) and Y(t), respectively.
Based on the previous work [40], this article focuses on the long-time behaviors of a stochastic mussel-algae model incorporating both white noises and telephone noises. To the best of our knowledge, it is the first attempt to analyze the above stochastic mussel-algae model with two types of environmental noises. We establish two sets of sufficient criteria to verify the extinction and the existence of a unique ESD. And some crucial influences of environmental noises on mussel evolution are revealed. Based on the above theoretical analysis, this paper draws the followings conclusions.
● We can conclude, by Theorem 4.1, that the mussel will be extinct if
After further observation, if the conversion rate of ingested algae to mussel production a and the rate of consumption of the algae by a mussel b(k) are suitably small while the amplitude noises σ1(k) is relatively large, then the mussel species die out (see Figure 1), which means environmental fluctuations can be used as a strategy to control the growth of mussel.
● We have demonstrated, according to Theorem 5.1, that model (2.4) has a unique ESD if R0>1 and ˆl>ˇσ22 which can be used to estimate the probability of mussel outbreaks. This fact implies that large conversion rate of ingested algae to mussel production a, large rate of consumption of the algae by a mussel b(k), while small maximal per capita mussel death rate c(k) and small noises amplitude σi(k)(i=1,2) are advantage for the persistence of species mussel (see Figure 2).
Looking forward, some challenging topics deserve further research studies. We may investigate more complex model by introducing different noises perturbations into mussel-algae model (2.4), such as Ornstein-Uhlebeck process, seim-Markov, Lévy jump[43]. Taking into account the long-term memory of population systems, further consideration of fractional-order modeling methods has important theoretical and practical implications[44,45,46,47,48]. We will leave these tasks in future.
Acknowledgements
We would like to thank the referees and editor for their careful reading of the original manuscript and their valuable comments and suggestions that greatly improved the presentation of this work. The work is supported by the NNSFs of China (Nos. 11871201, 12101211).
Conflict of interest
The authors declare there is no conflict of interest.
Appendix
The necessary theory listed as below, which is fundamental for our analysis result.
Lemma A.1System (1.5) is ergodic and positive recurrent if the following assumptions hold ([42])
(i) μij>0 for any i≠j;
(ii) given k∈S, P(z,k)=(pij(z,k)) is symmetric satisfying
ϱ|ξ|2≤ξTP(q,i)ξ≤ϱ−1|ξ|2,ξ∈Rn,
with some constant ϱ∈(0,1];
(iii) there exists a non-negative function W(⋅,k):Dc→R, for each k∈S, such that W(⋅,k) is twice continuously differentable and for some α>0,
LW(q,k)≤−α,forany(q,k)∈Dc×S,
where Dc is the complement of a bounded open subset D∈Rn with a smooth boundary. Moreover, for any Borel measurable function ς(⋅,⋅):Rn×S→R, there is a unique stationary density π(⋅,⋅) such that
H. J. MacIsaac, Potential abiotic and biotic impacts of zebra mussels on the inland waters of North America, Am. Zool., 36 (1996), 287–299. https://doi.org/10.1093/icb/36.3.287 doi: 10.1093/icb/36.3.287
[2]
P. Dolmer, Algal concentration profiles above mussel beds, J. Sea Res., 43 (2000), 113–119. https: //doi.org/10.1016/S1385-1101(00)00005-8 doi: 10.1016/S1385-1101(00)00005-8
[3]
Q. H. Huang, H. Wang, M. A. Lewis, A hybrid continuous/discrete-time model for invasion dynamics of zebra mussels in rivers, SIAM J. Appl. Math., 77 (2017), 854–880. https://doi.org/10.1137/16M1057826 doi: 10.1137/16M1057826
[4]
D. J. Wildish, D. D. Kristmanson, Importance to mussels of the benthic boundary layer, Can. J. Fish. Aquat. Sci., 41 (1984), 1618–1625. https://doi.org/10.1139/f84-200 doi: 10.1139/f84-200
[5]
G. Øie, K. I. Reitan, O. Vadstein, H. Reinertsen, Effect of nutrient supply on growth of blue mussels (Mytilus edulis) in a landlocked bay, Hydrobiologia, 484 (2002), 99–109. https://doi.org/10.1007/978-94-017-3190-4_9 doi: 10.1007/978-94-017-3190-4_9
[6]
J. Widdows, J. S. Lucas, M. D. Brinsley, P. N. Salkeld, F. J. Staff, Investigation of the effects of current velocity on mussel feeding and mussel bed stability using an annular flume, Helgol. Mar. Res., 56 (2002), 3–12. https://doi.org/10.1007/s10152-001-0100-0 doi: 10.1007/s10152-001-0100-0
[7]
H. M. Page, D. M. Hubbard, Temporal and spatial patterns of growth in mussels Mytilus edulis on a offshore platform: relationships to water temperature and food availability, J. Exp. Mar. Biol. Ecol., 111 (1987), 159–179. https://doi.org/10.1016/0022-0981(87)90053-0 doi: 10.1016/0022-0981(87)90053-0
[8]
R. A. Cangelosi, D. J. Wollkind, B. J. Kealy-Dichone, I. Chaiya, Nonlinear stability analyses of Turing patterns for a mussel-algae model, J. Math. Biol., 70 (2015), 1249–1294. https://doi.org/10.1007/s00285-014-0794-7 doi: 10.1007/s00285-014-0794-7
[9]
Z. L. Shen, J. J. Wei, Spatiotemporal patterns in a delayed reaction-diffusion mussel-algae model, Int. J. Bifur. Chaos. Appl. Sci. Engrg., 29 (2019), 1950164. https://doi.org/10.1142/S0218127419501645 doi: 10.1142/S0218127419501645
[10]
S. Djilali, B. Ghanbari, S. Bentout, A. Mezouaghi, Turing-Hopf bifurcation in a diffusive mussel-algae model with time-fractional-order derivative, Chaos Solitons Fractals, 138 (2020), 109954. https://doi.org/10.1016/j.chaos.2020.109954 doi: 10.1016/j.chaos.2020.109954
[11]
M. Holzer, N. Popović, Wavetrain solutions of a reaction-diffusion-advection model of musselalgae interaction, SIAM. J. Appl. Dyn. Syst., 16 (2017), 431–478. https://doi.org/10.1137/15M1040463 doi: 10.1137/15M1040463
[12]
J. V. D. Koppel, M. Rietkerk, N. Dankers, P. M. J. Herman, Scale-dependent feedback and regular spatial patterns in young mussel beds, Am. Nat., 165 (2005), E66–E77. https://doi.org/10.1086/428362 doi: 10.1086/428362
[13]
L. Arnold, W. Horsthemke, J. W. Stucki, The influence of external real and white noise on the Lotka-Volterra model, Biom. J., 21 (1979), 451–471. https://doi.org/10.1002/bimj.4710210507 doi: 10.1002/bimj.4710210507
[14]
M. Liu, Optimal harvesting of stochastic population models with periodic coefficients, J. Nonlinear Sci., 32 (2022), 1–14. https://doi.org/10.1007/s00332-021-09758-6 doi: 10.1007/s00332-021-09758-6
[15]
X. W. Yu, S. L. Yuan, T. H. Zhang, Survival and ergodicity of a stochastic phytoplankton-zooplankton model with toxin-producing phytoplankton in an impulsive polluted environment, Appl. Math. Comput., 347 (2019), 249–264. https://doi.org/10.1016/j.amc.2018.11.005 doi: 10.1016/j.amc.2018.11.005
[16]
X. W. Yu, S. L. Yuan, T. H. Zhang, Asymptotic properties of stochastic nutrient-plankton food chain models with nutrient recycling, Nonlinear Anal. Hybrid Syst., 34 (2019), 209–225. https://doi.org/10.1016/j.nahs.2019.06.005 doi: 10.1016/j.nahs.2019.06.005
[17]
P. P. Mathai, J. H. Bertram, S. K. Padhi, V. Singh, I. E. Tolo, A. Primus, et al., Influence of environmental stressors on the microbiota of zebra mussels (dreissena polymorpha), Microb. Ecol., 81 (2021), 1042–1053. https://doi.org/10.1007/s00248-020-01642-2 doi: 10.1007/s00248-020-01642-2
J. Hu, Q. M. Zhang, A. Meyer-Baese, M. Ye, Stationary distribution of a stochastic Alzheimer's disease model, Math. Methods Appl. Sci., 43 (2020), 9706–9718. https://doi.org/10.1002/mma.6642 doi: 10.1002/mma.6642
[20]
M. Y. Song, W. J. Zuo, D. Q. Jiang, T. Hayat, Stationary distribution and ergodicity of a stochastic cholera model with multiple pathways of transmission, J. Franklin Inst., 357 (2020), 10773–10798. https://doi.org/10.1016/j.jfranklin.2020.04.061 doi: 10.1016/j.jfranklin.2020.04.061
[21]
H. K. Qi, X. Z. Meng, Threshold behavior of a stochastic predator-prey system with prey refuge and fear effect, Appl. Math. Lett., 113 (2021), 106846. https://doi.org/10.1016/j.aml.2020.106846 doi: 10.1016/j.aml.2020.106846
[22]
X. R. Mao, S. Sabanis, E. Renshaw, Asymptotic behaviour of the stochastic Lotka-Volterra model, J. Math. Anal. Appl., 287 (2003), 141–156. https://doi.org/10.1016/S0022-247X(03)00539-0 doi: 10.1016/S0022-247X(03)00539-0
[23]
M. Liu, M. L. Deng, Analysis of a stochastic hybrid population model with Allee effect, Appl. Math. Comput., 364 (2020), 124582. https://doi.org/10.1016/j.amc.2019.124582 doi: 10.1016/j.amc.2019.124582
[24]
L. J. Ouyang, X. W. Guo, The effection of temperature on the feeding and growth of fish, Stud. Mar. Sin., 49 (2009), 87–95.
[25]
X. Pang, Research on the Fish Eco-physiological (Thermal Tolerance, Metabolism and Swimming) Response Based on Temperature Changes, Ph.D thesis, College of Resources and Environmental Science of Chongqing University, 2020.
[26]
P. J. Colby, H. Lehtonen, Suggested causes for the collapse of zander, Stizostedion lucioperca (L.), populations in northern and central Finland through comparisons with North American walleye, Stizostedion vitreum (Mitchill), Aqua. Fenn., 24 (1994), 9–20.
[27]
H. Lehtonen, Potential effects of global warming on northern European freshwater fish and fisheries, Fish. Manage. Ecol., 3 (1996), 59–71. https://doi.org/10.1111/j.1365-2400.1996.tb00130.x doi: 10.1111/j.1365-2400.1996.tb00130.x
[28]
J. Lappalainen, Effects of Environmental Factors, Especially Temperature, on the Population Dynamics of Pikeperch (Stizostedion lucioperca (L.)), Ph.D. thesis, University of Helsinki, Finland Academic, 2001.
[29]
E. M. Griebeler, A. Seitz, Effects of increasing temperatures on population dynamics of the zebra mussel dreissena polymorpha: implications from an individual-based model, Oecologia, 151 (2007), 530–543. https://doi.org/10.1007/s00442-006-0591-0 doi: 10.1007/s00442-006-0591-0
[30]
M. Slatkin, The dynamics of a population in a Markovian environment, Ecology, 59 (1978), 249–256. https://doi.org/10.2307/1936370 doi: 10.2307/1936370
[31]
A. Settati, A. Lahrouz, Stationary distribution of stochastic population systems under regime switching, Appl. Math. Comput., 244 (2014), 235–243. https://doi.org/10.1016/j.amc.2014.07.012 doi: 10.1016/j.amc.2014.07.012
[32]
R. McVinish, P. K. Pollett, Y. S. Chan, A metapopulation model with Markovian landscape dynamics, Theor. Popul. Biol., 112 (2016), 80–96. https://doi.org/10.1016/j.tpb.2016.08.005 doi: 10.1016/j.tpb.2016.08.005
[33]
D. Li, S. Q. Liu, Threshold dynamics and ergodicity of an SIRS epidemic model with Markovian switching, J. Differ. Equ., 263 (2017), 8873–8915. https://doi.org/10.1016/j.jde.2017.08.066 doi: 10.1016/j.jde.2017.08.066
[34]
Q. Liu, D. Q. Jiang, T. Hayat, A. Alsaedi, Dynamical behavior of a multigroup SIRS epidemic model with standard incidence rates and Markovian switching, Discrete Contin. Dyn. Syst. Ser. A., 39 (2019), 5683. https://doi.org/10.3934/dcds.2019249 doi: 10.3934/dcds.2019249
[35]
Y. Zhao, S. L. Yuan, T. H. Zhang, The stationary distribution and ergodicity of a stochastic phytoplankton allelopathy model under regime switching, Commun. Nonlinear Sci. Numer. Simul., 37 (2016), 131–142. https://doi.org/10.1016/j.cnsns.2016.01.013 doi: 10.1016/j.cnsns.2016.01.013
[36]
X. B. Jiang, L. Zu, D. Q. Jiang, D. O'Regan, Analysis of a stochastic Holling type II predator-prey model under regime switching, Bull. Malays. Math. Sci. Soc., 43 (2020), 2171–2197. https://doi.org/10.1007/s40840-019-00798-6 doi: 10.1007/s40840-019-00798-6
[37]
Y. M. Cai, S. Y. Cai, X. R. Mao, Stochastic delay foraging arena predator-prey system with Markov switching, Stoch. Anal. Appl., 38 (2020), 191–212. https://doi.org/10.1080/07362994.2019.1679645 doi: 10.1080/07362994.2019.1679645
[38]
Y. X. Zhou, W. J. Zuo, D. Q. Jiang, M. Y. Song, Stationary distribution and extinction of a stochastic model of syphilis transmission in an MSM population with telegraph noises, J. Appl. Math. Comput., 66 (2021), 645–672. https://doi.org/10.1007/s12190-020-01453-1 doi: 10.1007/s12190-020-01453-1
[39]
C. Q. Xu, S. L. Yuan, T. H. Zhang, Average break-even concentration in a simple chemostat model with telegraph noise, Nonlinear Anal. Hybrid Syst., 29 (2018), 373–382. https://doi.org/10.1016/j.nahs.2018.03.007 doi: 10.1016/j.nahs.2018.03.007
[40]
D. X. Zhou, M. Liu, K. Qi, Z. J. Liu, Long-time behaviors of two stochastic mussel-algae models, Math. Biosci. Eng., 18 (2021), 8392–8414. https://doi.org/10.3934/mbe.2021416 doi: 10.3934/mbe.2021416
[41]
H. Liu, X. X. Li, Q. S. Yang, The ergodic property and positive recurrence of a multi-group Lotka-Volterra mutualistic system with regime switching, Syst. Control Lett., 62 (2013), 805–810. https://doi.org/10.1016/j.sysconle.2013.06.002 doi: 10.1016/j.sysconle.2013.06.002
[42]
X. R. Mao, C. G. Yuan, Stochastic Differential Equations with Markovian Switching, Imperial College Press, London, 2006.
[43]
T. T. Ma, X. Z. Meng, Z. B. Chang, Dynamics and optimal harvesting control for a stochastic one-predator-two-prey time delay system with jumps, Complexity, 2019 (2019), 5342031. https://doi.org/10.1155/2019/5342031 doi: 10.1155/2019/5342031
[44]
M. Yavuz, N. Sene, Stability analysis and numerical computation of the fractional predator-prey model with the harvesting rate, Fractal. Fract., 4 (2021), 35. https://doi.org/10.3390/fractalfract4030035 doi: 10.3390/fractalfract4030035
[45]
P. A. Naik, K. M. Owolabi, M. Yavuz, J. Zu, Chaotic dynamics of a fractional order HIV-1 model involving AIDS-related cancer cells, Chaos Solitons Fractals, 140 (2020), 110272. https://doi.org/10.1016/j.chaos.2020.110272 doi: 10.1016/j.chaos.2020.110272
[46]
P. A. Naik, M. Yavuz, S. Qureshi, J. Zu, Modeling and analysis of COVID-19 epidemics with treatment in fractional derivatives using real data from Pakistan, Eur. Phys. J. Plus., 135 (2020), 1–42. https://doi.org/10.1140/epjp/s13360-020-00819-5 doi: 10.1140/epjp/s13360-020-00819-5
[47]
S. Allegretti, I. M. Bulai, R. Marino, M. A. Menandro, K. Parisi, Vaccination effect conjoint to fraction of avoided contacts for a Sars-Cov-2 mathematical model, Math. Model. Num. Simul., 1 (2021), 56–66. https://doi.org/10.53391/mmnsa.2021.01.006 doi: 10.53391/mmnsa.2021.01.006
[48]
P. Kumar, V. S. Erturk, Dynamics of cholera disease by using two recent fractional numerical methods, Math. Model. Num. Simul., 1 (2021), 102–111. https://doi.org/10.53391/mmnsa.2021.01.010 doi: 10.53391/mmnsa.2021.01.010
This article has been cited by:
1.
Shuqi Zhai, Qinglong Wang, Ting Yu,
Fuzzy optimal harvesting of a prey-predator model in the presence of toxicity with prey refuge under imprecise parameters,
2022,
19,
1551-0018,
11983,
10.3934/mbe.2022558
2.
Peng Zhu, Min Xiao, Xia Huang, Fuchen Zhang, Zhen Wang, Jinde Cao,
Spatiotemporal dynamics optimization of a delayed reaction–diffusion mussel–algae model based on PD control strategy,
2023,
173,
09600779,
113751,
10.1016/j.chaos.2023.113751
Yan Xie, Zhijun Liu, Ke Qi, Dongchen Shangguan, Qinglong Wang. A stochastic mussel-algae model under regime switching[J]. Mathematical Biosciences and Engineering, 2022, 19(5): 4794-4811. doi: 10.3934/mbe.2022224
Yan Xie, Zhijun Liu, Ke Qi, Dongchen Shangguan, Qinglong Wang. A stochastic mussel-algae model under regime switching[J]. Mathematical Biosciences and Engineering, 2022, 19(5): 4794-4811. doi: 10.3934/mbe.2022224
The conversion rate of ingested algae to mussel production
b
The rate of consumption of the algae by a mussel
c
The maximal per capita mussel death rate
Yup
The concentration of algae in the upper water layer
l
The rate of exchange between the lower and upper water layers
p
The height of the lower water layer
β
The value of X(t) at which mortality is half maximal
Figure 1. a. The sample paths of the Markov chain r(t). b. The mussel population density X(t)
Figure 2. a. The sample paths of the Markov chain r(t), the mussel population density X(t) and algae population density Y(t), respectively. b. The probability density function of r(t), X(t) and Y(t), respectively