# Analysis of the Noise-Induced Regimes in Ricker Population Model with Allee Effect via Confidence Domains Technique

Authors: Irina Bashkirtseva; Lev Ryashko
Journal: BioMed Research International
 (2014)

---

## Abstract

We consider a discrete-time Ricker population model with the
Allee effect under the random disturbances. It is shown that noise can cause various
dynamic regimes, such as stable stochastic oscillations around the equilibrium, noise-induced extinction, and a stochastic trigger. For the parametric analysis of these
regimes, we develop a method based on the investigation of the dispersions and
arrangement of confidence domains. Using this method, we estimate threshold
values of the noise generating such regimes.


---

## Body

## 1. Introduction

Environmental noise is an inevitable attribute of any living system. Investigations of noise-induced phenomena in biological systems attract the attention of many researchers [1–4]. Obviously, noise changes the quantitative properties of the system dynamics. Moreover, even small deterministic disturbances and stochastic fluctuations can cause abrupt catastrophic qualitative shifts in ecosystems [5–7]. Theoretically, such shifts can be attributed to the multistability of the corresponding nonlinear mathematical models. Due to nonlinearity, these dynamic models exhibit coexisting attractors, nonuniformity of phase portraits, and high sensitivity of boundaries of basins of attraction. Under the random disturbances, a phase trajectory can cross a separatrix between basins of the attraction of coexisting attractors and exhibit new dynamical regimes [5, 8]. In spatial population models, noise can generate pattern formations [9–11].A classic example of noise-induced phenomena in ecosystems is an extinction of the population in the stochastic models with the Allee effect. The Allee effect means that there is a threshold population level below which the population goes to extinction. Deterministic population models with Allee effect are fairly well studied [12–18]. An analysis of stochastic population continuous-time models with the Allee effect is presented in [19–23].In the present paper, we study Allee effect in the discrete-time population Ricker model forced by additive and parametric noises. For discrete-time systems, an exhaustive mathematical description of the stochastic dynamics in terms of probabilistic distributions requires a solution of the Frobenius-Perron equation [24, 25]. An analytical solution of this functional equation even for one-dimensional systems is possible only in very special cases. In these circumstances, for the description of stochastic attractors, a method of direct numerical simulation is widely used. This method requires a lot of computing power, so a development of the analytical approximations is a highly relevant area of research [26–28].In the present paper, for the general one-dimensional discrete-time systems with parametric noise, we develop a new analytical method for the approximation of the dispersions of random states around stochastically forced equilibria. Mathematical background of this method is shortly presented in Section2.In Section3, for the modified Ricker population model with Allee effect, we study probabilistic mechanisms of the noise-induced extinction and generation of stochastic trigger regime. Here, we demonstrate constructive abilities of the new approach based on confidence domains technique.


## 2. Analysis of Stochastic Equilibrium

