# The Bifurcation and Control of a Single-Species Fish Population Logistic Model with the Invasion of Alien Species

Authors: Yi Zhang; Qiaoling Zhang; Jinghao Li; Qingling Zhang
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

The objective of this paper is to study systematically the bifurcation and control of a single-species fish population logistic model with the invasion of alien species based on the theory of singular system and bifurcation. It regardsSpartina anglica as an invasive species, which invades the fisheries and aquaculture. Firstly, the stabilities of equilibria in this model are discussed. Moreover, the sufficient conditions for existence of the trans-critical bifurcation and the singularity induced bifurcation are obtained. Secondly, the state feedback controller is designed to eliminate the unexpected singularity induced bifurcation by combining harvested effort with the purification capacity. It obviously inhibits the switch of population and makes the system stable. Finally, the numerical simulation is proposed to show the practical significance of the bifurcation and control from the biological point of view.


---

## Body

## 1. Introduction

Singular system, also known as differential-algebraic system, is relatively accurate model description of some problems in the practical applications. The theory of singular system has been widely applied to various fields [1–4]. Both the Leontief dynamic input-output model and Hopfield neural network model are differential algebraic-systems [5, 6] and so are dynamic model of multiple robot body coordination operation and power system model with nonlinear load. Moreover, the study of the various bifurcations in the power system has produced great significance to reveal the mechanism of voltage instability [7–9]. It is known through [10–12] that the bifurcation theory of singular system already has had a solid theoretical foundation. Recently, Zhang and other scholars studied several kinds of biological dynamical systems [13–18] by utilizing the theory of singular system, which generates a profound impact on the better description of development situation for biological resources and the more accurate grasp of the real law for exploitation of biological resources.This paper is a study on the aquatic plantsSpartina anglica. According to the situation of our country, the introduction of aquatic species abroad accelerates so that the concern over the security of aquatic invasive species has also grown with the rapid development of aquaculture industry and the need of propagation and release for resource.Spartina anglica, which is one of the aquatic plants introduced abroad by the water conservancy department, is used to promote the sedimentation and reclamation. Now, it has become a big disaster. A few years ago, the shoals of Ling-kun Island were overrun with theSpartina anglica which occupied the tens of thousands of acres of tidal flats so that the fish and shellfish of the tidal flats vanished. It took a lot of manpower and material resources to burn and cut it [19]. In consequence, the further study on the influence of alien species on the singular ecological economic system has become an irresistible trend. In order to protect the ecological environment and guarantee the economic development simultaneously, many biological economic models are established by some mathematical researchers, and some valuable results are obtained. For example, the problems of the optimal harvesting in the biological economic model are studied, respectively, in [20–23], which provide the theory basis for reasonable development of biological resources. The existence and stability of equilibrium solution and positive periodic solution in the biological economy model are discussed, respectively, in [24, 25]. However, the singular biological economic system which further considers the purification capacity and harvested effort in the case of invasive species, according to the author’s knowledge, has never been studied systematically.The innovation of this paper considers purification capacity for theSpartina anglica, which is first proposed by us in this field. SinceSpartina anglica can promote the sedimentation and reclamation, increasing interest has been stimulated to study it. On one hand decreasing the density ofSpartina anglica can facilitate survival of fish and other aquatic organisms; on the other hand it cannot be too low to play a favorable part in promoting the sedimentation and reclamation. The character of soil can be improved so that the environment can be protected. In consequence, it is essential to control the density ofSpartina anglica. Besides, the purification capacity plays an important role in the economically sustainable harvesting capacity of the fish population, which at least provides space and adequate nutrients for the growth of fish. In fact, due to the appearance of the purification capacity, some complex dynamic characteristics will arise, such as bifurcation, impulse, and other characteristics, which will affect the stability of the system. Specially, the model in this paper exhibits two bifurcation phenomena when the economic revenue is zero. One is trans-critical bifurcation which changes the stability of the system, and the other is singular induced bifurcation which indicates that zero economic revenue brings impulse, and from the biological aspect, the impulse means a rapid expansion of the population. There are both theoretical significance and practical value to deeply study the problems of bifurcation and impulse, which are restrained with the purification capacity for invasion of alien species and can further reveal the instability mechanism of the whole biological economic system. Therefore, the proposed model is significantly different from the existing models in the field of research which use differential algebraic equations approach and needs separate investigation.


