Next Article in Journal
Ulam Stabilities and Instabilities of Euler–Lagrange-Rassias Quadratic Functional Equation in Non-Archimedean IFN Spaces
Next Article in Special Issue
Mathematical Model of Pest Control Using Different Release Rates of Sterile Insects and Natural Enemies
Previous Article in Journal
A Deterministic Learning Algorithm Estimating the Q-Matrix for Cognitive Diagnosis Models
Previous Article in Special Issue
“Realistic Choice of Annual Matrices Contracts the Range of λS Estimates” under Reproductive Uncertainty Too
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Properties of Solutions to Delay Differential Equations Describing Plankton—Fish Interaction

by
Maria A. Skvortsova
Laboratory of Differential and Difference Equations, Sobolev Institute of Mathematics, 4, Acad. Koptyug Avenue, 630090 Novosibirsk, Russia
Mathematics 2021, 9(23), 3064; https://doi.org/10.3390/math9233064
Submission received: 19 October 2021 / Revised: 18 November 2021 / Accepted: 24 November 2021 / Published: 28 November 2021
(This article belongs to the Special Issue Advances in the Mathematics of Ecological Modelling)

Abstract

:
We consider a system of differential equations with two delays describing plankton–fish interaction. We analyze the case when the equilibrium point of this system corresponding to the presence of only phytoplankton and the absence of zooplankton and fish is asymptotically stable. In this case, the asymptotic behavior of solutions to the system is studied. We establish estimates of solutions characterizing the stabilization rate at infinity to the considered equilibrium point. The results are obtained using Lyapunov–Krasovskii functionals.

1. Introduction

At present, there exist a large number of works devoted to the study of biological models described by delay differential equations (see, for example, the monographs [1,2,3,4] and bibliography therein). Among studied models, the models of population dynamics are widespread, in particular, models of predator–prey type, based on the classical predator–prey model that was introduced independently by A.J. Lotka [5] and V. Volterra [6]. An overview of some results for predator–prey models with delay is contained, for example, in [7,8].
In particular, predator–prey models are used when describing plankton–fish interaction (see, for example, [9,10,11,12,13,14]). Taking into account a model proposed in [14], in the present paper, we study the model of the following form:
d d t x ( t ) = r x ( t ) 1 x ( t ) K c 1 x ( t ) y ( t ) , d d t y ( t ) = d 1 y ( t ) + e 1 c 1 x ( t τ 1 ) y ( t τ 1 ) c 2 y ( t ) z ( t ) , d d t z ( t ) = d 2 z ( t ) + e 2 c 2 y ( t τ 2 ) z ( t τ 2 ) ,
which can be also considered a model of plankton–fish interaction (note that, in [14], it was considered the case that τ 2 = 0 , but nonlinear terms had a more general form). In this system, x ( t ) is the amount of phytoplankton, y ( t ) is the amount of zooplankton, and z ( t ) is the number of fish. It is assumed that phytoplankton is the favorite food of zooplankton, which serves as the favorite food of fish. The delay parameter τ 1 0 is responsible for time required for the appearance of new zooplankton, and the delay parameter τ 2 0 is responsible for time of fish maturation. The coefficients of the system have the following meaning: r > 0 is intrinsic growth rate of phytoplankton, K > 0 is environmental carrying capacity of phytoplankton, d 1 > 0 is mortality rate of zooplankton, d 2 > 0 is mortality rate of fish, c 1 0 is predation rate of zooplankton, c 2 0 is predation rate of fish, e 1 = b 1 e c 1 τ 1 , b 1 0 is birth rate of zooplankton, e 2 = b 2 e c 2 τ 2 , b 2 0 is birth rate of fish.
We consider system (1) for t > 0 , assuming that the initial conditions are given on the segment θ [ τ max , 0 ] , τ max = max { τ 1 , τ 2 } :
x ( θ ) = φ ( θ ) 0 , θ [ τ 1 , 0 ] , x ( + 0 ) = φ ( 0 ) > 0 , y ( θ ) = ψ ( θ ) 0 , θ [ τ max , 0 ] , y ( + 0 ) = ψ ( 0 ) , z ( θ ) = η ( θ ) 0 , θ [ τ 2 , 0 ] , z ( + 0 ) = η ( 0 ) ,
where φ ( θ ) , ψ ( θ ) , η ( θ ) are continuous functions. It is well known that a solution to the initial value problems (1) and (2) exists and is unique. Moreover, by analogy with [14], it is not difficult to show that the solution is defined on the entire right half-axis { t > 0 } , has non-negative components, and each component of the solution is a bounded function. In other words, there exist constants M 1 , M 2 , M 3 > 0 such that, for all t > 0 , the inequalities are valid
0 x ( t ) M 1 , 0 y ( t ) M 2 , 0 z ( t ) M 3 ,
i.e., the amount of plankton and fish cannot increase indefinitely.
One of the stationary solutions to system (1) is the equilibrium point ( x ( t ) , y ( t ) , z ( t ) ) T = ( K , 0 , 0 ) T , corresponding to the presence of only phytoplankton in the system and the absence of zooplankton and fish. We say that equilibrium point ( K , 0 , 0 ) T is asymptotically stable, if
1 . ε > 0 δ > 0 : max θ [ τ 1 , 0 ] | φ ( θ ) K | + max θ [ τ max , 0 ] | ψ ( θ ) | + max θ [ τ 2 , 0 ] | η ( θ ) | < δ
| x ( t ) K | + | y ( t ) | + | z ( t ) | < ε t > 0 ;
2 . ρ > 0 : max θ [ τ 1 , 0 ] | φ ( θ ) K | + max θ [ τ max , 0 ] | ψ ( θ ) | + max θ [ τ 2 , 0 ] | η ( θ ) | < ρ
| x ( t ) K | + | y ( t ) | + | z ( t ) | 0 a s t + .
In [14], it was noted that, for τ 2 = 0 , a sufficient condition for the asymptotic stability of this equilibrium point is the condition
e 1 c 1 K < d 1 ,
meaning a sufficiently high mortality of zooplankton. By analogy with [14], it is not difficult to establish that, for τ 2 > 0 , condition (3) is also a sufficient condition for the asymptotic stability of the considered equilibrium point. Moreover, it can be shown that all solutions to system (1) with initial conditions of the form (2) are stabilized at infinity to this equilibrium point.
The aim of the present paper is, under condition (3), to establish estimates for all components of the solution to the initial value problems (1) and (2) characterizing the stabilization rate at infinity to equilibrium point ( K , 0 , 0 ) T . To achieve our goal, we will use Lyapunov–Krasovskii functionals. Note that such functionals are actively used to obtain estimates of solutions to various classes of systems with delay (see, for example, [15,16,17,18], where systems of differential equations of delayed type were considered, and [19,20,21,22,23,24,25,26,27,28,29], where systems of differential equations of neutral type were considered). For specific biological models, the results based on the use of Lyapunov–Krasovskii functionals are contained, for example, in [30,31,32,33,34,35].