Consider a deterministic system(1)xt+1=f(xt),
where f(x) is a smooth scalar function. It is supposed that the system (1) has exponentially stable equilibrium x-. It means that f(x-)=x- and |f′(x-)|<1.Along with the deterministic system (1), consider a stochastically forced system
(2)xt+1=f(xt)+σ(xt)ξt,
where σ(x)=(σ1(x),…,σn(x)) is a smooth n-vector function and ξt=(ξ1,t,…,ξn,t)⊤ is uncorrelated n-vector random process with parameters Eξt=0, E(ξtξt⊤)=V, t=0,1,…. Here, n×n-matrix V defines second moments of the coordinates of the vector ξt. The function σ(x) describes a dependence of the intensity of random disturbances on the state of the system.For the deviationszt=xt-x- of the system (2) states xt from the equilibrium x-, the following first approximation system holds:
(3)zt+1=f′(x-t)zt+(σ(x-)+σ′(x-)zt)ξt,
where σ′(x)=(σ1′(x),…,σn′(x)). It follows from (3) that
(4)zt+12=[f′(x-)]2zt2+2f′(x-)[σ(x-)+σ′(x-)zt]ztξt+σ(x-)ξtξt⊤σ⊤(x-)+[σ(x-)ξtξt⊤σ′⊤(x-)+σ′(x-)ξtξt⊤σ⊤(x-)]zt+σ′(x-)ξtξt⊤σ′⊤(x-)zt2.Consider the dynamics of the first two momentsmt=Ezt, Dt=Ezt2 for the system (3). From (3) and (4), due to the noncorrelatedness of zt and ξt, it follows that
(5)mt+1=amt,Dt+1=bDt+μmt+α,
where
(6)a=f′(x-),b=a2+β,α=σ(x-)Vσ⊤(x-),β=σ′(x-)Vσ′⊤(x-),μ=σ(x-)Vσ′⊤(x-)+σ′(x-)Vσ⊤(x-).Suppose that the nonlinear randomly forced system (2) has a stochastic attractor with the stationary stable probabilistic distribution. For the approximation of the two first moments of the random states of this probabilistic distribution near equilibrium x-, we will use a stable stationary solution of the system (5).Due to the condition of stability|f′(x-)|<1, for any initial value m0, the sequence mt is stabilized:
(7)lim⁡t→∞mt=0.
A necessary and sufficient condition of the convergence of the sequence Dt to the constant value M is the inequality
(8)b<1.
Here,
(9)M=α1-b=σ(x-)Vσ⊤(x-)1-[f′(x-)]2-σ′(x-)Vσ′⊤(x-).Remark 1.
If all noises in system (2) are additive, the functions σi(x) are independent of x and, therefore, β=0, b=[f′(x-)]2, and condition (8) of the existence of stable stationary solutions of the system (5) is equivalent to the condition |f′(x-)|<1 of the stability of the equilibrium x-. If the intensity of noise depends on the system (2) state (β>0), then just the stability of the deterministic equilibrium is not sufficient. Here, condition (8) applies a restriction [f′(x-)]2+σ′(x-)Vσ′⊤(x-)<1 on the parameters σ′(x-) and V of noise. If this restriction is not satisfied, then the sequence of the second moments Dt given by the system (5) will increase indefinitely.Remark 2.
For the geometrical description of the scattering of random states, the confidence intervals are widely used. For scalar Gaussian random variable with mean valuex- and dispersion D, the confidence interval is (x--r,x-+r), where
(10)r=c2D,c=erf-1(P),erf⁡(x)=2π∫0x‍e-t2dt,
and P is a fiducial probability. Stationary distributed random states of the system (2) are localized in the neighborhood of the equilibrium x-. Using dispersion value M, one can construct a confidence interval (x--r,x-+r) where r=c2M.Remark 3.
Consider a stochastic system(11)xt+1=f(xt)+σ1(xt)ξ1,t+σ2ξ2,t
forced by only two noises: parametric noise with intensity σ1(x) and additive noise with intensity σ2. Here, ξ1,t and ξ2,t are uncorrelated scalar random processes with parameters
(12)E(ξi,t)=0,E(ξi,t2)=1(i=1,2);E(ξ1,tξ2,t)=0(t=1,2,…).
This system is a particular case of the system (2) with
(13)σ(x)=(σ1(x),σ2),ξt=[ξ1,tξ2,t],V=[1001],α=σ12(x-)+σ22,β=[σ1′(x-)]2,b=[f′(x-)]2+[σ1′(x-)]2.
It follows from (9) that, for this case,
(14)M=σ12(x-)+σ221-[f′(x-)]2-[σ1′(x-)]2.In the next section, we apply this theory to the study of noise-induced phenomena in Ricker model with Allee effect.


## 3. Analysis of Stochastic Ricker Model with Allee Effect

Consider a general one-dimensional population model governed by the following discrete-time equation:(15)Nt+1=g(Nt)Nt,
where N is a size of the population and g(N) is a per capita intrinsic growth rate function. Value N=0 is a trivial equilibrium of this system.The dynamics of the system (15) is defined by the features of the function g(N). The condition g>1 implies a growth of the population size; for g<1, the population decreases. Solutions of the equation g(N)=0 define the other equilibria of (15).The simplest examples of the functionsg(N) are plotted in Figure 1(a). As one can see, for sufficiently small values N0, due to g(N0)>1, the population grows. On the contrary, for large N0, the population decreases. The value N=K is a nontrivial equilibrium of the system (15). If this equilibrium is unstable, then the system (15) can exhibit different types of dynamics with periodic or chaotic oscillations around the equilibrium N=K.Intrinsic growth rate functiong(N): (a) without Allee effect; (b) with Allee effect.
(a)
(b)This type of dynamics is observed for well-known discrete-time models such as Verhulst equation withg(N)=1+r(1-(N/K)) (red line in Figure 1(a)) and Ricker system with g(N)=exp⁡[r(1-(N/K))] (blue line in Figure 1(a)). Such type of the function g(N) considered here adequately reflects a survival law: the more the population size, the less the growth rate.For small size of the population, the functiong(N) has to describe a birth rate mainly. But there are many ecological situations when at a low population level the function g(N) is significantly less than unit, and g(N) vanishes as the size of the population tends to zero. If we assume that the birth rate is proportional to the population size as the population level is low, one gets another type of the function g(N) (see Figure 1(b)).As one can see, this new function now provides another type of the population dynamics for small values ofN0. Indeed, one more equilibrium N=A for (15) appears. If N0<A, then the sequence Nt monotonically tends to zero. This means an extinction of the population. If N0 is slightly more than A, then the population grows and its future behavior is defined by features of the function g(N) for N>A. So, the value N=A is a threshold value separating two different types of population dynamics. Such a phenomenon of the existence of a threshold population level, below which the population goes to extinction, is called an Allee effect [12, 16].To incorporate an Allee effect into the classical Ricker model, we use the following modification [29] of the growth rate function:
(16)g(N)=NAexp⁡[r(1-NA)].
For this function, due to g(A)=1, (15) has an equilibrium N=A corresponding to the Allee threshold. Another nontrivial equilibrium N=K can be found from the following equation:
(17)Kexp⁡[r(1-KA)]=A.
For dimensionless new variable x=N/A, (15) can be rewritten as
(18)xt+1=xtφ(xt)=f(xt),φ(x)=xexp⁡[r(1-x)].
For the zone 0<r<1, the system (18) has three equilibria: x-0=0,   x-1=1, and x-2(r)>x-1. The equilibrium x-0 is stable, the equilibrium x-1 is unstable, and the equilibrium x-2(r) is stable for 0.1788<r<1. For r<0.1788, the system performs a Feigenbaum scenario of period-doubling bifurcations and transition to chaos. The unstable equilibrium x-1 separates basins of attraction of the stable equilibrium x-0 and attractors (both regular and chaotic) arranged in the zone x>x-1. Here, x-1 plays a role of the threshold population level of Allee effect: if the initial density x0 of the population is below x-1, then it goes to the extinction: lim⁡t→∞xt=0.In what follows, we focus on the case whenx-2 is stable and fix r=0.5. In Figure 2, for r=0.5, the dynamics of the system (18) for the different initial values x0 is shown. Here, stable equilibria x-0=0 and x-2=3.513 are plotted by green circles; the unstable equilibrium x-1=1 is plotted by red circle. A red square marks a right border of the interval of attraction of the equilibrium x-2.Figure 2
Dynamics of the deterministic Ricker model (18) with Allee effect for r=0.5.For the Ricker model in this deterministic case, there are only two variants of the dynamics: the population goes either to extinction or to the positive equilibrium. Such a strong separation of dynamical regimes can be destroyed by a noisy environment.Along with the deterministic model (18), we will consider the stochastic system
(19)xt+1=xt[φ(xt)+σ1ξ1,t]+σ2ξ2,t=xt2exp⁡[r(1-xt)]+σ1xtξ1,t+σ2ξ2,t,
where ξ1,t and ξ2,t are uncorrelated scalar random processes with parameters
(20)E(ξi,t)=0,E(ξi,t2)=1(i=1,2);E(ξ1,tξ2,t)=0.
Here, σ1 is an intensity of the parametric noise modeling the disturbances of the rate function φ(x), and σ2 is an intensity of the external additive noise.Under stochastic disturbances, the solutions of (19) leave the deterministic equilibria and form some stationary probabilistic distributions around points x-0 and x-2. For the analysis of the dispersions of the random states, the approximations (14) are used. For dispersions in the system (19), around the equilibrium x-0 we have
(21)M0=σ221-σ12,
and around the equilibrium x-2 we have
(22)M2=σ12x-22+σ221-(f′(x-2))2-σ12.When the noise intensity is quite small, random states leaving a stable deterministic equilibrium are concentrated around it with the small dispersion. As the noise intensity increases, a dispersion of random states increases too, and the system (19) can exhibit qualitative changes of stochastic dynamics.Upon reaching a certain critical value of the noise intensity, iterations of the stochastic system (19) with a high probability pass through the unstable equilibrium x-1 into the basin of attraction of the stable equilibrium x-0 and perform small-amplitude stochastic oscillations near x-0. Biologically, this phenomenon can be interpreted as a noise-induced extinction of the population. Evaluation of the critical noise intensity corresponding to the beginning of these transitions may be obtained on the base of confidence domains technique. Functions (21) and (22) give us explicit parametrical description for confidence intervals I0=(x-0-r0,x-0+r0) and I2=(x-2-r2,x-2+r2), where ri=erf-1(P)2Mi.In Figures3–5, the dependence of the stochastic attractors on the intensity σ1 of parametric noise for the fixed additive noise intensity σ2 is presented. Here, random states (grey color) of attractors have been found by the direct numerical simulation. The stable equilibria x-0, x-2 of the deterministic system are plotted by green lines, the unstable equilibrium x-1 is plotted by red line, and the nearest borders x=x-0+r0(σ1) and x=x-2-r2(σ1) of confidence intervals I0 and I2 are plotted by dashed blue lines. Here, the representative time series are also shown.Noise-induced transitions forσ2=0.1.
(a)
(b)Noise-induced transitions forσ2=0.
(a)
(b)Noise-induced transitions forσ2=0.3.
(a)
(b)As one can see in Figure3(a) for the fixed σ2=0.1, when σ1 exceeds threshold value σ1*≈0.2, random trajectories with high probability demonstrate noise-induced transitions from the basin of attraction of the nontrivial equilibrium x-2 to the small vicinity of the equilibrium x-0.An estimation of the critical value ofσ1* can be derived from the mutual arrangement of the separatrix x-1 and confidence intervals borders. In our case, the noise intensity that corresponds to the intersection of the lower border x=x-2-r2(σ1) of the confidence interval I2 with the separatrix x=x-1 can be used as an estimation of the threshold value σ1*≈0.17. Here, the fiducial probability is P=0.999.As one can see, the valueσ1* quite accurately localizes a qualitative change in the stochastic dynamics of the studied system. Note that, after transition to the neighborhood of the trivial equilibrium, a dispersion of random states is well described by the confidence interval found for x-0=0. In Figure 3(b), a solution of the system (19) for σ1=0.2, σ2=0.1 illustrates the transition from the neighborhood of x-2 to the neighborhood of x-0.Similar results can be obtained for the case when the system (19) is forced by parametric perturbations only (σ1≠0,σ2=0). In this case (see Figure 4), solutions of (19) with high probability cross the separatrix x-1=1 and quickly vanish. Biological interpretation of this transition is a noise-induced extinction of the population.As noise grows, both confidence intervals around thex-0 and x-2 expand beyond the boundary x=x-1 of the basins of attraction of these stable equilibria and begin to overlap each other (see Figure 5(a) for σ2=0.3). As a consequence, the system exhibits repeated random transitions between these basins. Figure 5(b) for σ1=0.2, σ2=0.3 demonstrates two types of stochastic oscillations: oscillations near equilibria and noise-induced transitions between the neighborhoods of equilibria. In this case, the population can be characterized as a stochastic trigger.So, in the presence of noise, the Ricker population model with Allee effect represents various types of stochastic dynamics: long-term random oscillations with the small dispersion near the positive equilibrium, noise-induced extinction, and stochastically generated trigger. The new mathematical technique presented here provides a useful tool for the analytical estimation of the noise intensities corresponding to these regimes.


## 4. Conclusion

In this paper, we propose a new approach for the constructive study of noise-induced phenomena in the Ricker population model with the Allee effect. Theoretical basis of this approach is an analysis of dispersions of random states near deterministic equilibria of this bistable model. Using this theory, we construct confidence intervals around the stable equilibria. The mutual arrangement of these intervals and Allee threshold are used in the geometrical analysis of the various stochastic regimes in this model: stable stochastic oscillations around nontrivial equilibrium, noise-induced extinction, and stochastic trigger. Note that the elaborated method is readily applicable to more complicated models of multiple interacting populations forced by parametric noises.



---

*Source: 346239-2014-05-28.xml*