## 2. Modeling

Introducing a density restriction factor(1-N/K) to Malthus model leads to the famous Logistic equation on ecology:
(1)dNtdt=rNt1-NtK,
where N(t) is the density or quantity of population at time t, r is the intrinsic growth rate of the population, and K is the environmental capacity.To be more reasonable, we introduce alien species to this model. Assuming thaty(t) is the density ofSpartina anglica at time t, the original water nutrients, water, oxygen, and living space are used due to the invasion of alien species, and we get the following model:
(2)dxtdt=rxt1-xt+ytK-Etxt,dytdt=βyt-hE1t,dE1tdt=θyt-wE1t-αxt,
where x(t), E1(t), E(t) are, respectively, the density of fish population, the purification capability forSpartina anglica, and the harvested effort to fish at time t. β is the intrinsic growth rate of theSpartina anglica. -hE1(t) is the quantity of purification toSpartina anglica. θy(t) is the growing degree ofSpartina anglica. w is the cost of a unit of purification effort. -αx(t) represents that the purification capability is stimulated by the reduction of fish to a certain extent. The constants mentioned above are supposed to be positive.Gordon founded the open or public fishery economic theory in 1954 [26]. Sustainable Economic Profit = Sustainable Total Revenue − Sustainable Total Cost, when the harvested effort E(t) is given. Therefore, when the harvested effort E(t) switches with time t, and we get the following algebraic equation:
(3)Etxtp-c=mt,
where p is the market price of the captive population. c is the cost of a unit of capture effort. m(t) is the net economic revenue. Based on (2) and (3), the following singular biological economic system is established:
(4)x˙t=rxt1-xt+ytK-Etxt,y˙t=βyt-hE1t,E1˙t=θyt-wE1t-αxt,0=Etxtp-c-mt.Equation (4) also can be expressed as the following matrix form:
(5)AtX˙t=Gxt,yt,E1t,Et,
where
(6)Xt=xt,yt,E1t,Et,At=1000010000100000,Gxt,yt,ut,Et=rxt1-xt+ytK-Etxtβyt-hE1tθyt-wE1t-αxtEtxtp-c-mt.The objective of this paper is to study the dynamic behavior of the system (4) and propose the singular state feedback controller to stabilize the bifurcation. The structures of the reminder of this paper are arranged as follows. In Section 3, the stability of equilibria is discussed. In Section 4, the dynamic behaviors around equilibria are studied. In Section 5, state feedback controller is designed. Finally, a real-world example is given to discuss the practical significance of the bifurcation and control from biological perspective.


## 3. Equilibria and Stability

