The stochastic dynamics of a population growth model with size-dependent carrying capacity is considered. The effect of a fluctuating environment on population growth is modeled as a multiplicative dichotomous noise. At intermediate values of population size the deterministic counterpart of the model behaves similarly to the Von Foerster model for human population, but at small and very large values of population size substantial differences occur. In the stochastic case, an exact analytical solution for the stationary probability distribution is found. It is established that variation of noise correlation time can cause noise-induced transitions between three different states of the system characterized by qualitatively different behaviors of the probability distributions of the population size. Also, it is shown that, in some regions of the system parameters, variation of the amplitude of environmental fluctuations can induce single unidirectional abrupt transitions of the mean population size.
1. Introduction
The dynamical stability of populations and ecosystems is known to govern their responsiveness to fluctuating environmental conditions, determining with what reliability natural resources provide life-sustaining services to society. Therefore, population and ecosystem dynamics is a major structuring theme in ecology [1]. The most productive abstraction of a self-regulation mechanism for population growth is the Malthus-Verhulst (logistic) model with a constant carrying capacity, which (or its various modifications like Richards model, Gompertz model, etc.) describes well a large class of populations [2, 3]. However, for some populations, such as the human one, empirical data show that the logistic-type models are irrelevant to describe the growth process [4–6]. For example, the historical estimates for the past 2000 years of human population growth indicate that human numbers have been growing faster than exponentially for most of the known history [4, 5]. In 1960 Von Foerster et al. suggested that the classical logistic growth in the case of the human population should be changed, because the individuals in the system can form coalitions that, due to technological progress, can grow faster [4]. In the light of the empirical data and the positive feedback mechanism between innovation and population size, in [4–7] the dynamics of human community was modeled by a growth model with a size-dependent “carrying capacity” K (called Von Foerster model), assuming a simple power-law relationship K~x1+α, where α>0 and x is the population size. Although this model fits the population data for humans during the past 2000 years quite well, in an evolutionary time scale it seems, however, that for the growth of an initially small human population a logistic model is more appropriate. Moreover, as the growth of any expanding population must eventually be limited by a shortage of resources and habitats, for very large values of the population size the Von Foerster model should be modified. Ecological investigations indicate that population dynamics is sensitive to fluctuations of environmental parameters (noise) [8, 9]. Noise, via its interaction with nonlinearity in physical, chemical, and biological systems, has given rise to several counterintuitive phenomena: stochastic resonance [10–12], noise-enhanced stability [13, 14], stochastic transport in ratchets [15, 16], noise-induced phase transitions [17–20], noise-induced transitions [21, 22], memory generated resonance [23–25], and so forth. In particular, the authors of [26] have shown that, in some symbiotic ecosystems, colored fluctuations of the carrying capacities of populations can induce bistability and produce abrupt changes between stable states. The fact that external multiplicative noise can induce multistability as well as discontinuous transitions in some complex systems (see, e.g., [27, 28]) inspired us to apply an analogous approach to analyze single population models. It is of interest, both from theoretical and practical viewpoints, to know whether environmental-noise-induced transitions and abrupt changes of mean population sizes also occur in population growth models with a size-dependent carrying capacity.
Thus motivated, we consider a stochastic population growth model with a size-dependent carrying capacity K. For the abovementioned reasons, the dependence of K on the population size x is assumed such that for small values of x the model behaves similarly to a logistic model; for intermediate and large values of x it behaves like Von Foerster model, and for very large values of population size the carrying capacity tends to a finite constant value. The influence of fluctuating environmental parameters on the growth of the population is modeled as colored fluctuations of the carrying capacity. For the sake of mathematical simplicity it is practicable to use such types of colored noise that enable evaluation of stationary probability density of the population size exactly for any value of the noise correlation time. The corresponding analytical solution would be of major interest in itself, since it can serve as a standard situation for comparison with approximate solutions and/or numerical simulations in the more realistic cases of bounded colored noises. The simplest noise of this kind is the symmetric dichotomous Markovian noise, also known as the random telegraph signal [29]. Thus motivated, in this work the fluctuations of the carrying capacity are modeled as dichotomous Markovian noise.
The main purpose of this paper is to provide exact formulas for the analytical treatment of the stationary probability density of the population size. Moreover, we will demonstrate the presence of noise-induced transitions of the probability distribution and single unidirectional abrupt transitions for the mean population size in the considered population growth model.
To avoid misunderstanding, let us mention that we use the term “noise-induced transitions” in the sense of [17, 22], meaning that this corresponds to a qualitative change in the probability law characterizing the random variable in response to noise parameters. The number and position of the extrema of the stationary probability density are, from this point of view, merely a practical way to monitor such a qualitative change.
2. Deterministic Population Growth Model
As was mentioned in the introduction, the present model is based on a generalization of the logistic equation for population growth (cf. [5, 7])
(1)dx(t)dt=rx[K(x)-x(t)],
where x(t) is the population size at time t and the constant r is the growth rate parameter. Note that in the case of a constant carrying capacity K this equation is the same as the usual logistic equation [2]
(2)dx(t)dt=r~x(1-x(t)K),r~=rK.
In this work the carrying capacity K(x) is assumed to be of the form
(3)K(x)=K0+ϵx1+α1+ϵγx1+α,
where K0 (henceforth, natural carrying capacity) is the carrying capacity without any influence of the population on the value of K. For example, for the human population the exponent α>0 and the parameter ϵ>0 characterize an increase of the carrying capacity due to the technological progress such as the use of tools and fire, the development of agriculture, and the introduction of fossil fuels and fertilizers as well as expansion into new habitats and the removal of limiting factors by the development of vaccines, pesticides, antibiotics, and so forth [5, 6]. The saturation parameter γ>0 takes into account that the growth of any expanding population must eventually be limited by a shortage of resources and habitats.
In the rest of the paper, we will always assume that the parameters α, ϵ, K0, and γ are restricted to the inequalities
(4)0<α<1,γ≪ϵ1/α,γ≪K0-1.
These conditions seem to be relevant for the human population as we bear in mind several datasets expressing the development of mankind on Earth in terms of size and economic impact represented in [5]. For the sake of simplicity, we assume that α is of the form α=m/n, where n>m>0 are integers and m/n is irreducible (i.e., not equal to some other k/s, where k<m and s<n are integers).
By applying the change of variables
(5)y=x1/n,
we rewrite (1) with (3) as
(6)dydt=rγϵyn(1+γϵyn+m)F(y),y≥0,
where
(7)F(y)=-y2n+m+(1γ+K0)yn+m-1γϵyn+K0γϵ.
Depending on the natural carrying capacity K0 the solution of (6) has either four or two equilibria. The unstable trivial equilibrium, y=0, occurs for all biologically possible values of the system parameters.
If the natural carrying capacity is below a critical value K0c, there are, in addition to the trivial equilibrium, two stable and one unstable equilibria, but if K0>K0c, only one nontrivial stable equilibrium occurs. This fact is illustrated in Figure 1 with dependencies of the polynomial F(y) on variable y at various values of K0. By conditions (4) the critical carrying capacity K0c is approximately given by the equation
(8)K0c≈αϵ1/α(11+α)1+1/α.
It is notable that in the case K0<K0c an explosive increase of the population size to a very large limit value (xmax~1/γ) is not possible if the initial population size x(0) is lower than the nontrivial unstable equilibrium value. In this case x saturates at the smaller stable equilibrium point, which is of the same size order as the natural carrying capacity K0. This circumstance and (8) give an ecological interpretation of the parameter ϵ. The quantity ϵ-1/α is a measure of the typical value of the population size above which an explosive increase of the population size is possible. This interpretation is in accordance with intuition: a significant growth of a population due to technological progress is possible only if the population is large enough.
Sketch of the typical dependencies of the polynomial F versus y by conditions (4) at various values of the natural carrying capacity K0 (see also (7)). Solid line, K0<K0c; dashed line, K0=K0c; dotted line, K0>K0c.
3. Stochastic Model
Random interaction with the environment (climate, diseases, etc.) is taken into account by introducing a colored noise in K(x). From now on we will use fluctuations of the natural carrying capacity
(9)K0=K~0(1+aZ(t)),
where colored noise Z(t) is assumed to be a dichotomous Markovian stochastic process [29]. The dichotomous process is a random stationary Markovian process consisting of jumps between two values z=-1,1. The jumps follow in time according to a Poisson process, while the values occur with the stationary probability 1/2. The mean value of Z(t) and the correlation function are
(10)〈Z(t)〉=0,〈Z(t)Z(t′)〉=e-ν|t-t′|,
where the switching rate ν is the reciprocal of the noise correlation time ν=1/τc. Obviously, model (1) with (3) and (9) is biologically meaningful only if
(11)a<1;
that is, the natural carrying capacity is nonnegative. To proceed further with the analytical examination of the stochastic version of model (1) it is reasonable to introduce the following notations:
(12)K01=K~0(1+a),K02=K~0(1-a),Fj(y)=-y2n+m+(1γ+K0j)yn+m-1γϵyn+K0jγϵ,j=1,2.
Now the stationary behavior of the stochastic system described with (6), (7), and (9) can be analytically studied by means of the corresponding composite master equation
(13)∂Pj(y,t)∂t=-∂∂y[rγϵyn(1+γϵyn+m)Fj(y)Pj(y,t)]+∑k=12UjkPk(y,t),
where Pj(y,t) denotes the probability density for the combined process (y,zj,t); j,k=1,2; -z2=z1=1; and
(14)(Ujk)=ν(-121212-12).
The stationary probability density in the y space, P(y), is then evaluated via the stationary probability densities Pjs(y) for the states (y,zj):
(15)P(y)=P1s(y)+P2s(y).
To establish the boundaries of the stationary process y(t) we consider the deterministic equation (6), where the natural carrying capacity K0 is replaced by K0j. Two cases should be discerned (cf. also Figure 1): (i) for K01>K0c (see (8)) there is just one positive solution of the polynomial equation
(16)F1(y)=0,
which is the stable fixed point y=y1 of the deterministic equation (6) at the noise value z=1. Although at the noise value z=-1 (6) may have, depending on the value of K02, two (or one) stable fixed points, in the long time limit (i.e., in a stationary regime) because of the switching of noise values the greater fixed point y=y2 is eventually realized (cf. also Figure 2(a)). Thus, in this case all trajectories y(t) in a stationary regime satisfy the following inequalities:
(17)y1≥y(t)≥y2.
(ii) In the case of K01<K0c the polynomial equation (16) has three positive solutions of y11<y12<y13, of which y11 and y13 are the stable fixed points of (6) with z=1. Obviously, in the noise phase z=-1 (6) also has two stable fixed points, y=y21 and y=y23. Therefore, depending on an initial value y(0), we can single out the following alternatives for a stationary regime (cf. also Figure 2(b)).
A schematic plot of the dependence of the polynomials F1 and F2 on y and attracting intervals. Because of the influence of the noise Z(t) in (6), (7), and (9) the dependence of the function F(y) on y fluctuates with a switching rate ν between two forms F1(y) and F2(y), which correspond to noise values z=1 and z=-1 and are illustrated by curves 1 and 2, respectively. The bold lines depict the attracting intervals. (a) The case of K01>K0c; the stable fixed points of (6) are denoted by y1 and y2. (b) The case of K01<K0c; the fixed points for (6) by noise states z=1 and z=-1 are y1i and y2i, i=1,2,3, respectively.
(1) If y(0)<y12, then the trajectories y(t) are restricted in the interval
(18)y11≥y(t)≥y21.(2) For y(0)>y22 it follows that
(19)y13≥y(t)≥y23.(3) In the case of y12<y(0)<y22 both attracting intervals (18) and (19) are possible; that is, with probabilities w and 1-w we can find the trajectories y(t) in the intervals (18) and (19), respectively.
4. Stationary Probability Distribution
For a stationary state we can solve (13), taking as the boundary condition that there is no probability current at the boundary (17) (for the case of (18) and (19), the results are analogous). After quite simple calculations one can find
(20)P(y)=C(1+γϵyn+m)2yF1(y)F2(y)×exp{-νnrγϵ∫y2yF(ξ)(1+ϵγξn+m)ξF1(ξ)F2(ξ)dξ},
where the constant C is determined by the normalization condition
(21)∫y2y1P(y)dy=1.
Taking into account the transformation (5), it is easy to establish that for the original variable x the stationary probability density P~(x) can be written as
(22)P~(x)=y(x)nxP(y(x)).
In the following, we assume that all 2n+m zeros of the polynomials Fj(y) are distinct; that is, we assume that K0j≠K0c. In this case it is convenient to use the partial fraction expansion
(23)1yFj(y)=1yFj(0)+[Fj′(yj)]-11(y-yj)yj+∑k=12n+m-1[Fj′(yj(k))]-11(y-yj(k))yj(k),
where Fj′(y)≡dFj(y)/dy and yj(k) are zeros of the polynomial Fj(y), which are distinct from yj (yj is also a zero of Fj(y)). By denoting
(24)Gj(y):=∑k=12n+m-1[Fj′(yj(k))]-11(y-yj(k))yj(k),
one can get
(25)P~(x)=C~M(x)(1+γϵx1+α)2x1+ν/rK~0(1-a2)|x1/n-y1|1-A1|x1/n-y2|1-A2,iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiix∈(y2n,y1n),
where
(26)M(x)=∏j=12{1∏k=12n+m-1(x1/n-yj(k))×exp[ϵγ(ξn+m-1-yjn+m-1)(ξ-yj)Fj′(yj))dξ-νn2rγϵ×∫y2x1/n(ϵγ(ξn+m-1-yjn+m-1)(ξ-yj)Fj′(yj)Gj(ξ)(1+ϵγξn+m)+ϵγ(ξn+m-1-yjn+m-1)(ξ-yj)Fj′(yj))dξ]},(27)Aj=νn(1+ϵγyjn+m)2rγϵyj|Fj′(yj)|,j=1,2,
and C~ is the normalization constant. It should be noted that the function M(x) is bounded in the interval x∈[y2n,y1n] including also the boundary points y2n and y1n. Formulas (25)–(27) are also applicable for the attracting intervals represented by inequalities (18) and (19). In this case y2 and y1 in (25)–(27) should be replaced by y21 and y11 (or by y23 and y13 for (19)), respectively.
In Figure 3, the probability density P~(x) is plotted as a function of the population size for various values of the noise switching rate ν. An important observation here is that three types of curves are distinguishable in Figure 3, namely, U-shaped ones, bell-shaped ones, and curves with either two or no local extrema. We interpret these three qualitatively different shapes of P~(x) as different states of the stationary system (see also [17, 22]). From (25) and (27) it follows that the critical noise switching rates νcj, j=1,2, at which noise-induced transitions occur can be expressed as
(28)νcj=2rγϵyj|Fj′(yj)|n(1+γϵyjn+m),j=1,2.
The stationary probability density P~(x) versus the population size x, computed from (25)–(27). The system parameter values are K~0=100, ϵ=0.03, γ=10-5, and α=0.5. (a) The case of K01<K0c with (18), a=0.4. Solid line, ν=45; dashed line, ν=55; dotted line, ν=150. (b) The case of K01>K0c, a=0.8. Solid line, ν=73000; dashed line, ν=74000; dotted line, ν=75000.
As the noise switching rate ν increases (or noise correlation time τc=1/ν decreases) the system undergoes two transitions: for ν<νc2, the probability density exhibits a U-shaped form; if νc2<ν<νc1, then P~(x) is characterized with two or no local extrema; finally, in the case of ν>νc1, P~(x) approaches a bell-shaped form. So the value of noise correlation time has a crucial influence on the position of the most probable population size in the attracting intervals. It should be noted that although in the case of 0<ν<νc1 the probability density either at one or both boundaries (x1=y1n and x2=y2n) approaches infinity (see Figure 3), the corresponding singularities are integrable; that is, the probability distribution P~(x) can be normalized. This follows from (25) and (27) using the fact that the quantities A1 and A2 are positive.
5. Single Unidirectional Transitions
In this section we consider the dependence of the mean population size 〈X〉 on the amplitude of environmental fluctuations. It is important to note that in a stationary regime the value of
(29)〈X〉=∫x2(a)x1(a)xP~(x)dx
depends strongly on the boundary points x1(a) and x2(a) for an attracting interval. If K~0<K0c<2K~0 and the initial population size x(0) is sufficiently small (x(0)≤K~0), then by increasing the noise amplitude a the stochastic population growth model ((5)–(9)) predicts the following scenario. For small values of the noise amplitude the mean population size grows in time and saturates finally at a value determined by the attracting interval (18). As a increases the mean population size slowly increases, but as the noise amplitude approaches the critical value ac, K01(ac)=K0c, the attracting interval (18) disappears and the system undergoes a transition to the state characterized by the attracting interval (17), where the population size is very much larger than in the interval (18). The further increase of a causes a monotonic and very slow decrease of 〈X〉 (cf. Figure 4).
The stationary mean value of the population size 〈X〉 versus the noise amplitude a. The system parameter values are K~0=100, ϵ=0.03, γ=10-5, α=0.5, and ν=75. Note the jump of the mean population size 〈X〉 at the critical value of the noise amplitude, ac≈0.66.
We emphasize that although an increase of the noise amplitude a causes an abrupt transition from a stationary state of a lesser number of individuals to one of a very much bigger number, by a decrease of a the opposite transition cannot occur. Let us look at the latter statement more closely on the assumption that at a>ac the initial population size x(0) is in the attracting interval (17). The population size evolves by the influence of an “effective flipping potential,” where the effective potential switches between two configurations Uefi (cf. Figure 5 and also Figure 2)
(30)Uefi(x)=-rγϵ∫xξFi(ξ1/n)1+γϵξ1+αdξ,i=1,2.
In the case of a>ac one potential configuration is monostable and the other is bistable (see Figure 5(a)) and the values of x are locked in the interval (x2,x1), that is, in the attracting interval (17). As the noise amplitude a decreases the monostable potential profile acquires a bistable profile at a≤ac, as a result of which the effective potential switches between two bistable configurations (see Figure 5(b)). But also in this case the population size x remains locked between the positions x12 and x22 of the minima of the configurations 1 and 2, respectively (i.e., in interval (19)). A similar phenomenon, called noise-induced single unidirectional transition, was previously considered in the context of symbiotic ecosystems [26], but the mechanism of that phenomenon is of a qualitatively different nature from the effect presented here.
A schematic illustration of the effective potential configurations Uefi versus the population size x, i=1,2. Because of the influence of the noise Z(t) the effective potential fluctuates with a switching rate ν between two configurations Uef1(x) and Uef2(x), which correspond to the noise values z=1 and z=-1 and are illustrated by curves 1 and 2, respectively. (a) The case of a>ac; the positions of the minima (stable fixed points of (1) with (3)) of Uef1(x) and Uef2(x) are denoted by x1 and x2, respectively. (b) The case of a<ac; the stable fixed points by the noise states z=1 and z=-1 are x1j and x2j, j=1,2, respectively. The initial population size is denoted by x(0) in (a).
6. Conclusions
Inspired by the Von Foerster growth model for human population [4], a generalization of the logistic growth model with a population size-dependent carrying capacity is given. Random influence of the environment on the dynamics of the population is taken into account by introducing a dichotomous noise in the carrying capacity. The exact expressions for the stationary probability density P~(x) of population sizes x are derived. One of our major results is the establishment of noise-induced transitions between three qualitatively different forms of a stationary population size distribution. At small noise correlation times the stationary distribution is characterized with a bell-shaped form, but as the correlation time of environmental fluctuations increases the system undergoes a transition to a state characterized with a probability density with either two or no local extrema, while a further increase of the correlation time brings another transition—the system goes to a state where the probability density exhibits a U-shaped form. Another result, perhaps the most important in the ecological context, is the existence of noise-induced single unidirectional transition for the mean population size (i.e., an increase in noise amplitude can cause an explosive growth of the size of population, while by decreasing the noise amplitude no opposite transitions can occur). It should be noted that noise-induced unidirectional transitions have been previously considered by investigations of symbiotic ecosystems [26], where an increase in noise amplitude can cause a catastrophic fall in the size of populations, but the mechanisms of those phenomena in [26] and in the current paper are qualitatively different.
Finally, we believe that the model and the results discussed here can be useful to describe the long-term behavior of the human population along the lines considered in [4–7, 30] and can also be applied to other types of processes that are explosive in their nature due to a positive feedback. A further detailed study is, however, necessary, especially an investigation of transient processes. Our exact results obtained in the case of dichotomous noise can be a good starting point for investigations, by numerical simulations, of the dynamics of the population size in the presence of more realistic models for environmental colored noise.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The work was supported by the Estonian Science Foundation under Grant no. 9005, by the Ministry of Education and Research of Estonia under Grant no. SF0130010s12, by the International Atomic Energy Agency under Grant no. 14797, and by European Union through the European Regional Development Fund (Centre of Excellence “Mesosystems: Theory and Applications,” TK114).
TurchinP.200335Princeton, NJ, USAPrinceton University PressMR2080584MurrayJ. D.198919Berlin, GermanySpringerMR1007836McGladeJ. M.1999Oxford, UKBlackwell ScienceVon FoersterH.MoraP. M.AmiotL. W.Doomsday: Friday, 13 November, A.D. 202619601323436129112952-s2.0-0003079808JohansenA.SornetteD.Finite-time singularity in the dynamics of the world population, economic and financial indices20012943-44655022-s2.0-003587314110.1016/S0378-4371(01)00105-4StrzałkaD.Connections between von Foerster coalition growth model and Tsallis q-exponential200940141472-s2.0-61449159772CohenJ. E.Population growth and earth's human carrying capacity199526952223413462-s2.0-0029148702BjørnstadO. N.GrenfellB. T.Noisy clockwork: time series analysis of population fluctuations in animals200129355306386432-s2.0-0035958646LaaksoJ.KaitalaV.RantaE.Population dynamic consequences of adaptive growth rate in fluctuating environment20061941–31321402-s2.0-3364466380810.1016/j.ecolmodel.2005.10.030BenziR.SuteraA.VulpianiA.The mechanism of stochastic resonance19811411L453L4572-s2.0-2064447059810.1088/0305-4470/14/11/006MR633455LaasK.MankinR.RekkerA.Constructive influence of noise flatness and friction on the resonant behavior of a harmonic oscillator with fluctuating frequency20097952-s2.0-6754910680710.1103/PhysRevE.79.051128051128GammaitoniL.HänggiP.JungP.MarchesoniF.Stochastic resonance19987012232872-s2.0-0039065101MankinR.LaasK.LaasT.ReiterE.Stochastic multiresonance and correlation-time-controlled stability for a harmonic oscillator with fluctuating frequency20087832-s2.0-5264914697810.1103/PhysRevE.78.031120031120MankinR.RekkerA.Memory-enhanced energetic stability for a fractional oscillator with fluctuating frequency20108142-s2.0-7795153774910.1103/PhysRevE.81.041122041122ReimannP.Brownian motors: noisy transport far from equilibrium20023612–4572652-s2.0-003617955710.1016/S0370-1573(01)00081-3MR1895137HänggiP.MarchesoniF.Artificial Brownian motors: controlling transport on the nanoscale20098113874422-s2.0-6434910502710.1103/RevModPhys.81.387García-OjalvoJ.SanchoJ. M.1999New York, NY, USASpringer10.1007/978-1-4612-1536-3MR1697197SaguésF.SanchoJ. M.García-OjalvoJ.Spatiotemporal order out of noise20077938298822-s2.0-3454722104510.1103/RevModPhys.79.829van den BroeckC.ParrondoJ. M. R.ToralR.KawaiR.Nonequilibrium phase transitions induced by multiplicative noise1997554408440942-s2.0-0000257575MangioniS.DezaR.WioH. S.ToralR.Disordering effects of color in nonequilibrium phase transitions induced by multiplicative noise19977913238923932-s2.0-0001523853HänggiP.MarchesoniF.GrigoliniP.Bistable flow driven by coloured gaussian noise: a critical study19845643333392-s2.0-000093938410.1007/BF01306642HorsthemkeW.LefeverR.198415Berlin, GermanySpringerMR724433SoikaE.MankinR.AinsaarA.Resonant behavior of a fractional oscillator with fluctuating frequency20108112-s2.0-7634909869110.1103/PhysRevE.81.011141011141MankinR.LaasK.SaugaA.Generalized Langevin equation with multiplicative noise: temporal behavior of the autocorrelation functions20118362-s2.0-7996104968210.1103/PhysRevE.83.061131061131MankinR.LaasK.LumiN.Memory effects for a trapped Brownian particle in viscoelastic shear flows201388404214210.1103/PhysRevE.88.042142SaugaA.MankinR.Addendum to ‘colored-noise-induced discontinuous transitions in symbiotic ecosystems’20057162-s2.0-2794444739810.1103/PhysRevE.71.062103062103KimS.ParkS. H.RyuC. S.Colored-noise-induced multistability in nonequilibrium phase transitions1998586799479972-s2.0-0001412496MüllerR.LippertK.KühnelA.BehnU.First-order nonequilibrium phase transition in a spatially extended system1997563265826622-s2.0-0031224179KlyatskinV. J.2011London, UKElsevierCabellaB. C. T.RibeiroF.MartinezA. S.Effective carrying capacity and analytical solution of a particular case of the Richards-like two-species population dynamics model20123914128112862-s2.0-8465517644210.1016/j.physa.2011.11.018