2. Main Results

As already noted above, throughout the paper, we assume that condition (3) is fulfilled, which guarantees the asymptotic stability of equilibrium point ( K , 0 , 0 ) T of system (1). We obtain estimates for all components of the solution to the initial value problems (1) and (2).
First, we formulate an auxiliary statement.
Lemma 1.
Let ( x ( t ) , y ( t ) , z ( t ) ) T be the solution to the initial value problems (1) and (2). Then, for the first component of the solution x ( t ) , the estimate holds
0 < x ( t ) K + ( φ ( 0 ) K ) e r t , t > 0 .
Proof. 
Due to the initial condition (2), we have x ( + 0 ) = φ ( 0 ) > 0 . Hence, from the first equation of system (1), it follows that x ( t ) > 0 for all t > 0 . Further, since y ( t ) 0 for all t > 0 , from this equation, we obtain the inequality
d d t x ( t ) r x ( t ) 1 x ( t ) K .
Changing the variables
x ¯ ( t ) = x ( t ) K ,
we establish the following estimate:
d d t x ¯ ( t ) r K x ¯ ( t ) x ¯ ( t ) + K r x ¯ ( t ) .
From here, it is not difficult to obtain the inequality
x ¯ ( t ) x ¯ ( 0 ) e r t .
Taking into account that x ¯ ( t ) = x ( t ) K , we establish (4).
Lemma is proved. □
Now, we establish estimates for the second component of the solution y ( t ) .
Lemma 2.
Let ( x ( t ) , y ( t ) , z ( t ) ) T be the solution to the initial value problems (1) and (2). Then, for the second component of the solution y ( t ) , the estimate holds
0 y ( t ) α , α = ψ ( 0 ) + e 1 c 1 τ 1 0 φ ( s ) ψ ( s ) d s , t [ 0 , τ 1 ] .
Proof. 
Let t [ 0 , τ 1 ] . Taking into account the initial conditions (2), and considering y ( t ) 0 and z ( t ) 0 , from the second equation of system (1), it is not difficult to obtain the estimate
d d t y ( t ) d 1 y ( t ) + e 1 c 1 φ ( t τ 1 ) ψ ( t τ 1 ) .
This estimate can be rewritten in an equivalent form:
d d t e d 1 t y ( t ) e 1 c 1 e d 1 t φ ( t τ 1 ) ψ ( t τ 1 ) ,
from which follows the inequality
y ( t ) ψ ( 0 ) e d 1 t + e 1 c 1 0 t e d 1 ( t s ) φ ( s τ 1 ) ψ ( s τ 1 ) d s
ψ ( 0 ) + e 1 c 1 0 τ 1 φ ( s τ 1 ) ψ ( s τ 1 ) d s ,
which coincides with (5).
Lemma is proved. □
Let us proceed to obtaining estimates for the second component of the solution y ( t ) for t > τ 1 . To do this, consider the Lyapunov–Krasovskii functional of the following form
V ( t , y ) = y 2 ( t ) + t τ 1 t e k 1 ( t s ) m 1 ( s ) y 2 ( s ) d s ,
where k 1 > 0 is such a number that the inequality is satisfied
e 1 c 1 K e k 1 τ 1 / 2 < d 1
(by virtue of condition (3), such k 1 > 0 exists),
m 1 ( t ) = e 1 c 1 e k 1 τ 1 / 2 K + ( φ ( 0 ) K ) e r t .
The following theorem is valid.
Theorem 1.
Let condition (3) be satisfied, and let ( x ( t ) , y ( t ) , z ( t ) ) T be the solution to the initial value problems (1) and (2). Then, for the second component of the solution y ( t ) , the estimate holds
0 y ( t ) V ( τ 1 , α ) e J 1 / 2 e ε ( t τ 1 ) / 2 , t > τ 1 ,
where α is defined in (5),
J 1 = 0 , i f 0 < φ ( 0 ) K , e 1 c 1 e k 1 τ 1 / 2 ( φ ( 0 ) K ) ( 1 + e r τ 1 ) r , i f φ ( 0 ) > K ,
ε = min 2 d 1 e 1 c 1 K e k 1 τ 1 / 2 , k 1 > 0 .
Proof. 
Let t > τ 1 . Using inequality (4), and considering y ( t ) 0 and z ( t ) 0 , from the second equation of system (1), we obtain the estimate
d d t y ( t ) d 1 y ( t ) + e 1 c 1 K + ( φ ( 0 ) K ) e r ( t τ 1 ) y ( t τ 1 ) .
Now we consider Lyapunov–Krasovskii functional (6). Differentiating it along the solution to the initial value problems (1) and (2), we establish validity of the following inequality
d d t V ( t , y ) = 2 y ( t ) d d t y ( t ) + m 1 ( t ) y 2 ( t ) e k 1 τ 1 m 1 ( t τ 1 ) y 2 ( t τ 1 )
k 1 t τ 1 t e k 1 ( t s ) m 1 ( s ) y 2 ( s ) d s
( 2 d 1 m 1 ( t ) ) y 2 ( t ) + 2 e 1 c 1 K + ( φ ( 0 ) K ) e r ( t τ 1 ) y ( t ) y ( t τ 1 )
e k 1 τ 1 m 1 ( t τ 1 ) y 2 ( t τ 1 ) k 1 t τ 1 t e k 1 ( t s ) m 1 ( s ) y 2 ( s ) d s .
Hence, using the inequality
2 e 1 c 1 K + ( φ ( 0 ) K ) e r ( t τ 1 ) y ( t ) y ( t τ 1 ) e k 1 τ 1 m 1 ( t τ 1 ) y 2 ( t τ 1 )
e 1 2 c 1 2 K + ( φ ( 0 ) K ) e r ( t τ 1 ) 2 e k 1 τ 1 m 1 ( t τ 1 ) y 2 ( t )
and considering definition (7) of function m 1 ( t ) , we obtain the estimate
d d t V ( t , y ) 2 d 1 e 1 c 1 K e k 1 τ 1 / 2 y 2 ( t ) + e 1 c 1 e k 1 τ 1 / 2 ( φ ( 0 ) K ) ( 1 + e r τ 1 ) e r t y 2 ( t )
k 1 t τ 1 t e k 1 ( t s ) m 1 ( s ) y 2 ( s ) d s .
By virtue of definition (10) of ε , from here, the inequality follows:
d d t V ( t , y ) ε V ( t , y ) + e 1 c 1 e k 1 τ 1 / 2 ( φ ( 0 ) K ) ( 1 + e r τ 1 ) e r ( t τ 1 ) y 2 ( t ) .
First, let 0 < φ ( 0 ) K ; then,
d d t V ( t , y ) ε V ( t , y ) ;
therefore,
V ( t , y ) V ( τ 1 , y ) e ε ( t τ 1 ) .
From definition (6) and inequality (5), the estimates follow: y 2 ( t ) V ( t , y ) , V ( τ 1 , y ) V ( τ 1 , α ) , from which we establish (8).
Now, let φ ( 0 ) > K . Then, from inequality (11), we obtain the estimate
d d t V ( t , y ) ε + J 1 r e r ( t τ 1 ) V ( t , y ) ,
where J 1 is defined in (9). Hence, the inequality follows
V ( t , y ) V ( τ 1 , y ) exp τ 1 t ε + J 1 r e r ( s τ 1 ) d s V ( τ 1 , y ) e J 1 e ε ( t τ 1 ) .
This estimate directly implies (8).
Theorem is proved. □
Now, we obtain estimates for the third component of the solution z ( t ) . First, we consider the case t [ 0 , τ 2 ] .
Lemma 3.
Let ( x ( t ) , y ( t ) , z ( t ) ) T be the solution to the initial value problems (1) and (2). Then, for the third component of the solution z ( t ) , the estimate holds
0 z ( t ) β , β = η ( 0 ) + e 2 c 2 τ 2 0 ψ ( s ) η ( s ) d s , t [ 0 , τ 2 ] .
Proof. 
Let t [ 0 , τ 2 ] . Taking into account the initial conditions (2), from the third equation of system (1), it is not difficult to obtain the representation
z ( t ) = η ( 0 ) e d 2 t + e 2 c 2 0 t e d 2 ( t s ) ψ ( s τ 2 ) η ( s τ 2 ) d s ;
hence, estimate (12) follows.
Lemma is proved. □
Now, suppose that t [ τ 2 , τ 1 + τ 2 ] . Consider the Lyapunov–Krasovskii functional
U ( t , z ) = z 2 ( t ) + e 2 c 2 α t τ 2 t z 2 ( s ) d s ,
where α is defined in (5).
The following statement takes place.
Lemma 4.
Let ( x ( t ) , y ( t ) , z ( t ) ) T be the solution to the initial value problems (1) and (2). Then, for the third component of the solution z ( t ) , the estimate holds
0 z ( t ) γ , γ = U ( τ 2 , β ) e ω τ 1 / 2 , t [ τ 2 , τ 1 + τ 2 ] ,
where β is defined in (12),
ω = max { 2 ( e 2 c 2 α d 2 ) , 0 } .
Proof. 
Let t [ τ 2 , τ 1 + τ 2 ] . Using inequality (5), from the third equation of system (1), we obtain the estimate
d d t z ( t ) d 2 z ( t ) + e 2 c 2 α z ( t τ 2 ) .
Now, we consider the Lyapunov–Krasovskii functional (13). Differentiating it along the solution to the initial value problems (1) and (2), it is not difficult to establish the following inequality:
d d t U ( t , z ) = 2 z ( t ) d d t z ( t ) + e 2 c 2 α z 2 ( t ) e 2 c 2 α z 2 ( t τ 1 )
( e 2 c 2 α 2 d 2 ) z 2 ( t ) + 2 e 2 c 2 α z ( t ) z ( t τ 2 ) e 2 c 2 α z 2 ( t τ 1 )
2 ( e 2 c 2 α d 2 ) z 2 ( t ) ω U ( t , z ) .
From this inequality, we obtain the estimate
U ( t , z ) U ( τ 2 , z ) e ω ( t τ 2 ) U ( τ 2 , z ) e ω τ 1 .
From definition (13) and inequality (12), the estimates follow: z 2 ( t ) U ( t , z ) , U ( τ 2 , z ) U ( τ 2 , β ) , from which we establish (14).
Lemma is proved. □
To obtain estimates for the third component of the solution z ( t ) for t > τ 1 + τ 2 , we consider the Lyapunov–Krasovskii functional of the following form:
W ( t , z ) = z 2 ( t ) + t τ 2 t e k 2 ( t s ) m 2 ( s ) z 2 ( s ) d s ,
where k 2 > 0 is arbitrary,
m 2 ( t ) = e 2 c 2 e k 2 τ 2 / 2 V ( τ 1 , α ) e J 1 / 2 e ε ( t τ 1 ) / 2 ,
V ( t , y ) is defined in (6), α is defined in (5), J 1 is defined in (9), and ε is defined in (10).
The following theorem is valid.
Theorem 2
Let condition (3) be satisfied, and let ( x ( t ) , y ( t ) , z ( t ) ) T be the solution to the initial value problems (1) and (2). Then, for the third component of the solution z ( t ) , the estimate holds
0 z ( t ) W ( τ 1 + τ 2 , γ ) e J 2 / 2 e σ ( t τ 1 τ 2 ) / 2 , t > τ 1 + τ 2 ,
where γ is defined in (14),
J 2 = e 2 c 2 e k 2 τ 2 / 2 V ( τ 1 , α ) e J 1 / 2 ( 1 + e ε τ 2 / 2 ) ( ε / 2 ) ,
σ = min { 2 d 2 , k 2 } .
Proof. 
Let t > τ 1 + τ 2 . Using inequality (8), from the third equation of system (1), we obtain the estimate
d d t z ( t ) d 2 z ( t ) + e 2 c 2 V ( τ 1 , α ) e J 1 / 2 e ε ( t τ 1 τ 2 ) / 2 z ( t τ 2 ) .
Now, we consider the Lyapunov–Krasovskii functional (15). Differentiating it along the solution to the initial value problems (1) and (2), we establish the validity of the following inequality:
d d t W ( t , z ) = 2 z ( t ) d d t z ( t ) + m 2 ( t ) z 2 ( t ) e k 2 τ 2 m 2 ( t τ 2 ) z 2 ( t τ 2 )
k 2 t τ 2 t e k 2 ( t s ) m 2 ( s ) z 2 ( s ) d s
( 2 d 2 m 2 ( t ) ) z 2 ( t ) + 2 e 2 c 2 V ( τ 1 , α ) e J 1 / 2 e ε ( t τ 1 τ 2 ) / 2 z ( t ) z ( t τ 2 )
e k 2 τ 2 m 2 ( t τ 2 ) z 2 ( t τ 2 ) k 2 t τ 2 t e k 2 ( t s ) m 2 ( s ) z 2 ( s ) d s .
Hence, using the inequality
2 e 2 c 2 V ( τ 1 , α ) e J 1 / 2 e ε ( t τ 1 τ 2 ) / 2 z ( t ) z ( t τ 2 ) e k 2 τ 2 m 2 ( t τ 2 ) z 2 ( t τ 2 )
e 2 c 2 V ( τ 1 , α ) e J 1 / 2 e ε ( t τ 1 τ 2 ) / 2 2 e k 2 τ 2 m 2 ( t τ 2 ) z 2 ( t )
and considering definition (16) of function m 2 ( t ) , we obtain the estimate
d d t W ( t , z ) 2 d 2 z 2 ( t ) + e 2 c 2 e k 2 τ 2 / 2 V ( τ 1 , α ) e J 1 / 2 ( 1 + e ε τ 2 / 2 ) e ε ( t τ 1 τ 2 ) / 2 z 2 ( t )
k 2 t τ 2 t e k 2 ( t s ) m 2 ( s ) z 2 ( s ) d s .
By virtue of definition (19) of σ and definition (18) of J 2 , the inequality follows:
d d t W ( t , z ) σ + J 2 ε 2 e ε ( t τ 1 τ 2 ) / 2 W ( t , z ) .
Therefore,
W ( t , z ) W ( τ 1 + τ 2 , z ) exp τ 1 + τ 2 t σ + J 2 ε 2 e ε ( s τ 1 τ 2 ) / 2 d s
W ( τ 1 + τ 2 , z ) e J 2 e σ ( t τ 1 τ 2 ) .
Since z 2 ( t ) W ( t , z ) and W ( τ 1 + τ 2 , z ) W ( τ 1 + τ 2 , γ ) , where γ is defined in (14), this estimate directly implies (17).
Theorem is proved. □
Finally, we obtain estimates for the first component of the solution x ( t ) .
The following theorem is valid.
Theorem 3.
Let condition (3) be satisfied, and let ( x ( t ) , y ( t ) , z ( t ) ) T be the solution to the initial value problems (1) and (2). Then, for the first component of the solution x ( t ) , the estimate holds
K x ( t ) K e c 1 Y min K , φ ( 0 ) X 1 e r t + X 2 ( t τ 1 ) e min { r , ( ε / 2 ) } ( t τ 1 ) , t > τ 1 ,
where
Y = α τ 1 + 2 ε V ( τ 1 , α ) e J 1 / 2 ,
X 1 = max { K φ ( 0 ) , 0 } + φ ( 0 ) c 1 α ( e r τ 1 1 ) r ,
X 2 = φ ( 0 ) c 1 V ( τ 1 , α ) e J 1 / 2 ,
α is defined in (5), V ( t , y ) is defined in (6), J 1 is defined in (9), and ε is defined in (10).
Proof. 
Let t > τ 1 . From the first equation of system (1), it is not difficult to obtain the representation
K x ( t ) = K K φ ( 0 ) + φ ( 0 ) c 1 0 t y ( s ) e r s exp c 1 0 s y ( ξ ) d ξ d s
× K + φ ( 0 ) r 0 t e r s exp c 1 0 s y ( ξ ) d ξ d s 1 .
Hence, the inequality follows:
K x ( t ) K max { K φ ( 0 ) , 0 } + φ ( 0 ) c 1 0 t y ( s ) e r s d s
× K + φ ( 0 ) r 0 t e r s exp c 1 0 s y ( ξ ) d ξ d s 1 .
By virtue of the estimate
K + φ ( 0 ) r 0 t e r s exp c 1 0 s y ( ξ ) d ξ d s
K + φ ( 0 ) r 0 t e r s d s exp c 1 0 t y ( ξ ) d ξ
min K , φ ( 0 ) e r t exp c 1 0 t y ( ξ ) d ξ ,
we obtain the inequality
K x ( t ) K min K , φ ( 0 ) exp c 1 0 t y ( ξ ) d ξ
× max { K φ ( 0 ) , 0 } e r t + φ ( 0 ) c 1 0 t y ( s ) e r ( t s ) d s .
Now, we estimate the integrals. By virtue of the inequalities (5) and (8), we have
0 t y ( ξ ) d ξ 0 τ 1 α d ξ + τ 1 V ( τ 1 , α ) e J 1 / 2 e ε ( ξ τ 1 ) / 2 d ξ = α τ 1 + 2 ε V ( τ 1 , α ) e J 1 / 2 = Y ;
hence,
K x ( t ) K e c 1 Y min K , φ ( 0 ) max { K φ ( 0 ) , 0 } e r t + φ ( 0 ) c 1 0 t y ( s ) e r ( t s ) d s .
Next,
0 t y ( s ) e r ( t s ) d s 0 τ 1 α e r ( t s ) d s + τ 1 t V ( τ 1 , α ) e J 1 / 2 e ε ( s τ 1 ) / 2 e r ( t s ) d s
α ( e r τ 1 1 ) r e r t + V ( τ 1 , α ) e J 1 / 2 ( t τ 1 ) e min { r , ( ε / 2 ) } ( t τ 1 ) .
By virtue of this inequality, taking into account notations (22) and (23), from estimate (24), we obtain inequality (20).
Theorem is proved. □
From Theorem 3 and Lemma 1, the statement follows.
Corollary 1.
Let condition (3) be satisfied, and let ( x ( t ) , y ( t ) , z ( t ) ) T be the solution to the initial value problems (1) and (2). Then, for t > τ 1 , for the first component of the solution x ( t ) , the estimate holds
K K e c 1 Y min K , φ ( 0 ) X 1 e r t + X 2 ( t τ 1 ) e min { r , ( ε / 2 ) } ( t τ 1 ) x ( t ) K + ( φ ( 0 ) K ) e r t ,
where Y, X 1 , and X 2 are defined in (21)–(23), and ε is defined in (10).