Whenm=0, (4) can be written as
(7)x˙t=rxt1-xt+ytK-Etxt,y˙t=βyt-hE1t,E1˙t=θyt-wE1t-αxt,0=Etxtp-c.
It is clear that (7) has three equilibria:
(8)P10,0,0,0,P2Khθ-βωM,hKαM,KαβM,0,P3cp,chαphθ-βω,cαβphθ-βω,r-crMKphθ-βω,
where M=hα+hθ-βω. Next, we consider the stability in the neighborhood of each equilibrium.Theorem 1.
(1) Equation (7) is unstable at the equilibrium P1(0,0,0,0).
(2) Ifω-β>0 and (c-(Kp(hθ-βω)/(hα+hθ-βω)))≠0, then (7) is stable at the equilibrium
(9)P2Khθ-βωM,hKαM,KαβM,0.Proof.
The Jacobian matrix of (7) at P1 is
(10)JP1=-rxK+r1-x+yK-E-rxK0-x0β-h0-αθ-ω0Ep00xp-cP1=r0000β-h0-αθ-ω0000-c.
Then the characteristic polynomial is
(11)det⁡⁡λA-JP1=c-r+λ×hθ-βλ+λ2-βω+λω=0.
The characteristic roots are
(12)λ1=r>0,λ2=12β-ω-β2-4hθ+2βω+ω2,λ3=12β-ω+β2-4hθ+2βω+ω2.
It can be judged that P1 is an unstable equilibrium.
Similarly, forP2(K(hθ-βω)/M,hKα/M,Kαβ/M,0), in order to guarantee the fish population,Spartina anglica and purification capability all exist, and the following inequalities need to be satisfied:
(13)Khθ-βωM>0,hKαM>0,KαβM>0,thatishθ-βω>0,
where
(14)M=hα+hθ-βω=hα+hθ-βω>0.
In order to study the stability of P2, we firstly obtain the Jacobian matrix of (7) at P2:(15)JP2=-hrθ+rβωhα+θ-βω-hrθ+rβωhα+θ-βω0-hKθ+Kβωhα+θ-βω0β-h0-αθ-ω0000-c+Kphθ-βωhα+hθ-βω. Based on det⁡(λA-Jp2)=0, then
(16)c-Kphθ-βωM×-h-θλ-hrαθM-hrθ2M+rαβωM+rβθωM+-β+λλ+ωλ--hrθ+rβωM-θλ-hrαθM-hrθ2M+rαβωM+rβθωM=0,
where M=hα+hθ-βω.
It is a cubic equation ofλ. When (c-(Kp(hθ-βω)/(hα+hθ-βω)))≠0, it can be simplified as follows:
(17)λ3+Δ1λ2+Δ2λ+Δ3=0,
where
(18)Δ1=-β+ω+r(hθ-βω)h(α+θ)-βω,Δ2=(hθ-βω)(h(α+θ)-βω+r(-β+ω))h(α+θ)-βω,Δ3=rhθ-βω.
Since ω-β>0, hθ-βω>0, Reλ1<0, Reλ2<0, Reλ3<0.
WhenΔ1>0, Δ2>0, Δ3>0, Δ1×Δ2-Δ3>0. Then based on the Routh-Hurwitz [27] theorem, (7) is stable at P2.It is explained through this theorem that the system keeps stable under certain condition that fish population,Spartina anglica, and purification capacity all exist while harvested effort does not exist. Since the economic interests are obtained through the harvest effort, we always pay more attention to the dynamic character of the positive equilibrium. Straightforward computation shows that P3 is the singular point [28] of (7), and the dynamics is relatively complex. The discussion will be shown in Section 5.


## 4. Trans-Critical Bifurcation

The topological structure of (4) will be switched when the biological phenomenon of economic balance appears. We consider the situation in the neighborhood of P1(0,0,0,0).Theorem 2.
Ifm=0, (hθ-βω)≠0, α≠(-hθ+βω), then (4) undergoes the trans-critical bifurcation at P1(0,0,0,0), and r=0 is the bifurcation value.Proof.
For (4), let
(19)fx,y,E1,E,m=rxt1-xt+ytK-Etxtβyt-hE1tθyt-ωE1t-αxt,gx,y,E1,E,m=Etxtp-c-mt.
We chooser as the bifurcation parameter. If m=0, then det⁡(DEg)p1=-c≠0 (DE=∂/∂E, similarly hereinafter). Applying the implicit function theorem, in certain neighborhood of P1(0,0,0,0), (4) can be simplified to the ordinary differential equation X˙=fR(X,r), where X=(x,y,E1)T, fR(X,r)=(fR1(X,r),fR2(X,r))T, which are adequately smooth two-dimensional vector functions. fR(0,r)=0→ and the matrix
(20)DXfR0,0=DXf-DEfDEg-1DXg=-E000β-h-αθ-ω-1xp-c-x00Ep00P1=0000β-h-αθ-ω.
It is clear that the characteristic roots of the matrix above are
(21)λ1=0,λ2=12β-ω-β2-4hθ+2βω+ω2,λ3=12β-ω+β2-4hθ+2βω+ω2.
There are one simple zero characteristic root and two nonzero characteristic roots. Through computation, the right eigenvector of the zero characteristic root is
(22)Ψ=--hθ+βωαβhβ1,
and the left eigenvector is Φ=100. Consider
(23)Φ∂2fR∂r∂XΨ=100100000000--hθ+βωαβhβ1=--hθ+βωαβ≠0,Φ∂2fR∂X2Ψ,Ψ=100ΨTΛ1ΨΨTΛ2ΨΨTΛ3Ψ=hr-hθ+βωKαβ2-rhθ-βωKαβ×2hθ-βωαβ+hβ=2hr-hθ+βωKαβ21--hθ+βωα≠0,
where
(24)Λ1=-2rK-rK0-rK00000,Λ2=000000000,Λ3=000000000.
Based on the definition of trans-critical bifurcation in [29], (4) undergoes the trans-critical bifurcation at P1(0,0,0,0).Obviously, the trans-critical bifurcation changes the stability of the system. In fact, the phenomenon of trans-critical bifurcation means that, for some reason, such as, the shortage of nutrition, water, and oxygen with the invasion ofSpartina anglica, fish population loses the natural fertility. The ecological interpretation of trans-critical bifurcation is that the phenomenon of economic overfishing [26] appears so that the extinction of population will arise in the biological economy system at the above depicted moment. The phenomenon admonishes people simultaneously that althoughSpartina anglica can promote the sedimentation and reclamation and improve the soil, its overproduction can also lead to a breakdown in local ecosystem. That means the density of the fish population,Spartina anglica, harvested effort, and purification capacity are all zero. In consequence, it is necessary to guarantee the fish species reproduce normally to maintain the sustainable survival of population when introducing the alien species to promote silting reclamation and protect the environment.


## 5. The Singularity Induced Bifurcation

We pay more attention to local dynamic characteristics near the positive equilibrium of the system in the actual situation of biological economics. Therefore, we obtain the following theorem.Theorem 3.
If(hθ-βω)>0, (r-(cr(hα+hθ-βω)/Kp(hθ-βω)))>0, then (4) will show the phenomenon of singularity induced bifurcation at
(25)P3cp,chαphθ-βω,cαβphθ-βω,r-crMKphθ-βω,
and m is a bifurcation parameter so that the system loses stability, where M=hα+hθ-βω.Proof.
In order to guarantee that each component exists atP3, the following inequalities need to be satisfied:
(26)cp>0,chαphθ-βω>0,cαβphθ-βω>0,r-crMKphθ-βω>0,
which means (hθ-βω)>0, r-(crM/Kp(hθ-βω))>0. Consider m as the bifurcation parameter, X=(x,y,E1)T, P3=(x*,y*,E1*,E*), and the singular surface
(27)S:=X,E,m∈Rn+m+p:gX,E,m=0,ΔX,E,m=0,
where Δ(X,E,m):=det⁡[DEg(X,E,m)].
Whenm=0, P3 passes through the singular surface S so that
(28)DEgP3=xp-cP3=0.
Therefore DEg has a zero characteristic root. 							(29)traceDEfadjDEgDXgP3=traceDEfadjDEgDxgDygDE1gP3,trace-xEp00000000P3=-cr-c2rMKphθ-βω=-cE*≠0.
Then dimker(DEgP3)=1 means that the dimension of  kernel is one, and
(30)det⁡DEgxm,ym,E1m,Em,m≠0whenm≠0,
where (x(m),y(m),E1(m),E(m)) represents the equilibrium curve of the initial value at P3(x(0),y(0),E1(0),E(0)).
Through computation, we get(31)DXfDEfDXgDEgP3=-r-crKp+crhα+hθ-βωKphθ-βω+r1-c/p+chα/phθ-βωK-crKp0-cp0β-h0-αθ-ω0pr-crhα+hθ-βωKphθ-βω000=-p-chθp+cβωpr-crhα+hθ-βωKphθ-βω=c(hθ-βω)E*≠0DXfDEfDmfDXgDEgDmgDXΔDEΔDmΔP3=-r-crKp+crhα+hθ-βωKphθ-βω+r1-c/p+chα/phθ-βωK-crKp0-cp00β-h00-αθ-ω00pr-crhα+hθ-βωKphθ-βω000-1p0000=chθ-cβω≠0,
where Δ(X,E,m)=xp-c, and
(32)B=-traceDEfadjDEgDXgP3=-trace-xEp00000000P3=cr-c2rhα+hθ-βωKphθ-βω,C=DuΔ-DXΔDEΔDXfDEfDXgDEg-1DmfDmg=-p000DXfDEfDXgDEgP3-1000-1=1r-crhα+hθ-βω/Kphθ-βω.
Clearly, B/C=c(r-(cr(hα+hθ-βω)/Kp(hθ-βω)))2=cE*2>0. Based on theorem 3 in [28], we know that one characteristic root of (4) moves from C- to C+ and along the real axis diverges to infinity when m goes from negative to positive. And
(33)dimkerDEgP3=1,traceDEfadjDEgDXgP3≠0,DXfDEfDXgDEgP3≠0,det⁡DEgxm,ym,E1m,Em,m≠0,
where (x(m),y(m),E1(m),E(m)) represents the equilibrium curve of the initial value P3(x(0),y(0),E1(0),E(0)).Based on the theorem  1.1 in [30], the other characteristic root of (4) in the neighborhood of m=0 is continuous, which means that it will not affect the switch of stability near the equilibrium. And through computation we know that this characteristic root has negative real part. In consequence, the system switches from being stable to being unstable as m moves from being negative to being positive through zero.Singularity induced bifurcation is a kind of bifurcation which arises in singular system, and it never occurs in an ordinary differential equation. The singularity induced bifurcation that the system (4) shows indicates that the real parts of some eigenvalues for the linear matrix (DXf-DEf(DEg)-1DXg) of the system (4) will change as bifurcation parameter and lead to the changes of system stability when DEg is singular. In addition, singularity induced bifurcation also has a great impact on the topological structure near the equilibrium point of the system, which leads to the emergence of pulse and even the collapse of the whole system [31].From the ecological prospective, singularity induced bifurcation and the system instability will cause the following effects on the ecological system. The collapse which singularity induced bifurcation leads to will be embodied by a surge increase of a population density within a short period in the ecological system and loss of the stability of the system. On one hand, the fierce competitions for the limited resources and survival occur within a population, which results that the young population in the inferior position would die in competition due to shortage of resources. There is no benefit to the sustainable development of this population. On the other hand, surge of a population may affect various food chains and lead to competition of limited living resource among populations. Some other populations in inferior position will face the extinction due to the shortage of resources, which is not beneficial to the maintenance of the population diversity in the ecological system.Accordingly, in this paper, the singularity induced bifurcation may cause an impulsive phenomenon of the single-species ecosystem especially with the invasion of alien species which can occupy the environmental and biological resources. On the other hand, the purification capacity plays an important role in the economically sustainable harvesting capacity of the fish population, which at least provides space and adequate nutrients for the growth of fish. The singularity induced bifurcation means that if fishers purify theSpartina anglica, the phenomenon of economic overfishing [26] appears when the total income is equal to the total cost. The rapid increase of the density ofSpartina anglica and the rapid decline of fish population can lead to the fact that one characteristic root of this system will tend to infinite and impulse phenomenon arises. It means the quantities of fish population and alien species switch in an instant even beyond the capacity of the environment. Surge increase of the population consumes a large number of nutrients in water so that nutrient concentration decreases. The decline of nutrient concentrations, which in turn affects breeding of populations, leads to mortality of a large number of populations and causes ecological imbalance. The existence of a singularity induced bifurcation phenomenon is verified at the equilibrium of the model system with zero economic revenue. It should be pointed out that the trapping behavior and purifying behavior cannot be unlimitedly accessed to the populations in the ecological system and they need to be controlled in a certain range. Otherwise the population resources would lose due to the excessive capture and purification, which affects the relevant populations in the food chain and the sedimentation and reclamation bySpartina anglica. It is not conducive to the sustainable development of the population diversity in the ecological system.To solve the above problem, we will introduce a state feedback controller to cure the ecological unbalance by some artificial methods in the next part.


## 6. Feedback Control of the Singularity Induced Bifurcation

From a long-term perspective, we should take some appropriate strategies to control harvested effort and increase purification capability so that gains keep in a proper range to eliminate singularity induced bifurcation.Let(34)ut=k1Et-E*+k2E1t-E1*,
where k1, k2 are feedback gains; then the closed-loop system is
(35)x˙t=rx1-x+yK-Ex,y˙t=βy-hE1,E1˙t=θy-wE1-αx+k2E1-E1*,0=Exp-c+k1E-E*.
It is easy to verify that (35) is locally controllable.Theorem 4.
Let(36)ut=k1Et-E*+k2E1t-E1*,
where k1, k2 are feedback gains and they satisfy(37)k1>cr2crβ+hKpθKphθ-βω-chα+θ-βω2c2r2β+K2p2β-β2+hθ+cKpr-2β2+hα+θhθ-βωk2<min⁡c-hrα+2rβω+Kpβ3-2hβθ+hθω2crβ+hKpθ,Kpβ4-hβ2θ-h2θ2-β3ω+2hβθω+cr2β3+h-αβ-2βθ+αωchrα-Kpβ(β2-2hθ),2c2hr2α-K2p2β2-hθ-hθ+βω+cKpr2hβθ+hαω-2β2ωKpcrhα-2β2+Kpβ-β2+hθ; then (35) is stable at P3(x*,y*,E1*,E*), where P3(x*,y*,E1*,E*) = (c/p,chα/p(hθ-βω),cαβ/p(hθ-βω),r-(cr(hα+hθ-βω)/Kp(hθ-βω))).Proof.
The Jacobian matrix of (35) at P3 is (38)Jp3=-r-crKp+crhα+hθ-βωKphθ-βω+r1-c/p+chα/phθ-βωK-crKp0-cp0β-h0-αθk2-ω0pr-crhα+hθ-βωKphθ-βω00k1. Then the characteristic polynomial is λ3+Δ1′λ2+Δ2′λ+Δ3′=0, where
(39)Δ1′=-ck1E*+crKp-β+ω-k2,Δ2′=cβk1E*-cωk1E*+hθ-βω+ck1E*k2-crβ-ω+k2Kp+βk2,Δ3′=-chθk1E*+cβωk1E*-cβk1E*k2+crhα+θ-βω+βk2Kp,E*=r-crhα+hθ-βωKphθ-βω.
Its characteristic roots are λ1, λ2, λ3. When Δ1′>0, Δ2′>0, Δ3′>0, Δ1′×Δ2′-Δ3′>0. Then based on the Routh-Hurwitz [27] theorem, (35) is stable at P3(x*,y*,E1*,E*). Thus, k1,k2 satisfy
(40)-ck1E*+rKx*-β+ω-k2>0,cβk1E*-cωk1E*+hθ-βω+ck1E*k2-rβ-ω+k2Kx*+βk2>0,-chθk1E*+cβωk1E*-cβk1E*k2+rhα+θ-βω+βk2Kx*>0,
where x*=c/p. Simplify the inequalities above, and we have (41)k1>cr2crβ+hKpθKphθ-βω-chα+θ-βω2c2r2β+K2p2β-β2+hθ+cKpr-2β2+hα+θhθ-βω,k2<min⁡⁡c-hrα+2rβω+Kpβ3-2hβθ+hθω2crβ+hKpθ,Kpβ4-hβ2θ-h2θ2-β3ω+2hβθω+cr2β3+h-αβ-2βθ+αωchrα-Kpββ2-2hθ,2c2hr2α-K2p2β2-hθ-hθ+βω+cKpr2hβθ+hαω-2β2ωKpcrhα-2β2+Kpβ-β2+hθ. At this point, we obtain Reλ1<0, Reλ2<0, Reλ3<0. Thus (35) is stable at P3.Notice. In practice, it is not enough to only control harvested effort. It is necessary to control the population density of alien species and release certain density of fish population. This paper studies the problem of alien species invasion. Controlling harvested effort and purification capability and taking k1, k2 appropriate values can eliminate the singularity induced bifurcation of the biological economic system. In other words, the biological economic system can be stabilized by increasing purification capability, controlling harvesting effort, and other strategies in the case of zero income temporarily.


## 7. Application Instance

Spartina anglica has played a positive role in protecting dam to shape earth and lessening wave at the beginning of the planting at the coast of eastern Fujian Province. However, some other characteristics ofSpartina anglica are highlighted at later period, such as strong fertility, developed root system, seeds fluttering around in the wind to spread rapidly, and the lack of natural enemies. At the present time, the tidal flats of 20 coastal towns are encroached in 4 counties (districts or cities) of Jiaocheng, Fuan, and Xiapu. The polluted area has been 110 thousand acres, which accounts for more than half of tidal flats of fish-farming.Spartina anglica, which is described as the cannibal grass, has already been overrun at the coast of Shandong province in China, especially in the Yellow River delta of the national conservation area at the Yellow River estuary of Dongying City.Spartina anglica today mainly aggregates in shellfish beds in south of estuary of fairy ditch, Xianhe Town, Hekou District, and Dongying City. It tends to overgrow at the speed of 6 times per year on the beach of 10 km2 approximately, which is a serious threat to marine life.Spartina anglica invades the coastal tidal flats of Dongying City, and the cover is more and more extensive so that the damaged area has been 86616 km2. The visible area ofSpartina anglica is about 3333 to 4000 km2, and the seed drifting area has been over 6666 km2. In addition, the scattered and small pieces ofSpartina anglica are also found in Binzhou City, Weifang City, and Yantai City [32].Based on the relevant data in Chinese alien marine species basic information database [32], we properly process the data (nondimensional transformation, equal ratio simplification, approximation, and so on) according to the parameters in the model of this paper. In sequence, the values of parameters are as follows:
(42)r=1,K=4,β=1.5,h=3,θ=2,ω=2,α=1,p=0.8,c=0.4,m=0.
Equation (4) can be written as
(43)x˙t=xt1-xt+yt4-Etxt,y˙t=1.5yt-3E1t,E1˙t=2yt-2E1t-xt,0=Et0.8xt-0.4.
And by evaluating, the positive equilibrium P3 is (0.5,0.5,0.25,0.75). Based on Theorem 3, (43) undergoes singularity induced bifurcation near P3. According to the discussion in Section 6, let
(44)ut=k1Et-0.75+k2E1t-0.25,
where k1, k2 are feedback gains to be determined.Then the closed-loop system is(45)x˙t=xt1-xt+yt4-Etxt,y˙t=1.5yt-3E1t,E1˙t=2yt-2E1t-xt+k2E1t-0.25,0=Et0.8xt-0.4+k1Et-0.75.
Based on Theorem 4, when k1, k2 satisfy
(46)k1>0.31k2<min⁡-0.35,-1.09,-1.71,(45) is stable at P3, and the singularity induced bifurcation of (43) is eliminated. The dynamic responses of (45) before and after control are illustrated in Figures 1 and 2.Figure 1
The dynamic responses of (45) near equilibrium before control.Figure 2
The dynamic responses of (45) near equilibrium (k1=4, k2=-1.75).Figure1 shows that the purification capacity fluctuates as the density of theSpartina anglica fluctuates. At the same time the density of the fish also fluctuates which is less choppy than the above two terms. The abundance ofSpartina anglica, the small number of fish, and the sustainable purification capacity together make the harvesting capacity nonexistent before controlling of the purification capacity and the harvesting capacity. In comparison with Figure 1, Figure 2 reflects the effectiveness of the feedback controller. It is noted that the four variations are sustainable at expected levels after ephemeral fluctuation in Figure 2. Furthermore, the density of fish andSpartina anglica are at the same level 0.5 and the harvesting capacity is much higher than the purification capacity, which ensures the economic revenue to be existent. It effectively shows the purpose of feedback control and demonstrates convincingly the theory proposed in this paper.


## 8. Conclusions

There exists broad recognition that the fisheries and aquaculture are under stress.The invasion ofSpartina anglica has been one of the most intimidating and dangerous factors to threaten the fisheries and aquaculture. Clearly, new management approaches or options must be considered to stem the damage and ensure that the fisheries and aquaculture ecosystems, as well as their unique features, are protected and restored. Current management of fish resources is insufficient. In this regard, purification capacity is increasingly proposed as major tools to relieve the stress on environmental resources and ecosystems. On the basis of this requirement and current trends, using differential algebraic systems theory, we analyze the dynamical behavior of a single-species fish population logistic model by treating the invasion ofSpartina anglica as a regulatory mechanism.This paper has mainly investigated the problems of bifurcation and control for a singular species fish population with logistic growth by incorporating invasion of alien species and further gotten the following conclusions. At the zero equilibrium, the system will raise trans-critical bifurcation so that the fish population is extinct. At this moment, the local ecological cycle of the system will breakdown and the respective component will not exist. It may be concluded that the ecosystem will collapse owing to the massive harvested effort of fish and disregarding of the control ofSpartina anglica for a long time. It is observed that the alien speciesSpartina anglicais the activator for the extinction. At the positive equilibrium, the system will raise singularity induced bifurcation. The quantity of fish population and alien speciesSpartina anglicaswitch at the instant, even beyond the capacity of the environment, eventually lead to ecological disaster and ecological imbalance. It is evident from the obtained results that the alien species can cause a stable equilibrium to become unstable, and even an impulsive phenomenon may occur when the economic revenue goes through its critical value. It is also clear that when the purification capacity is not adequate, both the fish and alien species populations reach immeasurable change around the equilibrium. As the density of the alien species increases, the fish population dies out. Since the environmental resources are finite, the alien species dies out eventually. It is obvious that both the purification capacity and the harvested effort are necessary to be considered.In consequence, theSpartina anglica, which is introduced for the purpose of improving soil and protecting environment, will lead to ecological destruction with the improper control of the density. Since this paper studies the alien speciesSpartina anglica, it is inevitable to introduce the purification capacity and controlling; moreover both too high and too low density ofSpartina anglica cannot attain the aim of the introduction. On the other hand, the purification capacity plays an important role in the economically sustainable harvested effort of the fish population which at least provides space and adequate nutrients for the growth of fish.It should be noted that the proposal of purification makes the work studied in this paper have some new and positive features. The feedback controller designed in this paper eliminates the singularity induced bifurcation and inhibits the switch of population to make the system stable by controlling purification capability and harvested effort. It is deemed important to undertake this type of study for the purpose of investigating the impact of purification so that sustainability of the ecosystem may be resumed through achieving the commercial purpose of the fishery. Hence, the results of our study are not only useful for assessing the biological, social, and economic impacts of existing resources but also provide appropriate measures to maintain long-term sustainability of the resource.In practice, the biological economic system can be stable by formulating some strategies according to the purification capacity and harvested efforts of fishers, such as taxation, license fees, lease of properties rights, seasonal harvesting, marine protected areas, and issuing purification subsidy. Therefore, not only ecological resources can be protected to promote silting reclamation, but also economics can be developed sustainably. The study has produced far-reaching influence to maintain the sustainable and favorable cycle of the ecosystem and improve the environment of water and soil and gain economic interest in the long term. Therefore, the theory proposed in this paper is significantly different from the existing models in the field of research which use differential algebraic equations and needs separate investigation.The entire study of the paper is primarily based on a deterministic framework. On the other hand, it will be more realistic if it is possible to incorporate ecological fluctuations and uncertainty in the model system. Thus, a future research problem would be to consider the system with interval parameter functions.



---

*Source: 548384-2014-11-16.xml*