3. Conclusions

In the present paper, we have considered a system of differential equations with two delays describing the interaction between fish, zooplankton, and phytoplankton. Provided that zooplankton mortality was sufficiently high, we have established estimates of solutions that characterize the decay rates of the amount of zooplankton and fish, as well as the change rate to a positive constant value of the amount of phytoplankton. The obtained estimates are constructive, and all the values responsible for the stabilization rates are indicated explicitly. To obtain the results, special Lyapunov–Krasovskii functionals were constructed.
Note that the results of this paper have been obtained under the assumption that the equilibrium point corresponding to the presence of only phytoplankton in the system is asymptotically stable. Under certain conditions on the coefficients of the system, there also exist other equilibrium points corresponding to the presence of all populations in the system. Finding the conditions of asymptotic stability of these equilibrium points, as well as obtaining estimates of the stabilization rate of solutions and estimates for attraction sets, is of interest both from a mathematical and biological points of view. The construction of special Lyapunov–Krasovskii functionals to study the asymptotic behavior of solutions in the case of asymptotic stability of equilibrium points corresponding to the presence of all populations in the system is the goal of further research.

Funding

The author was supported by the Mathematical Center in Akademgorodok (Agreement No. 075-15-2019-1613 with the Ministry of Science and Higher Education of the Russian Federation).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author is grateful to G.V. Demidenko and I.I. Matveeva for their attention to the research.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Gopalsamy, K. Stability and oscillations in delay differential equations of population dynamics. In Mathematics and its Applications (Dordrecht); Kluwer Academic: Dordrecht, The Netherlands, 1992; Volume 74. [Google Scholar]
  2. Kuang, Y. Delay differential equations: With applications in population dynamics. In Mathematics in Science and Engineering; Academic Press: Boston, MA, USA, 1993; Volume 191. [Google Scholar]
  3. Smith, H.L. Monotone dynamical systems: An introduction to the theory of competitive and cooperative systems. In Mathematical Surveys and Monographs; American Mathematical Society (AMS): Providence, RI, USA, 1995; Volume 41. [Google Scholar]
  4. Erneux, T. Applied delay differential equations. In Surveys and Tutorials in the Applied Mathematical Sciences; Springer: New York, NY, USA, 2009; Volume 3. [Google Scholar]
  5. Lotka, A.J. The Elements of Physical Biology; Williams&Wilkins Co.: Baltimore, MA, USA; Tindall&Cox: London, UK, 1925. [Google Scholar]
  6. Volterra, V. Variazioni e fluttuazioni del numero d’individui in specie animali conviventi. Mem. Della Accad. Naz. Dei Lincei Cl. Sci. Fis. Mat. Nat. 1927, 2, 31–113. [Google Scholar]
  7. Liu, S.; Chen, L.; Agarwal, R. Recent progress on stage-structured population dynamics. Math. Comput. Model. 2002, 36, 1319–1360. [Google Scholar] [CrossRef]
  8. Ruan, S. On nonlinear dynamics of predator-prey models with discrete delay. Math. Model. Nat. Phenom. 2009, 4, 140–188. [Google Scholar] [CrossRef]
  9. Pal, S.; Chatterjee, A. Dynamics of the interaction of plankton and planktivorous fish with delay. Cogent. Math. 2015, 2, 1074337. [Google Scholar] [CrossRef]
  10. Kloosterman, M.; Campbell, S.A.; Poulin, F.J. An NPZ model with state-dependent delay due to size-structure in juvenile zooplankton. SIAM J. Appl. Math. 2016, 76, 551–577. [Google Scholar] [CrossRef] [Green Version]
  11. Meng, X.Y.; Wu, Y.Q. Bifurcation and control in a singular phytoplankton–zooplankton–fish model with nonlinear fish harvesting and taxation. Int. J. Bifurc. Chaos Appl. Sci. Eng. 2018, 28, 1850042. [Google Scholar] [CrossRef]
  12. Zheng, W.; Sugie, J. Global asymptotic stability and equiasymptotic stability for a time-varying phytoplankton–zooplankton–fish system. Nonlinear Anal. Real World Appl. 2019, 46, 116–136. [Google Scholar] [CrossRef]
  13. Raw, S.N.; Tiwari, B.; Mishra, P. Analysis of a plankton–fish model with external toxicity and nonlinear harvesting. Ric. Mat. 2020, 69, 653–681. [Google Scholar] [CrossRef]
  14. Thakur, N.K.; Ojha, A. Complex dynamics of delay-induced plankton–fish interaction exhibiting defense. SN Appl. Sci. 2020, 2, 1114. [Google Scholar] [CrossRef]
  15. Demidenko, G.V.; Matveeva, I.I. Asymptotic properties of solutions to delay differential equations. Vestn. Novosib. Univ. Ser. Mat. Mekh. Inform. 2005, 5, 20–28. [Google Scholar]
  16. Khusainov, D.Y.; Ivanov, A.F.; Kozhametov, A.T. Convergence estimates for solutions of linear stationary systems of differential-difference equations with constant delay. Differ. Equ. 2005, 41, 1196–1200. [Google Scholar] [CrossRef]
  17. Mondié, S.; Kharitonov, V.L. Exponential estimates for retarded time-delay systems: LMI approach. IEEE Trans. Autom. Control 2005, 50, 268–273. [Google Scholar] [CrossRef]
  18. Demidenko, G.V.; Matveeva, I.I. Stability of solutions to delay differential equations with periodic coefficients of linear terms. Sib. Math. J. 2007, 48, 824–836. [Google Scholar] [CrossRef]
  19. Demidenko, G.V. Stability of solutions to linear differential equations of neutral type. J. Anal. Appl. 2009, 7, 119–130. [Google Scholar]
  20. Demidenko, G.V.; Matveeva, I.I. On estimates of solutions to systems of differential equations of neutral type with periodic coefficients. Sib. Math. J. 2014, 55, 866–881. [Google Scholar] [CrossRef]
  21. Demidenko, G.V.; Matveeva, I.I. Estimates for solutions to a class of nonlinear time-delay systems of neutral type. Electron. J. Differ. Equ. 2015, 2015, 34. [Google Scholar]
  22. Demidenko, G.V.; Matveeva, I.I. Estimates for solutions to a class of time-delay systems of neutral type with periodic coefficients and several delays. Electron. J. Qual. Theory Differ. Equ. 2015, 2015, 83. [Google Scholar] [CrossRef]
  23. Matveeva, I.I. On exponential stability of solutions to periodic neutral-type systems. Sib. Math. J. 2017, 58, 264–270. [Google Scholar] [CrossRef]
  24. Matveeva, I.I. On the exponential stability of solutions of periodic systems of the neutral type with several delays. Differ. Equ. 2017, 53, 725–735. [Google Scholar] [CrossRef]
  25. Demidenko, G.V.; Matveeva, I.I.; Skvortsova, M.A. Estimates for solutions to neutral differential equations with periodic coefficients of linear terms. Sib. Math. J. 2019, 60, 828–841. [Google Scholar] [CrossRef]
  26. Matveeva, I.I. Estimates for exponential decay of solutions to one class of nonlinear systems of neutral type with periodic coefficients. Comput. Math. Math. Phys. 2020, 60, 601–609. [Google Scholar] [CrossRef]
  27. Matveeva, I.I. Exponential stability of solutions to nonlinear time-varying delay systems of neutral type equations with periodic coefficients. Electron. J. Differ. Equ. 2020, 2020, 20. [Google Scholar]
  28. Yskak, T. Estimates for solutions of one class of systems of equations of neutral type with distributed delay. Sib. Electron. Math. Rep. 2020, 17, 416–427. [Google Scholar] [CrossRef]
  29. Matveeva, I.I. Estimates for solutions to a class of nonautonomous systems of neutral type with unbounded delay. Sib. Math. J. 2021, 62, 468–481. [Google Scholar] [CrossRef]
  30. Skvortsova, M.A. Asymptotic properties of solutions to a system describing the spread of avian influenza. Sib. Electron. Math. Rep. 2016, 13, 782–798. [Google Scholar]
  31. Skvortsova, M.A. Estimates for solutions in a predator–prey model with delay. Izv. Irkutsk. Gos. Univ. Ser. Mat. 2018, 25, 109–125. [Google Scholar]
  32. Skvortsova, M.A. Asymptotic properties of solutions in a model of antibacterial immune response. Sib. Electron. Math. Rep. 2018, 15, 1198–1215. [Google Scholar]
  33. Skvortsova, M.A. On estimates of solutions in a predator–prey model with two delays. Sib. Electron. Math. Rep. 2018, 15, 1697–1718. [Google Scholar] [CrossRef] [Green Version]
  34. Skvortsova, M.A. Asymptotic properties of solutions in a model of interaction of populations with several delays. Math. Notes NEFU 2019, 26, 63–72. [Google Scholar]
  35. Skvortsova, M.A.; Yskak, T. Asymptotic behavior of solutions in one predator-prey model with delay. Sib. Math. J. 2021, 62, 324–336. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Skvortsova, M.A. Asymptotic Properties of Solutions to Delay Differential Equations Describing Plankton—Fish Interaction. Mathematics 2021, 9, 3064. https://doi.org/10.3390/math9233064

AMA Style

Skvortsova MA. Asymptotic Properties of Solutions to Delay Differential Equations Describing Plankton—Fish Interaction. Mathematics. 2021; 9(23):3064. https://doi.org/10.3390/math9233064

Chicago/Turabian Style

Skvortsova, Maria A. 2021. "Asymptotic Properties of Solutions to Delay Differential Equations Describing Plankton—Fish Interaction" Mathematics 9, no. 23: 3064. https://doi.org/10.3390/math9233064

APA Style

Skvortsova, M. A. (2021). Asymptotic Properties of Solutions to Delay Differential Equations Describing Plankton—Fish Interaction. Mathematics, 9(23), 3064. https://doi.org/10.3390/math9233064

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop